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level  density  of  states  and  <1  >  is  the  Fermi  level  mean  square 
electron-ion  matrix  element.  The  parameter  n  is  calculated  for 
the  f-band  metals  La,  Ce  and  Th  using  self-consistent  muffin-tin 
potentials.  The  corrections  within  the  rigid  muffin-tin  ap¬ 
proximation  (RMTA)  due  to  spin-orbit  coupling  derived  by  John 
and  Hamann  are  found  to  increase  n  by  1%  in  La  and  in  Th,  and 
by  10%  and  33%  in  the  a  and  y  phases  of  Ce,  respectively.  The 
RMTA  expression  is  reinterpreted  to  provide  an  understanding  of 
the  low  values  of  n  in  Ce  (0.7  to  l.leV/S2)  in  spite  of  an 
extremely  large  density  of  states  at  the  Fermi  level.  Using 
moments  of  the  measured  phonon  spectra,  the  superconducting 
transition  temperature  is  calculated  and  compared  with  experi¬ 
ment.  The  RMTA  appears  to  overestimate  n  by  approximately  253 
in  these  metals,  although  spin  fluctuations  may  be  contributing 
to  this  apparent  discrepancy  with  experiment. 
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The  expression  for  the  McMillan-Hopfield  constant  r'=N(Ep) 
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<1  >  within  the  rigid-muf in-tin  approximation  (RMTA)  is  rein¬ 
terpreted  in  terms  of  single-channel  electron-ion  matrix  ele- 
2 

ments  I  ^  and  angular-momentum-character  fractions  f  ^  of 

the  electronic  states  at  the  Fermi  level.  Reexamination  of  Nb- 
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and  V-based  A15  compounds  suggests  that  I  +1  is  more  nearly 
an  atomic  property,  and  thus  transferable  from  system  to  system, 
than  other  commonly  used  quantitities .  The  fractions  f  are 
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increasing  n  within  RMTA  are  discussed. 
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A  microscopic  treatment  of  the  consequences  for  supercon¬ 
ductivity  of  a  nonconstant  electronic  density  of  states  is 
presented.  Generalized  Ellashberg  gap  equations  valid  for  a 
varying  density  of  states  are  presented,  from  which  the  change 
of  T  with  static  or  thermal  disorder  can  be  calculated.  The 
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temperature  dependence  of  the  effective  mass  is  shown  to  be 
altered  by  disorder.  Use  of  these  results  provides  a  possible 


experimental  approach  for  deducing  the  energy  variation  of  the 
density  of  states  of  superconductors. 


4.  "INFLUENCE  OF  ELECTRONIC  STRUCTURE  ON  SUPERCONDUCTING  PROP¬ 
ERTIES  OF  COMPLEX  CRYSTALS:  THEORY  AND  APPLICATION  TO  NbjSn. " 
Solid  State  Commun.  38,  95  (1981).  With  B.M.  Klein. 

A  thermodynamic  theory  valid  for  complex  crystalline  super¬ 
conductors  is  applied  to  Nb^Sn.  It  is  shown  how  the  structure 
of  the  Eliashberg  equation  and  the  solution  is  altered  by  den¬ 
sity  of  states  fine  structure.  One  important  implication  is 
that  the  current  method  of  inversion  of  tunneling  data  can  be 
only  approximate. 

5.  "PARAMETER-FREE  CALCULATION  OF  THE  ENHANCED  SPIN  SUSCEPTI¬ 
BILITY  OF  Nb^Sn  INCLUDING  ELECTRON-PHONON  EFFECTS."  Physica 
107B,  703  (1981). 

A  parameter-free  calculation  of  the  3pin  susceptibility 

XC„(T)  of  Nb-.Sn,  which  includes  static  and  dynamic  disorder  as 
sp  3 

well  as  exchange-correlation  enhancement,  is  presented.  It  Is 
found  that  Xgp  Is  only  15%  of  the  measured  susceptibility,  and 
its  small  temperature  dependence  cannot  account  for  the  experi¬ 
mental  findings. 

6.  "GENERALIZATION  OF  THE  THEORY  OF  THE  ELECTRON-PHONON  INTER¬ 
ACTION:  THERMODYNAMIC  FORMULATION  OF  SUPERCONDUCTING-  AND 
NORMAL-STATE  PROPERTIES."  Phys.  Rev.  B26,  1186  (1982). 

A  thermodynamic  formulation  for  the  electron  self-energy 
is  given  which  is  applicable  when  the  electronic  spectrum 
posseses  structure  on  the  scale  of  phonon  frequencies,  provided 
only  that  the  ratio  of  phonon  phase  velocity  to  electron  Fermi 
velocity  is  small.  Electron-phonon,  Coulomb,  and  electron- 
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defect  Interactions  are  included  on  an  equal  footing  and  it  is 
shown  that  their  different  frequency  dependencies  lead  to 
specific  effects  on  the  Eliashberg  self -energy:  (a)  The  Coulomb 
interaction  contributes  nothing  of  essence  to  the  normal-state 
self-energy  (in  this  isotropic  approximation)  but  retains  its 
usual  depairing  effect  upon  the  superconducting  gap  function, 
(b)  defects  affect  superconducting  properties  primarily  through 
a  broadening  of  the  electronic  spectrum,  and  (c)  phonons  con¬ 
tribute  a  thermal  shift  and  broadening  as  well  as  the  mass 
enhancement.  A  generalization  to  intensive  electron-phonon, 
electron-electron,  and  electron-defect  interaction  constants  is 
necessary  to  redevelop  an  intuition  into  the  effects  of  these 
interactions.  The  change  in  the  structure  of  the  Eliashberg 
equation  due  to  a  nonconstant  density  of  states  (DOS)  and  the 
consequent  Interplay  between  static  and  thermal  disorder  is 
analyzed  In  detail,  with  a  central  feature  being  the  change  in 
frequency  dependence  of  the  self -energy  compared  to  a  constant 
DOS  solution.  The  effect  of  DOS  structure  on  the  superconduct¬ 
ing  transition  temperature  T  ,  which  is  manifested  In  the 

C 

defect  dependence  of  T  ,  Is  analyzed  in  detail.  Further  it  is 

V 

proposed  that  an  extension  of  the  self-consistent  Eliashberg 
approach  be  extended  above  Tq  to  determine  the  normal-state 
self-energy  and  thereby  the  electronic  contribution  to  thermo¬ 
dynamic  quantities.  Phonon  broadening  is  shown  to  affect 
the  spin  susceptibility  at  finite  temperature.  Reinterpretation 
of  several  of  the  anomalous  properties  of  A15  compounds  in  terms 
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of  the  present  theory  is  suggested.  Several  aspects  of  the 
theory  are  compared  to  experimental  data  for  Nb^Sn. 

7.  "RENORMALIZED  THERMAL  DISTRIBUTION  FUNCTION  IN  AN  INTERACT¬ 
ING  ELECTRON-PHONON  SYSTEM."  Phys.  Rev.  Lett.  48,  1548  (1982). 

The  electron-phonon  interaction  is  used  to  demonstrate  the 
important  effect  of  interactions  on  the  electronic  distribution 
function  at  finite  temperature.  It  is  shown  that  the  usual 
picture  of  "thermal  (Fermi)  smearing"  is  a  greatly  oversimpli¬ 
fied  one.  The  distribution  function  resulting  from  an  Einstein 
spectrum  with  various  coupling  strengths  is  presented  and 
interpreted,  and  an  exact  expression  for  the  spin  susceptibility 
is  used  to  illustrate  the  utility  of  this  novel  viewpoint  for 
thermodynamics. 

8.  "THEORY  OF  THE  NORMAL  STATE  HEAT  CAPACITY  OF  NbjSn.  " 
Superconductivity  in  d-  and  f-band  metals.  Edited  by  W.  Buchel 
and  W.  Weber  (Kernf orschungszent rum  Karlsruhe  GmbH,  1982), 
p. 97.  With  B.M.  Klein. 

2 

The  experimental  a  F  of  Wolf  et  al.  is  applied  to  calcu¬ 
late  the  T-dependence  of  the  electronic  heat  capacity.  It  Is 
found  that.  If  the  phonon  spectrum  Is  Debye-like  (which  is  not 
the  case  in  Nb^Sn),  this  T-dependence  results  in  an  extrap¬ 
olated  value  of  N(Ep)  (1  +  x)  which  is  overestimated  by  an 
amount  which  Is  proportional  N(Ep)x,  equal  to  13?  for  Nb^Sn. 
It  Is  further  shown  that  the  change  In  slope  of  C/T  found  by 
Stewart,  Cort  and  Webb  can  be  modeled  by  a  combined  Debye  and 
Einstein  spectrum  with  9^=267  ^  and  having  1.5?  of  the  acous¬ 
tic  modes  at  an  Einstein  frequency  Q=40  K. 
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The  electron-ion  scattering  parameter  rj  is  calculated  for  the  f-band  metals  La.  Ce  and  Th  using  self-consistent  muilin-tin 
potentials.  The  corrections  within  the  rigid  muffin-tin  approximation  (RMTA)  due  to  spin-orbit  coupling  derived  by  John 
and  Hamann  are  found  to  increase  n  by  1%  in  La  and  4%  in  Th,  and  by  10%  and  33%  in  the  a  and  y  phases  of  Ce, 
respectively.  The  RMTA  expression  is  reinterpreted  to  provide  an  understanding  of  the  low  values  of  ij  in  Ce  (0.7  to 
1.1  eV/A2)  in  spite  of  an  extremely  large  density  of  states  at  the  Fermi  level.  Using  moments  of  the  measured  phonon 
spectra,  the  superconducting  transition  temperature  is  calculated  and  compared  with  experiment.  The  RMTA  appears  to 
overestimate  17  by  approximately  23%  in  these  metals,  although  spin  fluctuations  may  be  contributing  to  this  apparent 
discrepancy  with  experiment. 


1.  Introduction 

In  his  classic  work  on  the  theory  of  the  elec¬ 
tron-phonon  (EP)  interaction  McMillan  [1] 
showed  that  the  EP  coupling  constant  A  can  be 
expressed  as 


A  » 


N(E,Xft 

\i(o> 


V 

M/co1)  ’ 


(1) 


where  S(Er)  is  the  electronic  (single  spin)  den¬ 
sity  of  states  at  the  Fermi  energy  Er,  </*>  is  a 
mean  square  electron-ion  matrix  element 
averaged  over  the  Fermi  surface,  Af  is  the  ionic 
man  and  (a*2)  is  an  appropriately  defined  [1] 
mean  square  phonon  frequency.  Specifically  (I1) 
is  given  by 


h  ' 


sv 

SR 


1 


x  5(£ta  -  Ef)8(Et  „  -  Ef)/N(EFf ,  (2) 


where  k,  n  are,  respectively,  the  wavevector  and 
band  index  of  the  state  |k,  n)  with  E»,„  and 
8V(r)/6R  is  the  change  in  crystal  potential  per 
unit  displacement  SR  of  an  ion  at  position  R.  Up 


to  this  point  the  expression  for  A  is  rigorous. 
However,  for  all  but  weak  scattering  ions  SV/SR 
is  very  difficult  to  compute  precisely,  and  further 
progress  in  the  understanding  of  (I2)  or  the 
McMillan-Hopfield  [2]  parameter  17  necessitates 
approximations  of  SV/SR. 

Noting  that  the  muffin-tin  approximation  for 
V(r)  has  been  very  successful  in  describing  the 
electronic  structure  of  transition  metals,  Gaspari 
and  Gyorffy  [3]  suggested  that  the  local  environ¬ 
ment  of  a  vibrating  strongly  scattering  ion  would 
be  modeled  realistically  by  rigidly  displacing  the 
muffin-tin  potential  for  that  ion.  Within  this  rigid 
muffin-tin  approximation  (RMTA),  Gaspari  and 
Gyorffy  obtained  the  relation  (in  Rydberg  units 
h  »  2m  »  e2/2  -  1  throughout  this  paper) 

x  sin2(5(  -  <5,4.i)v(V/*i  ,  (3) 

where  81  is  the  phase  shift  for  the  /th  partial 
wave  scattering  from  the  muffin-tin  potential. 
The  crystalline  enhancement  ratio  = 
N,(Ef)/N,(1)(Ef)  of  the  density  of  states  is  defined 
in  terms  of  the  single  scatterer  partial  density  of 
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states 

Ntx\Ef)  =  (2/  +  1)  f  drrRfc  EF) 

0 

=  ^(2/+l)rirr(rs,EF)|yt|.  (4) 

Here  r(  is  the  muffin-tin  radius,  R,  is  the  radial 
function  and  yt  is  the  energy  derivative  of  the 
logarithmic  derivative  of  the  radial  function 
evaluated  at  r,  and  Ef.  Eq.  (3)  has  been  applied 
widely  [3-5]  in  the  study  of  elemental  transition 
metals.  Often  it  is  found  that  the  results  of  the 
RMTA  in  this  form  seem  to  agree  with  empir¬ 
ically  determined  values  to  within  10-20%  and  to 
reproduce  trends  with  pressure  [6]  rather  well. 

Aside  from  the  approximation  of  rigidly  dis¬ 
placed  muffin-tin  potentials,  eq.  (3)  contains  two 
simplifying  assumptions.  The  first  is  the  so-called 
“spherical  band”  approximation  [3]  in  which  the 
angular  dependence  of  the  bands  is  ignored.  This 
is  not  an  essential  approximation,  and  recent 
studies  [5, 7]  have  generalized  eq.  (3)  ap¬ 
propriately  and  shown  the  corrections  to  be 
small  in  elemental  systems  with  cubic  symmetry. 
The  second  approximation  in  eq.  (3)  is  that  of 
ignoring  the  spin-orbit  interaction  (denoted  by 
the  superscript  ‘0’).  John  and  Harmann  [8]  have 
shown  how  this  restriction  can  be  removed 
(within  the  spherical  band  approximation)  but  to 
date  no  numerical  test  has  been  performed  to 
test  the  size  of  corrections. 

In  this  paper  the  RMTA  will  be  applied  to  the 
fee  phases  of  La,  Ce  and  Th.  each  of  which  has  f 
bands  [9-11]  near  or  at  EF.  Since  spin-orbit 
contributions  to  jj  may  not  be  negligible  in  sys¬ 
tems  with  large  atomic  number  (Z  *  57.  58  and 
90,  respectively,  for  these  atoms),  the  expression 
of  John  and  Hamann  will  be  evaluated  in  section 
2  and  compared  to  the  zero  spin-orbit  limit  in 
eq.  (3).  Contrary  to  the  expectations  of  Butler 
[5],  the  contribution  to  tj  from  d-to-f  scattering 
in  Ce  is  surprisingly  small  in  spite  of  a  huge  /  =  3 
partial  density  of  states  NKEf)-  This  finding  leads 


to  a  reinterpretation  [12]  of  the  contributions  to 
tj  in  terms  of  a  rearrangement  of  terms  in  eq.  (3). 
This  reinterpretation  and  a  discussion  of  its  im¬ 
plications  for  states  which  are  nearly  confined  to 
the  muffin-tin  sphere  are  given  in  section  3.  The 
results  for  the  superconducting  transition  tem¬ 
perature  To  which  is  found  to  be  overestimated, 
are  presented  in  section  4. 

2.  Spin-orbit  effects  in  the  RMTA 

The  expression  for  tj  given  by  John  and 
Hamann  [8],  which  includes  spin-orbit  effects,  is 

X  {^33  sin2(5f  -  SUi)vt vU\ 

+  2T+i  sin2(5<’  ~  vJ+i 

+  (21  -r  1)(2/  +  3)  “  $M)vf  ■  (5) 

The  superscripts  a  indicate  the  corresponding 
quantity  for  total  angular  momentum  j  »  /  ±  j 
and  the  “single  scatterer”  density  of  states  is 
given  by 

NP(EF)  =  (/  +  1)  j  dr  H|R  r(r.EF)P ,  (6a) 

0 

N?^(Ef)  =  /  J  dr  r2|/?r(rEP)|2 .  (6b) 

0 

The  large  component  Rf  of  the  radial  wave 
function  is  normalized  according  to 

R  f (r„  Ef)  =  ji(*Fr,)  cos  5?  -  n,(xFr,)  sin  8f  .  (7) 

where  *cF *  V£F.  in  terms  of  the  spherical  Bessel 
(ji)  and  Newmann  (n,)  functions.  Eqs.  (6)  can  be 
expressed  in  terms  of  logarithmic  derivatives  yf 
analogously  to  eq.  (4). 
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Eq.  (5)  has  been  evaluated  in  detail  for  La,  Ce 
and  Th  using  the  results  of  self-consistent,  spin- 
orbit,  linearized-augmented-plane-wave  band 
structure  calculations  [9-11].  To  compare  with 
the  corresponding  zero  spin-orbit  limit  the  fol¬ 
lowing  identifications  have  been  made: 

*i-[»r  +  (/+l)*r]/(2/+l).  (8a) 

X,{E,)-NT(Er)  +  NnEr).  (8b) 

,V|»(£p)  =  iVp(£F)  +  iVp(£p).  (8c) 

It  is  easily  shown  that  in  the  zero  spin-orbit  limit 
(Sf-*S,  and  vf-*  v,)  eq.  (5)  reduces  to  eq.  (3).  It 
is  instructive  also  to  introduce  the  quantities 

Ai  —  SJ  —  8J  ,  (9a) 

It,  *  [IvJ  +  (/  +  1K]/(2Z  +  1) ,  (9b) 

e,-vj-v*.  (9c) 

in  terms  of  which  the  relativistic  quantities  are 
given  by 

S'~S,~2iTTA''  (10a) 


s~i  =  s'+jnhA<'  (10b> 

and  similarly  for  vf  in  terms  of  v,  and  e,.  Res¬ 
pectively,  .If  and  €,  represent  the  spin-orbit 
splitting  of  the  phase  shift  and  the  crystalline 
enhancement  of  the  partial  density  of  states.  It 
has  not  been  emphasized  previously  that  the 
corrections  to  7?  in  eq.  (5)  which  are  first  order  in 
the  spin-orbit  corrections  1,  and  vanish  iden¬ 
tically.  This  may  be  due  to  the  traceless  nature  of 
the  spin-orbit  operator  and  it  is  in  keeping  with 
the  picture  that  quantitative  changes  in  observ¬ 
ables  tend  to  be  second  order  in  the  spin-orbit 
splitting. 

In  table  I  phase  shifts  5f,  crystalline 
enhancements  vf  and  partial  and  total  values  of 
Tf  (with  and  without  spin-orbit  effects)  are 
presented.  In  each  case  tj  is  increased  by  spin- 
orbit  corrections  (this  is  not  a  general  require¬ 
ment),  by  0.03,  0.08,  0.10  and  0.18eV/AJ  (or  1%. 
4%,  10%  and  33%,  respectively)  in  fee  La,  Th. 
a-Ce  and  y-Ce.  (This  is  a  fictitious  “paramag¬ 
netic  y-Ce”  which,  however,  illustrates  spin- 
orbit  effects  more  clearly  than  the  other  metals.) 
In  thorium  the  increase  arises  from  corrections 


Table  I. 

Relativistic  phase  shifts  5  f  (in  radians),  crystalline  density  of  states  enhancements  vf  and  the 
resulting  contributions  to  n(eV/Az);  the  s  sign  indicates  /  »  /  s ;  and  V0’  is  the  value  of  the 
McMillan-Hopfield  parameter  neglecting  spin-orbit  corrections 


La 

T-Ce 

or- 

-Ce 

Th 

«r.5f/-o 

_ 

-0.96 

-0.98 

_ 

-1.23 

- 

-1.14 

1 

-0.42 

-0.50 

-0.43 

-0.51 

-0.60 

-0.69 

-0.56 

-0.82 

2 

0.59 

0.54 

0.56 

0.51 

0.63 

0.58 

0.75 

0.61 

3 

0.026 

0.024 

0.142 

0.067 

0.187 

0.109 

0.095 

0.079 

i >T,v‘I  m  0 

- 

0.71 

- 

1.33 

- 

2.30 

- 

1.05 

1 

2.46 

2.65 

1.41 

1.50 

1.16 

1.17 

0.74 

1.01 

2 

1.20 

1.45 

0.63 

0.77 

0.48 

0.57 

0.58 

0.82 

3 

3.10 

3.61 

0.99 

7.01 

0.65 

2.75 

1.44 

2.24 

4 

7.86 

7.74 

6.89 

5.94 

4.30 

3.77 

3.12 

2.51 

lu* 

»P 

0.12 

0.12 

0.07 

0.07 

0.29 

0.29 

0.10 

0.10 

Pd 

1.58 

1.57 

0.24 

0.24 

0.41 

0.41 

0.93 

0.87 

df 

1.08 

1.08 

0.30 

0.14 

0.23 

0.15 

0.85 

0.83 

f* 

0.02 

0.02 

0.11 

0.09 

0.16 

0.15 

0.10 

0.10 

total 

2.81 

2.78 

0.72 

0.54 

1.10 

1.00 

1.99 

1.91 
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to  p-d  scattering  while  the  channels  involving  f 
states  are  unaffected.  In  the  two  phases  of 
cerium  the  increase  occurs  in  the  d-f  channel. 

The  spin-orbit  correction  is  seen  most  clearly 
in  y-Ce.  Using  an  obvious  notation  from  eq.  (5), 
=*  Vet  +  vi~  +  b£”.  can  first  be  noted  that, 
for  /3>2.  becomes  negligible  with  respect  to 
and  p'~  due  to  the  factor  [(2/  +  1)(2/  +  3)]“‘. 
The  large  increase  in  p  in  y-Ce  arises  from  the 
large  crystalline  enhancement  of  vf  at  the 
expense  of  vf.  This  leads  to  vi*  ~  0.28.  Va  = 
0.02  compared  to  their  zero  spin-orbit  counter¬ 
parts  0.08  and  0.05  (all  in  eV/A:).  Ultimately  the 
spin-orbit  increase  results  from  the  shifting  of 
states  at  Er  by  crystallinity  (solid  state  effects) 
from  one  member  of  the  pair  j  =  l±\  to  the 
other.  The  central  requirement  for  a  large  cor¬ 
rection  of  this  type  to  occur  is  that  v *  and  be 
as  dissimilar  as  possible  (i.e.,  that  one  be  maxi¬ 
mized  at  the  expense  of  the  other),  and  that  this 
skewing  be  in  the  same  direction  for  both  l  and 
/  +  1,  thereby  maximizing  either  or 

i/7v7.|.  To  obtain  the  maximum  ratio  i»7  requires 
that  the  crystal  density  of  states  (numerator)  and 
the  single  scatterer  density  of  states  (denomina¬ 
tor)  be  oppositely  skewed.  In  y-Ce  this  occurs 
only  for  /  =  3,  with  normal  ratios  (v‘i  =*  v~,)  for 
l  *  2  and  /  =  4.  This  opposite  skewing  of  the  f 
states  N7/JV7  =  13.3/7.0,  Np/Ni11*  =  1. 9/7.0 
leads  to  v\fv~t  =  7.08  and  a  factor  of  2.5  increase 
in  tj*. 

An  interesting  limiting  case  occurs  when  all 
the  crystal  states  (for  both  l  and  /  +  1)  are  shifted 
to  either  ;  =*  /  +  ;  or  /  =  /  -  i  Assuming  equal 
single  scatterer  densities  of  states  for  j  =  l±\,  the 
value  of  tju-i  is  increased  by  the  factor  (2/  + 
3)/(/+l).  Further  increase  may  result  from 
oppositely  skewing  the  single  scatterer  density  of 
states  as  noted  above,  which  will  require  a  nar¬ 
row  l  resonance  in  the  vicinity  of  £F. 

In  the  next  section  the  problem  of  under¬ 
standing  the  small  values  of  rj  in  Ce.  which 
possesses  an  extremely  large  /  density  of  states  at 
Ef.  is  addressed  in  terms  of  a  reinterpretation  of 
the  RMTA  expressions  eqs.  (3)  and  (5).  It  should 


be  noted  that  the  small  values  of  t)  calculated 
from  the  Ce  band  structure  persist  only  for 
values  of  Ef  lying  within  the  4f  bands.  If  Ef  is 
artificially  placed  above  the  f  bands  (say  0.05  to 
0.50  Ry  above),  larger  *■ alues  of  p  ~  3  to  5  eV/A; 
result,  just  as  obtained  from  a  similar  treatment 
using  the  La  band  structure. 

3.  Reinterpretation  of  the  RMTA  expression 

The  discussion  in  the  prc ,  ious  section  did  not 
address  the  cause  of  the  small  values  of  rj  in  Ce. 
compared  to  La,  in  spite  of  much  larger  values  of 
N(Ef).  The  following  rearrangement  and  rein¬ 
terpretation  of  the  RMTA  expression  for  17  or 
(I2)  provides  insight  into  this  puzzle. 

Begin  by  defining  a  rigid  muffin-tin  matrix 
element  for  the  /  -*  /  +  1  channel  by  (see  also 
Butler  [5]  for  the  non-relativistic  analog) 

ITT- 1  =  JRT  ^7  Kfl.rdr/rW:, ,  (11) 

0 

(r?)2  *  |  (RfYrdr,  <7  -  ±  .  (12) 

0 

Evidently  III-  \  is  independent  of  the  nor¬ 
malization  of  the  radial  functions  and  depends 
only  on  dV/dr  and  the  shapes  of  the  radial 
functions.  Using  the  normalization  in  eq.  (7)  the 
numerator  in  eq.  (11)  is  sin(5?  -  Sf-',).  and  the 
normalization  integral  can  be  expressed  as 

(Tlf=r;\RHr,.Eftf\yl\.  (13) 

where  y f  denotes  the  energy  derivative  of  the 
logarithmic  derivative  of  Rf.  Finally,  eq.  (11) 
becomes 

(Ilf.\Y  =  sin:(S7  -  Sr,)/[r7(RfR7:,):(yfy7:i)]- 

(14) 

In  terms  of  these  matrix  elements  <I:)  can  be 
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expressed  as 

</^*?<2/+lj(2/  +  3) 

x  {(2/  +  mtuifftfUi  +  (2/  +  3  Xnr-tffTfT+i 

+  T^T(/rr*1)2/r/r.,},  a  5) 

which  in  the  zero  spin-orbit  limit  becomes 

a*)'0'  =  £  (2/+  l)(2/  +  3)  •  (15’^ 

In  these  relations  (Z2)  is  expressed  in  terms  of  the 
fraction  /f  =  Nf(EF)/N(Er)  of  states  at  £f, 
and  //  =*  /f  +/r.  This  form  for  (Z2)  has  the  desir¬ 
able  property  [12]  that  the  atomic-like  quantity 
which  depends  only  on  the  togarithmic 
derivative  and  varies  smoothly  and  predictably 
with  Ef,  is  separated  from  the  crystalline  pro¬ 
perty  f.  Although  it  requires  a  full  band  struc¬ 
ture  calculation  to  ascertain  {f,},  these  fractions 
typically  vary  much  more  smoothly  [12]  than 
Ni(Ef)  and  N(Er)  separately  which  generally 
possess  Van  Hove  singularities  at  the  same 
energies.  This  behavior  is  strikingly  illustrated 
for  y-Ce  in  fig.  1  in  a  20  mRy  region  centered  on 
the  calculated  value  of  Ef.  In  this  region  there  is 
strong  hybridization  among  the  s,  p,  d  and  f 
bands  and  N(E)  varies  non-monotonically  by 
more  than  a  factor  of  2  whereas  the  fractions  /, 
vary  more  smoothly. 

The  form  for  Jy+i  given  by 

A/-i  *  sin(5,  -  S,^)/rm*i  (1 1') 

is  useful  for  understanding  the  contribution  of  f 
states  to  (I2)  as  EF  is  varied.  As  a  result  of  the 
boundary  condition  eq.  (7)  “fixing”  R,  at  r„  the 
normalization  integrals  rt  are  unbounded  (0< 

T i  <  °°)  and  they  become  the  principal  determin¬ 
ing  factors  in  the  magnitude  of  Zu„ ,.  Near  the 
bottom  of  the  bands  derived  from  the  atomic  / 


Fig.  1.  The  total  density  of  states  Oeft  scale)  and  partial 

wave  ratios  /  =*  0 . 4  (right  scale)  in  a  20  mRy  region 

centered  at  the  Fermi  energy  (0.450  Ry)  in  y-Ce.  In  this 
region  N(E)  has  fine  structure  and  varies  by  a  factor  of  two, 
whereas  the  ratios  ft  vary  much  more  smoothly. 


states,  Rt  peaks  at  or  near  r,  and  r,  will  be  small. 
Near  the  top  of  the  bands,  Rt  assumes  very  large 
values  for  r  <  r,  compared  to  R,(rt),  resulting  in 
large  values  of  r,  and  smaller  values  of  Ju*,  (and 
Ii-u).  Schematic  radial  functions  are  shown  in 
fig.  2.  The  boundary  condition  of  fixed  R,(rt) 
results  in  “the  dog  (r,)  being  wagged  by  the  tail 
[K/(r,)]”[13].  Physically,  the  small  values  of  /,;*• 
at  the  top  of  the  bands  reflects  the  poor  coupling 
to  plane  waves  of  states  which  are  too  strongly 
confined  to  the  muffin-tin  sphere.  (This  dis¬ 
cussion  is  somewhat  oversimplified:  the  “fixed” 
value  of  Ri(r,)  and  the  numerator  sin(5(  -  5,,i) 
both  vary  with  EF,  but  these  variations  usually 
are  secondary  effects.) 

This  description  of  r,  and  its  influence  on 
applies  to  wide  bands  such  as  s,  p  and  4d,  5d 
bands  in  transition  metals.  The  narrow  4f  bands 
in  La  and  Ce  (and  the  late  3d  bands)  present  a 
different  behavior  however.  For  these  metals  the 
f  bands  are  only  1-2  eV  wide  and,  even  at  the 
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Fig.  2.  Schematic  5d  radial  charge  densities  rR\(r,E)  illus¬ 
trating  the  bonding  behavior  at  the  bottom  of  the  band  and 
the  antibonding  behavior  at  the  top  of  the  band.  The  radial 
wavefunctions  (Ri)  have  been  normalized  to  the  same  value 
at  the  sphere  radius  (taken  as  3  a.u.). 

bottom  of  the  f  bands,  the  4f  radial  function  is 
strongly  confined  to  the  muffin-tin  sphere.  This  is 
illustrated  in  fig.  3  for  a-Ce,  where  £F  lies  near 
the  bottom  of  the  4f  bands.  Due  to  the  strong 
peaking  of  R,mJ  at  0.75  a.u.,  n.3  is  correspond¬ 
ingly  larger,  and  Iu*i  smaller,  than  for  the  s-p-d 
bands.  The  detrimental  effect  only  subsides  at  or 
below  about  2eV  below  the  4f  bands,  which  is 
where  £F  lies  in  La.  The  radial  charge  density  at 
£F  in  La  is  shown  in  comparison  with  that  of 
a-Ce  in  fig.  3,  normalized  such  that  the  radial 
functions  are  equal  at  r,.  Although  the  l  »  3 
radial  function  in  La  shows  a  mild  peak  at 
0.75  a.u.  due  to  the  incipient  4f  bands,  it  is 
insufficient  to  decrease  r3  significantly. 

The  squared  normalizations  at  £F  for  La,  Ce 
and  Th  are  included  in  table  II.  La  and  Ce  differ 
significantly  only  in  the  values  of  (rf)2,  which  for 
Ce  are  more  than  an  order  of  magnitude  larger 
than  those  of  La.  In  the  magnitude  of  r3,  as  in 
the  value  of  Te,  Th  is  much  more  similar  to  La 
than  to  the  isoelectronic  element  Ce.  This  of 


Fig.  3.  The  / »  3  radial  charge  densities  at  Ef  for  La  and 
a-Ce,  normalized  to  the  same  value  at  the  muffin-tin  sphere 
boundary.  In  La  it  is  clear  that  Er  is  far  off  the  4f  resonance, 
resulting  in  a  small  atomic  4f  character  (peak  at  0.75  a.u.). 
This  graph  clearly  illustrates  the  order  of  magnitude 
difference  in  normalization  integrals  t*  between  La  and  Ce. 

course  reflects  the  similar  position  of  the  f  bands 
with  respect  to  £F  in  La  [9]  and  Th  [11], 

The  small  df  contribution  to  17  in  Ce,  com¬ 
pared  to  both  La  and  Th,  can  be  attributed  to 
the  dramatically  decreased  value  of  lot  (see  table 
II).  This  small  value  is  illustrated  graphically  in 
fig-  4;  l\  is  two  orders  of  magnitude  smaller  than 
that  for  the  other  channels  in  a-Ce.  The 
extremely  small  value  results  from  £F  lying 
“near”  the  center  of  both  d  and  f  bands  giving 
large  values  of  r2  and  r3  as  well  as  a  reduced 
numerator  sin(5:  -  d3). 

A  second  important  feature  is  illustrated  in  fig. 
4.  The  quantities  Ih*i  are  seen  to  be  slowly  and 
smoothly  varying  with  energy,  while  r)u+,  reflects 
the  sharp  structure  which  may  occur  in  N(£F) 
and,  to  a  lesser  extent,  in  /».  The  fact  that 
can  be  written  solely  in  terms  of  the  phase  shifts 
(or  equivalently,  logarithmic  derivatives)  and 
their  energy  derivatives  [eqs.  (1 1)— (13)]  guaran¬ 
tees  their  smooth  variation  with  energy.  Fur- 
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1  able  II 

Relativistic  crystalline  density  of  stales  ratios  //.  normalization  integrals  ti  .  single  channel  matrix  elements  and  average  matrix  elements  (in 
eVJ/A2)  lor  La,  t'e  and  lb;  the  t  sign  indicates  /  =  ft }  and  (I2)0  is  the  value  of  (I2)  neglecting  spin-orhit  corrections 
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Fig  4.  (a)  The  squared  matrix  element  lij-i  lor  scattering 
from  /  to  / *  1  versus  energy  Note  that  l'*  is  smaller  by  a 
factor  of  100  than  for  the  other  channels,  and  also  that  the 
variation  with  energy  is  smooth  for  all  channels,  (b)  Con¬ 
tributions  iju- i  to  tj  (left  scale)  and  average  squared  elec¬ 
tron-ion  matrix  element  (I1)  (right  scale)  versus  energy.  Both 
</:i  and  (to  a  greater  extent)  n  show  rapid  variation  wuh 
energy  between  2  and  4  mRy  which  comes  solely  from  den¬ 
sity  of  states  factors. 


therniore,  its  atomic-like  nature  suggests  that  for 
a  given  atom  It* t  will  possess  a  transferability, 
frctii  element  to  compound  and  compound  to 
compound,  which  the  crystal  structure  related 
quantities  ft  and  v,  can  never  attain.  These  pro¬ 
perties  have  been  investigated  in  detail  else¬ 
where  [12]  for  the  case  of  niobium  and 
vanadium. 


4.  The  transition  temperature 

In  this  section  the  calculated  values  of  Te,  with 
and  without  the  spin-orbit  correction  discussed 
in  Section  2.  are  compared  with  experimental 
values.  For  Te  we  use  the  Alien-Dynes  [14] 


modification  of  the  McMillan  [l]  equation 

t  t  ,i  -,s  J  1 .04(1  +  A)  1 
7 «  =  (‘•W1  -)  exP\  "  A  -  M  *(1  +  0.62A )  J 

with  wio,  defined  as  usual  [14]  and  fi*  -  0.13.  For 
the  phonon  moments  coy*  and  (at1)  of  Th.  the 
Butler  [5]  prescription 

oiioi  *  0-60  0e , 

(u>yn  =  o.69  ec . 

has  been  applied,  with  the  high-temperature 
Debye  temperature  0C  taken  from  the  inelastic 
neutron  scattering  measurements  of  Reese  et  al. 
[IS].  The  same  prescription  has  been  used  [16] 
for  La  and  Ce. 

Table  III  summarizes  the  results.  (Note  that  no 
comparison  is  possible  between  the  theoretical 
“paramagnetic  y-Ce"  and  the  actual  high  tem¬ 
perature  magnetic  phase.)  In  each  of  the  metals 
Tc  is  substantially  overestimated.  It  is  not  the 
purpose  of  the  present  paper  to  investigate  these 
discrepancies  in  detail,  but  a  few  observations 
can  be  noted.  For  La  and  Th,  the  empirical 
values  of  A  which  give  the  experimental  value  of 
7*c  (1.15  and  0.60.  respectively)  suggest  that  the 
RMTA  overestimates  tj  by  0.5-0.6eV/A:  in 
these  similar  metals.  These  overestimates  have 

Table  III 

Experimental  and  theoretical  quantities  determining  the 
superconducting  transition  temperature  Tc.  The  superscript 
zero  indicates  quantities  calculated  disregarding  spin-orbit 
corrections  to  the  RMTA  expressions 


La 

y-Ce 

<*-Ce 

Th 

uho«(K) 

"4 

-o 

104 

86 

(w^'OC) 

86 

81 

119 

99 

■Vf'ur'.ieV  A:t 

t  .84 

0.92 

1.99 

2.28 

i.Av;  A:) 

2.65.2.62 

0.72.0.54 

1.10.1.00 

1.99.1.91 

1.44.1.42 

0  7-.0  58 

0.55.0.50 

0.87.0.84 

T.-.Tf'(K) 

84.8.3 

2.8. 1.0 

1.3.0. 8 

3.9,36 

up 

6.05 

- 

-001 

14 
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been  discussed  elsewhere,  for  La  by  Pickett  et  al. 
[9]  and  for  Th  by  Winter  [17].  For  a-Ce  a  good 
case  can  be  made  [10]  for  invoking  a  paramag- 
non  contribution  A,pin  detrimental  to  spin  pairing, 
which  need  only  be  of  the  order  of  0.1S  to 
explain  the  overestimate  of  Tc  by  two  orders  of 
magnitude.  Indeed,  recent  careful  estimates  by 
Rietschel  and  Winter  [18]  point  to  a  probable 
value  of  »  0.2  for  Nb  and  V.  so  this  effect 
may  be  contributing  to  the  overestimate  of  Tc  in 
La  and  Th  as  well. 


5.  Conclusions 

The  specific  conclusions  of  this  study  are  two¬ 
fold.  Firstly,  spin-orbit  corrections  to  rj  (within 
RMTA)  may  become  important  in  metals  with 
large  atomic  numbers  Z.  One  condition  which 
can  lead  to  an  appreciable  correction  is  for  the 
spin-orbit  interaction  to  be  large  enough  to 
result  in  a  net  displacement  of  ;'  =  /  +  ]  states 
with  respect  to  j  ~  l-  \  states,  viz.,  a  deviation  of 
the  ratio  fit  ft  from  the  zero  spin-orbit  value 
1/(1  +  1).  These  deviations  are  small  (see  table  II) 
in  La  (Z  *  57)  and  even  Th  (Z  =  90)  at  Ef.  Due 
to  the  proximity  of  Er  to  the  4f  bands,  deviations 
in  Ce  become  significant.  Perhaps  a  more  import¬ 
ant  criterion  for  spin-orbit  correction  to  r;  may 
be  that  the  spin-orbit  splitting  J£,  of  the  /- 
resonance  [<l£i  =  |E(*  -  Erl,  where  SJ(Et)  = 
S/(E r)],  be  comparable  to  the  distance  of  the  / 
band  center  from  Ef.  This  leads  to  a  variation  of 
the  ratio  r/lrt  from  unity  (table  II)  which  is 
large  for  Ce.  As  a  result  the  matrix  elements 
It/- 1  and  lull  differ.  In  the  systems  studied  here 
It/-i IJ/*u  and  ft  is  generally  increased  at  the 
expense  of  ft.  resulting  in  increased  values  of  tj. 

Secondly  for  each  of  the  “f-band  metals”  La. 
Ce  and  Th,  the  RMTA  appears  to  overestimate 
Tf  by  approximately  25%.  In  this  respect  these 
f-band  metals  resemble  d-band  metals,  where 
similar  overestimates  [4,  5, 19]  by  the  RMTA  are 
common. 
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Transferability  and  the  electron-phonon  interaction: 

A  reinterpretation  of  the  rigid-muffin-tin 
approximation 
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The  expression  for  the  McMillan-Hopfield  constant  tj  =  .V(EfK/J)  within  the  rigid- 
muffin-tin  approximation  (RMTA)  is  reinterpreted  in  terms  of  single-channel  electron-ion 
matrix  elements  Iu+\  and  angular-momentum-character  fractions//  of  the  electronic 
states  at  the  Fermi  level.  Reexamination  of  Nb-  and  V-based  ^15  compounds  suggests 
that  lu+\  is  more  nearly  an  atomic  property,  and  thus  transferable  from  system  to  sys¬ 
tem,  than  other  commonly  used  quantities.  The  fractions  //  are  dependent  on  bonding 
character  and  crystal  structure  but  tend  to  be  constant  within  a  class  of  compounds.  Cri¬ 
teria  for  increasing  rf  within  RMTA  are  discussed. 


I.  INTRODUCTION 

In  his  classic  study  of  the  electron-phonon  in¬ 
teraction  parameter  k  and  its  relation  to  the  super¬ 
conducting  transition  temperature  Te,  McMillan1 
showed  that  A.  can  be  written 


WiEF)U2) 

M(tit2) 


(1) 


in  terms  of  the  density  of  states  N(EF)  per  spin  at 
the  Fermi  energy  EF,  a  mean-square  electron-ion 
matrix  element  (.I2),  and  an  appropriately  defined 
lattice  stiffness  M(co2).  A  topic  of  particular  in¬ 
terest  in  the  ensuing  decade  has  been  the  discussion 
of  high  Tc  metals  in  terms  of  one  or  the  other  of 
the  factors  N(EF),  {I2),  or  Mia2)  as  most  im¬ 
portant  in  leading  to  high-temperature  supercon¬ 
ductivity.  With  one  of  these  factors  [say,  .V(£f )] 
as  the  principal  determinant  of  the  high  Te,  a  na¬ 
tural  procedure  which  could  be  tried  to  produce  a 
higher  temperature  superconductor  is  to  "transfer” 
this  large  value  of  N(EF)  to  a  metal  with  more 
favorable  values  of  </2)  and  M(cj2). 

The  limited  data  available  to  McMillan1  [all  em¬ 
pirical  except  for  a  few  calculated  values  of 
S(Ef )]  suggested  that  the  product  S(EF)(I2)  is 
roughly  constant  within  a  class  of  materials,  in 
which  case  the  search  for  high  Te  materials  should 
concentrate  on  soft  lattices  with  small  values  of 
M(ot2).  Although  this  approach  was  used  by 
several  investigators  in  the  following  years,  often 
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with  apparent  success  in  understanding  trends  in 
Te,  more  recently  it  has  been  called  into  question 
by  a  number  of  calculations.  Based  upon  the 
"rigid-muffin-tin  approximation”  (RMTA)  of 
Gaspari  and  Gyorffy,2  Papaconstantopoulos 
et  al.,3  and  Butler4  have  shown  that,  within  iso- 
structural  elemental  transition  metals,  N(EF)(I2) 
can  vary  by  as  much  as  a  factor  of  3.  Further¬ 
more,  these  estimates  appear  to  agree  rather  well 
with  more  recent  experimental  data.  Thus  varia¬ 
tions  of  N(Ef)  and  (I2)  separately  must  be  taken 
into  account. 

Another  early  attempt  to  correlate  superconduct¬ 
ing  properties  was  made  by  Hopfield,5  who  intro¬ 
duced  the  notational  convenience  rf  =  .V(£F)</:) 
(the  "McMillan-Hopfield  parameter”).  Analysis  by 
Hopfield  which  emphasized  only  the  p-d  scattering 
led  him  to  anticipate  that,  in  transition  metals 
which  had  a  large  d  partial  density  of  states,  tj 
would  be  essentially  an  atomic  parameter.  As 
such,  rf  would  be  transferable,  from  element  to  al¬ 
loy  and  within  a  class  of  compounds,  and  Hopfield 
used  this  idea  in  an  attempt  to  understand  super¬ 
conductivity  in  transition-metal  alloys  and  within 
A 1 5  compounds.  However,  the  startling  discovery 
by  Klein  and  Papaconstantopoulos6  that  the  d-f 
contribution7  to  17  is  not  only  appreciable  but  in 
fact  dominant  in  transition  metals,  and  the  ensuing 
realization  that  the  /  contribution  is  a  property  of 
the  environment  rather  than  atomic  in  nature,  has 
pointed  out  the  limited  usefulness  of  Hopfield's  ap¬ 
proach. 

There  have  also  been  suggestions  that,  within 
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limited  classes  or  subclasses  of  compounds,  there 
may  be  an  “atomic  k"  which  is  more  or  less 
transferable.  Phillips'  suggested  that,  within  the 
NaQ-structure  NbN  family  of  compounds,  the 
sum  of  atomic  contributions  k  =  kA+kB  should  be 
a  useful  concept.  However,  RMTA  calculations  by 
Klein,  Papaconstantopoulos,  and  Boyer9  indicate 
kc  for  carbon  varies  by  more  than  a  factor  of  5  in 
the  compounds  NbC,  TaC,  and  HfC.  This  idea  of 
“transferable  k"  was  then  further  restricted  to  the 
subclass  “NaCl-structure  compounds  with  phonon 
anomalies”  by  Haufe,  Kerker,  and  Benneman.10 
However,  at  such  a  specialized  level  the  concept  of 
transferability  loses  much  of  its  usefulness  in 
understanding  superconductivity  as  well  as  for 
predicting  higher  Tc  materials. 

In  view  of  these  largely  unfruitful  attempts  to 
identify  transferable  quantities  related  to  supercon¬ 
ductivity,  and  also  the  difficulty  in  calculating  the 
phonon  spectrum  and  hence  the  lattice  stiffness 
M(a)1),  emphasis  has  shifted  to  attempting  to 
understand  the  behavior  of  rj  within  rigid-ion 
models,  chiefly  the  RMTA.  Notably,  the  group  at 
the  Naval  Research  Laboratory  has  published  ex¬ 
tensive  RMTA  calculations  for  elements3-6  and  for 
NaCl-structure9"  compounds  and  A\5  (Ref.  12) 
compounds,  among  others.  ButleT4  has  provided  a 
detailed  RMTA  study  of  the  Ad  transition  metals, 
noting  in  particular  the  dependence  of  </2>  on 
atomic  number  (i.e.,  valence)  and  volume  in  this 
class  of  metals.  A  simultaneous  examination  of 
the  behavior  of  17  across  the  Ad  series  was  given  by 
Pettifor,13  who  related  RMTA  quantities  to  band- 
structure  parameters. 

In  this  paper  I  present  a  reinterpretation  of  the 
original  Gaspari-Gyorffy  expression  for  </2) 
which  results  from  a  simple  regrouping  of  terms. 

In  addition  to  providing  more  physical  insight  into 
the  quantities  which  determine  (I1),  this  rein¬ 
terpretation  has  two  other  favorable  consequences. 
One  is  that  atomiclike  quantities  are  identified, 
these  being  the  electron-ion  matrix  elements  lu  +  \ 
for  scattering,  from  partial  wave  /  to  partial  wave 
/  -V 1,  by  an  atomic  displacement.  The  atomiclike 
nature  of  7/,/  +  i  suggests  an  approximate  transfera¬ 
bility  which  seems  to  hold  at  least  in  systems  with 
common  bonding  characteristics.  The  other  is  that 
the  fraction  of  states  //  =.V,(£f  )/N(EF)  with  an¬ 
gular  momentum  /,  which  multiplies  I/j  +  t  is  a 
much  smoother  quantity  (in  several  senses  to  be 
discussed  below)  than  the  original  crystalline 
enhancement  factors.  Results  of  various  previous 
calculations  are  reviewed  in  terms  of  these  ideas. 


II.  REARRANGEMENT  OF  THE 
GASPARI-GYORFFY  EXPRESSION 


Gaspari  and  Gyorffy  derived  an  expression  for 
(I2)  which  can  be  written  (in  atomic  units 
fl-2 m  =el/2  —  1) 


<72> 


Ef 

tr.V[EF)2 


2 


2I/  +  1) 


Xsin2(5/— S,  +  1)v,v,  +  1  ,  (2) 

where  5;  is  the  Fermi-energy  phase  shift  for  partial 
wave  /.  The  “crystalline  enhancement”  of  the  den¬ 
sity  of  states  is  given  by 

V;  =  :V/  ( Efr  )  /Nj  U(Ef)  , 

where  .Vj  "(£f )  is  the  corresponding  density  of 
states  for  a  single  scatterer  rather  than  for  a  lattice 
of  scatterers.  Specifically, 

1 1  \  V  Ef 

iV;  (£f )= - H2/  +  1) 

TT 

X  fQSdrr2R2(r,EF)  ,  (3) 

where  Rt  is  the  radial  wave  function  and  Rs  is  the 
muffin-tin  radius.  The  trigonometric  factor  in  Eq. 
(2)  arises  from  the  RMTA  matrix  element 

Jo  drr2R,— £,  +  |=sin(5,-5,  +  1)  ,  (4) 

with  the  radial  function  normalization  chosen  to  be 

RiiRs,EF)  =  j,(KFRs  )cos5, 

— n/(Kf R$)sin6/  .  (5) 

Here  j;  and  nt  are  the  spherical  Bessel  and  Neu¬ 
mann  function  and  kf  =  Ef. 

Several  features  of  Eq.  (2)  should  be  noted. 

First,  the  factor  sin:(5;  — 5;t])  which  is  bounded 
between  0  and  1,  depends,  from  Eq.  (4),  on  the 
muffm-tin  potential  V,  on  its  overlap  with  radial 
functions,  and  on  the  radial  function  normaliza¬ 
tion,  which  is  not  bounded  [see  Eq.  (6)  below]. 
Secondly,  the  rather  unphysical  single  scatterer 
density  of  states  tends  to  obscure  the  physical  in¬ 
terpretation.  Thirdly,  the  factor  EF ,  which  seems 
to  imply  an  oven  dependence  on  the  zero  of  ener¬ 
gy,  is  in  fact  artificial,  being  canceled  by  the  V~Ef 
factors  from  Eq.  (3).  This  last  feature  points  to 
the  fact  that  the  factor  (t/£f/“)  (21  —  1)  was  in¬ 
troduced  into  the  expression  (2)  to  convert  the  nor¬ 
malization  integral 
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’  r*s , 

rj=  J0  drr-Rf  (6) 

to  a  density  of  states. 

Although  the  original  expression  (2)  may  be 
desirable  from  some  points  of  view,  the  following 
reinterpretation  in  many  ways  provides  a  simpler, 
and  thereby  more  illuminating  picture.  Begin  by 
defining  a  rigid-muffin-tin  matrix  element  for  the 
/ — /  -j- 1  channel  by  (see  also  Butler4! 

r&s  dV 

hi* i  =  J0  drr'-Ri  —  R^yArtr,^^.  (7) 


Evidently  Iu + j  is  independent  of  the  normaliza¬ 
tion  of  the  radial  functions  and  depends  only  on 
dV  /dr  and  the  shapes  of  the  radial  functions.  This 
form  of  matrix  element  has  been  used  previously 
by  Pickett  and  Gyorffy,14  Pettifor,13  and  to  some 
extent  by  Butler,4  and  was  also  preferred  by  Al¬ 
len.11  In  terms  of  these  matrix  elements  and  the 
fraction  //  of  states  of  angular  momentum  /,  (I1) 
can  be  expressed  as 


</J>*2 


(21+2) 


i  (21  + 1)(2/  4-3) 


hi*Jifi*i 


(8) 


This  expression  provides  a  more  straightforward 
and  useful  interpretation  for  primarily  two  reasons. 
The  first  reason  is  that  +  |  will  be,  to  a  degree  to 
be  discussed  in  Sec.  Ill,  an  atomic  property  as  it 
depends  only  on  the  muffin-tin  potential,  through 
its  derivative  dV  /dr  and  its  radial  functions,  but 
not  explicitly  on  the  crystalline  arrangement  of 
neighboring  atoms.  This  in  fact  is  essentially  the 
idea  of  Hopfield,1  but  restricted  to  the  muffin-tin 
region  and  making  no  attempt  to  simplify  the  ra¬ 
tios  /;  out  of  the  expression. 

The  second  desirable  feature  of  the  summand  in 
Eq.  (8)  is  its  dependence  on  energy  EF.  From  Eqs. 
(4)  and  (7)  and  the  expression 


where  Yi  is  the  energy  derivative  of  the  logarithmic 
derivative,  J3/+I  is  seen  to  be  very  smooth,  varying 
as  do  6/  and  yt  on  the  scale  of  the  /  and  (/  + 1 ) 
bandwidth*.  This  energy  variation  has  been  dis¬ 
cussed  previously  by  Pettifor.13  Although  it  is 
somewhat  less  obvious,  the  ratios  ft  are  much 
more  smoothly  varying  then  either  its  numerator 
or  denominator  alone.  This  is  the  result  of  (i) 

<V/(£)  and  N(E)  having  canceling  van  Hove  singu¬ 
larities  at  the  same  energies,  and  (ii)  ft  depending 
only  on  the  character  of  eigenstates,  rather  than  on 
the  density  of  states  which  may  vary  rapidly.  It 
follows  that  (I1)  is  slowly  varying  with  energy,  a 


characteristic  which  was  not  evident  in  the  original 
Gaspari-Gyorffy  expression. 

The  behavior  of  ft  in  La  (Ref.  16)  is  shown  in 
Fig.  1  and  it  can  be  compared  with  .V/  and  .V  in 
Fig.  2.  La  is  a  good  metal  in  which  to  illustrate 
this  behavior  because,  although  over  most  of  its 
spectrum  it  is  typical  of  a  5d  transition  metal,  it 
possesses  4/  (/=  3)  bands  centered  2.5  eV  above  £> 
which  provide  an  extreme  example  of  the  features 
mentioned  above.  Figures  1  and  2  show  dramati¬ 
cally  how  the  sharply  structured  behavior  of  .Vji  £' 
is  translated  into  a  smooth,  almost  resonance-like 
behavior  of  /j(£).  Likewise,  the  structured  d  spec¬ 
trum  .V:(£)  is  converted  into  a  smoother  form  in 
/>(£)  which,  however,  is  nearly  split  by  the  / 
bands.  [In  the  absence  of  the  /  bands,  as  in  more 
typical  transition  metals,  only  a  minor  minimum  in 
/»(£)  will  occur  in  the  low  density-of-states 
bonding-antibondina  “gap.”] 

Deviation  of  /;(£)  from  very  smooth  behavior 
results  primarily  from  hybridization,  as  mentioned 
above  for  the  d-f  interaction  effect  on  /,(£).  This 
effect  is  more  dearly  drawn  in  f0(E )  and  / ( ( £), 
each  of  which  show  structure  due  to  hybridization 
with  d  bands  (at  —2  eV)  and  /  bands  ( — 1.5  eV). 
However,  the  resulting  variation  with  energy  in 
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FIG.  1.  The  electronic  density-of-states  ratios  f/ 
versus  energy  for  fee  La,  /«0,  I,  2,  3,  4.  Successive 
plots  are  displaced  by  one  unit  (0$//  5 1).  The  zero  of 
energy  is  fixed  at  the  Fermi  level. 
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Lanthanum 


ENERGY  (eV) 


FIG.  2.  The  total  and  partial  densities  of  states  of  fee 
La.  Rescalings  of  the  various  partial  densities  of  states 
are  as  noted. 

each  case  remains  much  smaller  than  that  of 
X,(E). 

To  give  a  full  picture  of  the  present  approach  to 
</:),  the  quantities  8/t  rh  and  I2j+l  for  La  are 
shown  in  Fig.  3.  [For  8/  and  rt  the  <2y  +  D- 
weighted  averages  of  the  corresponding  fully  rela¬ 
tivistic  quantities  for  j  =l±  y  have  been  used  for 
/=  1,  2,  and  4.  For  /= 3  this  averaging  results  in 
unphysical  behavior  in  the  energy  region  between 
the  j  =  4  and  y  =  y  resonances,  so  only  the  j  =  y 
quantities  are  shown  in  Fig.  3.  The  full  relativistic 
treatment  of  </2>  will  be  presented  elsewhere.17] 

The  phase  shifts  for  /= 0,  1,  and  2  are  typical  of 
an  early  5 d  transition  metal,  but  the  sharp  /=  3 
resonance  (where  53=it/2)  at  2.1  eV  gives  rise  to 
flat  4/ bands  not  present  (or  rather,  fully  occupied) 
in  heavier  5 d  transition  metals.  The  /=4  phase 
shift  is  positive  but  less  than  7x  10” J  throughout 
this  range. 

The  normalizations  rt  calculated  from  Eq.  (6) 
and  shown  in  Fig.  3(b)  are  proportional  fexcept  for 
the  VEf  factor)  to  the  “single  scatterer”  density 
of  states  used  in  the  original  Gaspari-Gyorffy  ex¬ 
pression.  For  /= 0,  1,  and  4  these  are  monotonic 
in  the  interesting  energy  range.  For  1=2,  r2  peaks 
somewhat  below  the  d-band  center  due  to  the  clear 
nonresonant  form  of  82  in  Fig.  3(a).  The  sharp 
1= 3  resonance  results  in  a  huge  increase  in  r3  in 
the  4/-band  region.  The  effect  of  this  resonance 
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FIG.  3.  La  functions  entering  the  RMTA  expression 
for  17.  (a)  Phase  shifts  6/  relative  to  (b)  Semilog  plot 
of  the  normalization  integrals  r,.  (c)  Semilog  plot  of  the 
single-channel  matrix  elements  lit  *  i-  Note  the  zero  of 
l\,}  at  the  point  where  6;=6). 


on  Iu  +  \  is  evident  in  Fig.  3(c),  where  l\j  drops  by 
over  an  order  of  magnitude  in  this  region;  in  addi¬ 
tion,  / 2t3  vanishes  where  8:=8}.  However,  in  the 
absence  of  the  4/ bands,  Iu  +  \  are  smooth  over  the 
entire  d-band  region. 


III.  TRANSFERABILITY  OF  /,?/  +  l 

In  this  section  we  investigate  the  degree  to 
which  the  matrix  elements  Ifj  + 1  represent  atomic- 
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like,  transferable  quantities.  For  this  purpose  we 
focus  on  Nb  and  V,  for  which  Papaconstanto- 
poulos,  Klein,  and  Boyer  (Refs.  3,  6,  9,  11,  and  12) 
have  carried  out  self-consistent  band-structure 
determinations  and  calculations  of  r/  for  the  ele¬ 
mental  metals,  several  /115-structure  compounds, 
and  two  Nb- based  NaCl-structure  metals.  These 
compounds  are  listed  in  Table  I. 

A  feature  which  was  not  mentioned  in  the  previ¬ 
ous  section  is  the  dependence  of  the  various  quan¬ 
tities  on  the  sphere  radius  Rs ■  As  far  as  the 
RMTA  itself  is  concerned,  Rs  should  be  fixed  at 
the  radius  which  most  nearly  reproduces  the  actual 
scattering  properties.  However,  theoretical  ambi¬ 
guities  can  arise  if  spheres  from  distinct  atoms 
overlap,  and  in  the  NRL  calculations  which  will  be 
used  here,  maximal  nonoverlapping  spheres  have 
been  used.  Each  of  the  quantities  8/,  .V,,  and  .V,"1 
depend  on  Rs,  but  the  ratio  v;  =  iV,/.V,ul  is  in¬ 
dependent  of  Rs  to  lowest  order  and  the  Rs  depen¬ 


dence  of  8/  is  often  weak.  The  resulting  weak 
dependence  of  rj  on  Rs  is  central  to  the  usefulness 
of  the  RMTA  in  that  it  provides  a  relatively 
unambiguous  result  to  compare  with  experiment. 

The  quantities  /;  and  (,  the  intermediate 
quantities  in  the  interpretation  of  RMTA  being 
proposed  here,  are  separately  dependent  on  Rs. 

For  a  comparison  of  these  quantities  in  various 
systems,  renormalization  to  a  common  radius  Rs 
must  be  carried  out  first.  This  has  been  done  as 
follows.  For  1=0,  1,  and  2  the  radial  wave  func¬ 
tions  are  varying  slowly  near  Rs  and  the  normali¬ 
zation  correction  Arj  is  approximated  by 

g 

Arf=  f  Sr2drR,2(r) 
j  «s 

=Rf{Rs)(Rs-R^)/ 3  .  (101 

In  addition  the  radial  function  Rt  itself  is  renor¬ 
malized  to  the  appropriate  value  at  the  new  radius 
Rs.  Explicitly, 


r?(Rs)aE[rftRs)+Arf] 


j,(KFRs  )cos6/ — nt(KF  Rs  )sin  8/ 
j/(KFRs  )cos6] — nt(KFRs  )sin8/ 


(11) 


For  /  3*3  the  fact  that  R}(r)*ij}(KFr)  (inside  as 
well_as  outside  the  sphere)  was  used  to  recalculate 
t]{Rs)  directly.  No  correction  to  the  phase  shifts 
has  been  considered. 

The  radius  Rs  was  chosen  to  be  midway  between 


I 

the  elemental  value  of  Rs  and  the  mean  of  the 
/415-structure  values  given  in  Table  I;  these  were 
taken  as  Rs=  2.57  and  2.37  a.u.  for  Nb  and  V, 
respectively.  The  resulting  values  of  rji  Rs '  in  the 
.4 1 5  compounds  were  5  —  20  %  increases  over 


TABLE  I.  Sphere  radius  Rs  (a.u.),  Fermi  energy  EF  (Ry),  density  of  states  per  metal  atom  per  spin  .V  *(£f)  iRy-1), 
phase  shifts  8 f  (radians)  and  density-of-states  ratios  //,  and  enhancements  v,,  for  Nb-  and  V-based  compounds  discussed 
in  the  test. 


Rs 

Ef 

nhef) 

So 

5, 

Sj 

Ss 

fo 

/. 

/: 

ft 

v0 

V| 

Vv 

vj 

NbC 

2.313 

0.869 

4.76 

-1.11 

-0.51 

1.30 

0.0107 

0.0086 

0.0123 

0.383 

0.0288 

0.193 

0.150 

0.361 

2.571 

NbN 

2.308 

0.849 

6.43 

-1.01 

-0.47 

1.72 

0.0111 

0.0057 

0.0034 

0.412 

0.0234 

0.161 

0.056 

0.613 

2.945 

NbjAl 

2.450 

0.837 

16.60 

-1.13 

-0.50 

0.86 

0.0076 

0.0115 

0.0613 

0.530 

0.0123 

0.965 

3.257 

1.707 

4.248 

NbjGa 

2.443 

0.851 

15.98 

-1.14 

-0.50 

0.90 

0.0082 

0.0110 

0.0668 

0.529 

0.0121 

0.887 

3.341 

1.608 

3.827 

NbjSi 

2.457 

0.832 

6.97 

-1.12 

-0.49 

0.90 

0.0079 

0.0069 

0.0557 

0.560 

0.0120 

0.242 

1.241 

0.721 

1.747 

NbjGe 

2.438 

0.866 

8.89 

-1.15 

-0.51 

0.92 

0.0086 

0.0061 

0.0583 

0.563 

0.0126 

0.271 

1.597 

0.967 

2.108 

NbjSn 

2.495 

0.833 

13.22 

-1.11 

-0.49 

0.94 

0.0083 

0.0047 

0.0544 

0.594 

0.0102 

0.306 

2.367 

1.367 

2.73 1 

NbjSb 

2.482 

0.874 

4.33 

-1.15 

-0.51 

1.02 

0.0094 

0.0046 

0.0502 

0.594 

0.0137 

0.099 

0.687 

0.446 

1.052 

Nb 

2.685 

0.676 

9.71 

-0.93 

-0.36 

1.14 

0.0074 

0.0178 

0.0746 

0.724 

0.0128 

0.807 

2.923 

0.681 

3.846 

V 

2.477 

0.675 

12.70 

-0.69 

-017 

1.03 

0.0027 

0.0102 

0.0638 

0.846 

0.0071 

0.549 

2.841 

0.634 

6.936 

VjAl 

2.282 

0.819 

15.74 

-0.88 

-0.28 

0.95 

0.0026 

0.0047 

0.0543 

0.638 

0.0056 

0.365 

2.698 

1.003 

4.288 

VjGa 

2.276 

0.827 

24.63 

-0.88 

-0.28 

1.02 

0.0028 

0.0037 

0.0565 

0.652 

0.0044 

0.445 

4.333 

1.539 

5.173 

VjSi 

2.231 

0.889 

16.68 

-0.94 

-0.31 

0.98 

0.0031 

0.0030 

0.0568 

0.629 

0.0063 

0.254 

2.795 

1.224 

4.130 

VjGe 

2.253 

0.871 

9.56 

-0.92 

-0.J0 

1.02 

0.0031 

0.0038 

0.0415 

0.656 

0.0071 

0.183 

1.183 

0.670 

2.776 

V,Sn 

2.355 

0.803 

10.25 

-0.84 

-0.26 

1.07 

0.0029 

0.0028 

0.0317 

0.739 

0.0052 

0.133 

1.043 

0.617 

2.665 
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-rjiRs)  for  /=* 0,  1,  and  3  and  5—10%  decreases 
for  /=2.  The  corresponding  changes  in  elemental 
Mb  and  V  were  of  opposite  sign  in  each  case,  since 
8 s  <  Rs  in  these  cases.  It  is  notable  that  these  re¬ 
normalizations  resulted  in  bringing  ij  for  the  ^4 15 
compounds  closer  to  the  corresponding  elemental 
value  only  for  1=0  and  1.  For  1=2  the  renormali¬ 
zation  worsened  this  “correspondence,”  while  no 
clear  trend  appeared  for  1=2.  In  any  case  these  re¬ 
sults  seem  to  reflect  real  similarities  or  differences 
between  the  Nb-  and  V-atom  environments  in  the 
elements  and  the  compounds. 

In  Tables  I  and  II  an  extensive  listing  of  the 
quantities  which  enter  into  the  determination  of 
Nb-  or  V-  atom  contributions  to  tj  is  presented.18 
Both  iV/  and  Nt  have  been  renormalized  to  Rs  by 
multiplying  by  the  ratio  r](8s)/-rj{Rs),  hence  17  is 
unaffected.  Although  corrections  to  the  phase 
shifts  could  have  been  incorporated,  this  would 
amount  to  investigating  the  (small)  Rs  dependence 
of  77,  which  is  not  the  purpose  of  this  paper.  It 
should  be  noted  that  the  extrapolation  from  Rs  to 
Rs  is  somewhat  more  uncertain  in  NbC  and  NbN 
than  for  the  other  cases,  due  to  the  interval 
Rs  —  Rs  being  twice  as  large. 

The  principal  result  here  is  the  regularity  in  the 
values  of  Iu+ltNb)  and  //?/+ j(V)  from  system  to 
system,  greater  than  that  of  the  numerator 
sin2(6,—  6/  +  |)  or  denominator  rfr}+l  separately, 
each  of  which  is  dependent  on  radial  function  nor¬ 
malization.  The  deviations  of  Ifj  +  (  from  perfect 
regularity  within  the  A 15  compounds  are  generally 


5%  or  less,  whereas  sin“(6j +  and  rjrj^t  may 
vary  by  more  than  25%.  Comparing  the  .415  com¬ 
pounds  with  Nb  and  V,  it  is  evident  that  Ifj  + 1  of 
the  element  is  smaller  (for  /  =  1,  by  nearly  a  factor 
of  2)  than  that  of  the  .415  compound.  However, 
for  the  dominant  1=2  channel,  variations  in 
sin*(S(  — SJ  +  |)  and  rjrj+ ,  partially  cancel,  leaving 
7/‘/  + 1  as  the  most  transferable  quantitity.  The  rela¬ 
tive  constancy  of  {I2}  in  /415  compounds, 
described  in  the  form  17  <x  MEF),  has  been  pointed 
out  previously.12 

The  situation  for  the  B\  compounds  is  not  as 
clear.  As  mentioned  previously,  the  extrapolation 
to  the  radius  Rs  =  2.51  a.u.  may  be  introducing 
some  unphysical  irregularities,  so  only  a  pair  of 
compounds  is  not  sufficient  to  allow  an  evaluation 
of  constancy  of  /2/+(  within  this  class.  It  is  clear, 
however,  that  I2  3  is  significantly  larger  than  in  the 
A 15  compounds,  which  may  lead  to  comparable 
values  of  tj(Nb)  in  spite  of  a  lower  density  of 
states.  This  is  discussed  further  below  and  in  Sec. 
IV  C. 

In  addition  to  variations  in  72/  _ ,  two  possibili¬ 
ties  remain  for  producing  a  larger  Nb  or  V  contri- 
buton  to  17.  The  obvious  possibility  is  to  find  a 
compound  with  a  larger  value  of  N(EF )  per  atom. 
This  approach  leads  to  a  well-known  instability  to¬ 
ward  formation  of  a  material  with  a  smaller  value 
of  N(Ef).  Rather  than  placing  a  very  large  num¬ 
ber  of  electrons  at  the  highest  occupied  energy,  the 
compound  will  tend  to  (a)  distort  to  a  lower  sym¬ 
metry,  thereby  moving  some  electrons  to  lower  en- 


TABLE  II.  Normalization  integrals  rj  (a.u.),  sin:(6]  — 6i+i)=S2/+i,  I2J+i  (eV/A)!,  Vu~\  an d  tj1  (eV/A)2.  and  (I2)4 
(eV/A)2.  Note  the  definition  tia  =  .\(Ef)a(I2)a,  with  .V( Ef)A  defined  in  the  caption  to  Table  I. 


NbC 

0.715 

0.436 

3.40 

0.0257 

0.32 

0.95 

0.92 

675 

421 

7000 

0.02 

0.18 

4.62 

4.82 

NbN 

0.770 

0,449 

2.95 

0.0249 

0.26 

0.67 

0.98 

505 

333 

8850 

0.00 

0.06 

6.91 

6.97 

NbjAl 

0.676 

0.358 

3.54 

0.0235 

0.35 

0.95 

0.56 

966 

498 

4480 

0.55 

5.27 

6.09 

11.91 

Nb,Ga 

0.677 

0.363 

3.58 

0.0246 

0.36 

0.97 

0.61 

956 

495 

4580 

0.55 

5.48 

5.90 

11.93 

NbiSi 

0.688 

0.359 

3.73 

0.0235 

0.35 

0.97 

0.61 

929 

478 

4570 

0.12 

2.04 

2.69 

4.85 

NbjGe 

0.671 

0.366 

3.49 

0.0256 

0.36 

0.98 

0.62 

959 

508 

4600 

0.15 

2.91 

3.65 

6.71 

NbjSn 

0.694 

0.349 

3.98 

0.0242 

0.34 

0.98 

0.64 

932 

466 

4400 

0.15 

3.91 

4.44 

8.50 

NbjSb 

0.678 

0.355 

3.88 

0.0271 

0.35 

1.00 

0.72 

963 

480 

4490 

0.05 

1.22 

2.00 

3.26 

Nb 

0.819 

0.316 

7.90 

0.0176 

0.29 

1.00 

0.82 

743 

262 

3900 

0.47 

2.72 

4.42 

7.61 

V 

0.900 

0.363 

12.97 

0.0071 

0.25 

0.87 

0.73 

495 

122 

5230 

0.20 

1.64 

5.06 

6.90 

VjAl 

0.704 

0.366 

6.96 

0.0101 

0.32 

0.89 

0.66 

822 

231 

6220 

0.16 

2.47 

4.38 

7.01 

VjGa 

0.706 

0.370 

7.21 

0.0104 

0.32 

0.93 

0.72 

808 

230 

6400 

0.20 

4.09 

5.72 

10.02 

VjSi 

0.658 

0.376 

5.71 

0.0122 

0.34 

0.92 

0.69 

919 

284 

6500 

0.13 

3.32 

5.45 

8.90 

VjGe 

0.676 

0.376 

6.31 

0.0117 

0.33 

0.94 

0.73 

869 

262 

6480 

0.07 

1.34 

3.64 

5.05 

VjSn 

0.746 

0.362 

8.62 

0.0101 

0.30 

0.94 

0.77 

733 

200 

5840 

0.03 

0.94 

2.92 

3.90 
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ergy  and  lowering  ,V(  EF'<  (the  band  Jahn-Teller  ef¬ 
fect),  bi  crystallize  in  an  unrelated  structure  with  a 
smaller  value  of  .V(£f ),  or  (c)  phase  separate  into 
distinct  compounds  with  lower  values  of  ,V(£P). 

The  alternative  choice  is  to  increase  7  by  shift¬ 
ing  the  character  of  states  at  Er,  that  is,  altering 
the  fractions  /;  to  maximize  ( /  :).  An  extreme  ex¬ 
ample  of  this  effect  is  to  shift  'somehow)  the  s  and 
p  character  'or  an  equal  amount  of  d  character)  in 
Nb  to  /  character,  thereby  taking  advantage  of  the 
stronger  d-f  scattering  described  by  l\j  in  Table 
II.  and  in  this  particular  example  increasing  <  /2) 
by  a  factor  of  8.  It  is  not  immediately  clear  how 
such  a  shift  is  to  be  accomplished,  although  a 
change  of  crystal  structure  and/or  bonding  charac¬ 
ter  is  probably  required. 

A  possible  example  of  this  effect  may  be  provid¬ 
ed  bv  the  refractory  compounds  NbC  and  NbN, 
which  are  currently  understood  in  terms  of  a  much 
more  covalent  type  of  bonding14  20  than  the  pri¬ 
marily  metallic  bcc  elements  and  115  compounds. 
As  a  result  /j  is  twice  as  large,  and  much 
smaller  for  Nb  in  these  compounds  than  in  the 
other  in  Table  I.  Since  Nb  has  no  atomic  /states, 
f  i  results  primarily  from  the  Nb-site  decomposi¬ 
tion  of  neighboring  atom  states/13  and  covalent 
bonding  <i.e.,  strong  overlap  extending  into  the  Nb 
sphere)  with  C  or  Sp  states  evidently  enhances  the 
Nb  /=  3  character  of  states  of  EF  in  these  com¬ 
pounds.  Of  course,  ( Iz )  will  be  maximized  if  as 
much  weight  as  possible  is  shifted  into  the  channel 
r with  maximum  Ifj+l  such  that  //-=/jr+1,  with 
other  fractions  ft  vanishing.  However,  such  shifts 
may  well  lead  to  small  values  of  .V(£f ),  or  if  not, 
the  resulting  increase  in  the  electron-phonon  in¬ 
teraction  can  lead  to  a  “covalent  instability”  as  dis¬ 
cussed  previously  .21,22 


IV.  DISCUSSION  AND  CONCLUSIONS 
A.  Elemental  transition  metals 

Expressing  7  =  .V(£f)</2)  with  {I1)  given  by 
Eq.  (8)  leads  to  a  simple  understanding  of  the  vari¬ 
ation  of  Te  in  transition  metals  and  their  alloys. 
Both  ifj + 1  and  //  vary  smoothly  with  Fermi-level 
position  (equivalent  to  the  mean  valence  or  alloy 
concentration  in  the  rigid-band  model),  so  rapid 
changes  in  Te  such  as  occur  in  the  NbMo  alloy 
system23  result  from  variations  in  .V(£f),  with 
some  contribution  due  to  changes  in  («2).  On  the 
larger  scale,  however,  Butler4  has  found  (I1)  to 


vary  by  a  factor  of  40  within  the  4 d  transition 
series,  in  large  part  due  to  the  factor/;.  On  this 
point  there  is  little  to  add  to  his  thorough  discus¬ 
sion  except  to  recall  the  interesting  crystal- 
structure  effects  found  in  Ru.  Fur  Ru  (which  is 
actually  hep)  Butler  carrier  out  both  bcc  and  fee 
calculations,  finding  /;  and  (I2)  to  be  about  20% 
larger  in  the  fee  phase.  However,  .V  £f)  is  small¬ 
er  by  33%,  leading  to  a  smaller  value  of  7  in  the 
fee  phase.  The  differences  in  /_•  between  the  two 
phases  are  no  doubt  related  to  the  nearest-neighbor 
distances  and  coordination  numbers,  but  such  de¬ 
tailed  relationships  are  not  well  understood  at 
present. 

The  present  picture  provides  a  new  viewpoint  on 
the  differences  in  the  /  =  3  effects  between  V  and 
Nb.  In  the  important  df  channel  is  one-third 
larger  in  V  than  in  Nb,  with  the  smallness  of  ri  in 
V  more  than  compensating  for  the  less  favorable 
values  of  rt  and  sin:(S;  —  5;).  The  average  matrix 
element  (I2),  however,  is  only  two-thirds  that  of 
Nb.  Although  this  is  due  partly  to  large  values  of 
/0.1  and  /2  -  in  Nb,  it  is  the  larger  value  of  /3  in 
Nb  which  is  primarily  responsible.  This  illustrates 
that  the  presence  of  a  larger  /3  will  be  accom¬ 
panied  by  a  larger  r3,  and  hence  smaller  /tj  and 
that  the  relative  importance  of  these  effects  is  not 
necessarily  reflected  in  the  crystalline  enhancement 
Vj  (which  is  nearly  twice  as  large  in  V  as  in  Nb). 


B.  A 15  compounds 

It  can  be  noted  from  Table  II  that  within  each 
of  the  NbyY  and  VyX  classes,  the  values  of  /2/  +  l 
are  virtually  constant.  The  NbjY  class  includes  X 
atoms  with  valences  of  3.  4,  and  5  leading  to 
differing  positions  of  the  Fermi  level  and  varia¬ 
tions  of  .V(£f  >  by  a  factor  of  4.  NbjSn  and 
NbjSb.  with  Tc's  which  differ  by  2  orders  of  mag¬ 
nitude  (18.2  and  0.2  K,  respectively),  illustrate 
dramatically  how  their  differences  in  7  arise  solely 
from  the  differences  in  .V(£f).  More  to  the  point 
of  this  study,  however,  is  the  degree  of  regularity 
of  lh  + 1  and  //  which  is  not  apparent  in  the  v,  nor 
always  in  the  factors  of  sin:(5;  —  8|_|)  (see  Tables  I 
and  II). 

It  has  been  established  in  several  A 1 5  com¬ 
pounds24  that  Tc  is  sensitive  to  the  degree  of 
disorder  (as  measured  by  the  residual  resistivity). 
The  main  effect  of  disorder  in  the  low-disorder  re¬ 
gime  is  to  broaden  the  electron  states  (“lifetime  ef¬ 
fects").  The  effect  on  7  and  on  Te  if  the  electron- 


23 


752 


WARREN  E.  PICKETT 


25 


phonon  spectral  function  a2F  is  not  altered  by  the 
disorder,  can  be  accounted  for  by  using 
Lorentzian-broadened  values  of  N(EF)  and  (I1). 
As  was  pointed  out  above,  however,  </2(£)>  is  a 
smooth  function  which  will  be  insensitive  to 
broadening.  On  the  other  hand,  many  unusual 
properties  of  ,415  compounds  have  been  interpreted 
in  terms  of  sharp  structure  in  ,V(£)  and  such 
structure  has  been  verified  by  band-structure  calcu¬ 
lations.12-25  The  present  interpretation  of  the 
RMTA  makes  it  clear  that  the  extreme  sensitivity 
of  Tc  to  disorder  must  be  reflecting  the  fine  struc¬ 
ture. 


C.  NbC  AND  NbN 

These  SI  (NaCl-structure)  compounds  differ 
from  the  .415  compounds  in  having  a  more  strong¬ 
ly  covalent  (metal  d  state  with  nonmetal  p  state) 
and  ionic  (metal-to-nonmetal  charge  transfer  of  the 
order  of  one  electron)  bonding,  rather  than  pri¬ 
marily  metallic  bonding.  This  difference  is  reflect¬ 
ed  in  the  Fermi  energy  falling  nearer  the  d  reso¬ 
nance  (52  =  ir/2)  and  results  in  a  value  of  /2j  that 
is  80—  90%  larger  than  in  Nb  and  the  Nb-based 
,415  compounds.  Another  apparent  consequence 
(in  this  case)  of  the  strong  bonding  is  the  low  value 
of  NiEp),  which  leads  to  a  rather  unimpressive 
value  of  jjN(>  in  spite  of  (I2)  being  40  —  50% 
larger  than  in  the  A 15  compounds.  These  com¬ 
pounds  do  illustrate,  however,  that  a  change  in 
bonding  character  can  lead  to  substantial  increases 
in  (I2)  relative  to  those  in  A 15  compounds,  and 
there  remains  the  possibility  that  compounds  with 
larger  values  of  77  may  be  found. 

The  predominance  of  Nb-based  compounds  in 
high  Te  materials  is  still  not  completely  under¬ 
stood.  It  is  clearly  not  due  to  </2)\b,  since  </2) 
peaks  strongly  at  Mo  rather  than  at  Nb  in  the  4 d 
series.4  The  evidence  suggests25  this  predominance 
is  due  instead  to  the  tendency  of  Nb-based  tand 
similar)  compounds  to  form  phases  with  strongly 
split  bonding  and  antibonding  d  bands.  For  Nb, 
with  its  slightly  less  than  half-filled  d  shell,  the 
Fermi  level  is  left  in  a  favorable  region  for  large 
values  of  ,V(£f ).  The  half-filled  d  shell  Mo  in¬ 
stead  leaves  EF  in  the  low  density-of-states  “gap." 

The  bcc  transition-metal  alloys  and  the  A 1 5 
compounds  are  prime  examples  of  this  behavior. 
This  point  of  view  suggests  that  Mo  will  only  be 
useful  in  raising  Te  if  a  d  electron  is  transferred 
onto  another  atom,  transforming  Mo  into 


“pseudo-Nb.”  Exactly  this  behavior  occurs26  in 
the  Chevrel  phase  compounds  .VfMo6S8  and 
.VfMOftSeg  IM = metal  atom,  e.g.,  Pb,  Sn,...),  for 
which  d  bonding  between  neighboring  Mo  atoms  is 
strong.  An  electron  is  transferred  from  each  Mo 
to  the  chalcogen  atoms  and  EF  falls  in  a  region  of 
high  density  of  states  just  below  the  bonding¬ 
antibonding  gap.  For  PbMo6S8,  Tc  =  15  K  results. 

This  behavior  is  violated  in  Mo-chalcogen  com¬ 
pounds  in  which  d  —  d  bonding  is  less  dominant, 
and  such  compounds  tend  toward  semimetallic  or 
semiconducting  character.  It  is  also  violated  in  the 
SI -structure  compounds  such  as  those  discussed 
above,  where  metal  d,  nonmetal  p  bonding  dom¬ 
inates  and  no  d  —d  bonding-antibonding  gap  oc¬ 
curs.  Approximately  one  electron  is  transferred11 
off  the  metal  atom  in  the  Nb-based  compounds, 
and  similar  behavior  is  expected  in  their  Mo-based 
counterpans,  since  Nb  and  Mo  are  expected  to 
form  good  rigid-band  systems.  A  rigid-band  pic¬ 
ture  suggests  ,V(£f)  in  MoC  should  be  10% 
larger' 1  than  in  NbC,  and  the  results  of  Butler4 
suggest  (I2) Mo  will  be  somewhat  larger  than 
(f  ‘)>n>  (as  in  the  elements).  Indeed,  it  is  found 
that  Tc  is  30%  larger  in  MoC  (14.3  K  versus  —  1  1 
K). 

The  same  argument  suggests  Tc  of  MoN  should 
be  considerably  larger  that*  that  of  NbN  ‘Tc  =  16 
K).  Rigid-band  behavy*  .  cj-jgests  a  20% 
increase911*1  in  S(Ef)\  however,  </2)Mo  may  be 
slightly  less  in  MoN  since  EF  already  lies  slightly 
above  the  d  resonance  in  NbN  (see  Table  I).  The 
few  studies  of  molybdenum  nitrides  reported  in  the 
literature27  have  not  unambiguously  established  the 
existence  of  B\ -structure  MoN.  Although  negative 
results  often  go  unreported,  there  are  at  least  two 
reasons  why  an  extensive  search  for  this  material 
may  have  not  been  undertaken.  The  first  is  the  ap¬ 
parent  dominance  of  Nb-based  compounds  in  high 
Te  superconductors,  as  mentioned  above,  wljich 
makes  the  substitution  of  Mo  for  Nb  in  these  com¬ 
pounds  unappealing.  The  second  reason  is  the 
electron-atom  ratio  of  MoN  (e/a  =  5.5),  which 
violates  he  "Matthias  rule"28  that  high  Te  materi¬ 
als  cluster  round  e  /a  =  4.75  and  6.5  with  a  deep 
intermediate  valley.  The  Matthias  rule  can  be  un¬ 
derstood  in  terms  of  the  structure  in  .V(£)  in  ma¬ 
terials  dominated  by  d — d  bonding,  and  it  is  of 
considerable  interest  to  establish  whether  this 
guideline  is  violated  in  the  B\  structure.  Finally, 
the  lack  of  stability  of  MoN  itself  is  suggestive  of 
a  strong  electron-phonon  interaction  in  this  com¬ 
pound.  The  phonon  spectrum  will  also  affect  the 
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value  of  Te  but  such  considerations  are  beyond  the 
scope  of  the  present  paper. 

D.  A  caveat 

Allen  and  Dynes"  have  noted  that,  although 
soft  lattices  may  contribute  somewhat  to  high 
values  of  Tc,  known  metals  with  high  Tc  achieve 
this  primarily  through  a  large  value  of  77.  Within 
the  RMTA,  higher  values  of  tj  in  transition-metal 
compounds  thus  seem  to  rely  on  a  larger  SiEF)  or 
/ 2.3  or  on  larger  /3  ratios.  The  /= 3  character 
arises  from  tails  of  states  on  neighboring  atoms, 
and  an  /=3  expansion  of  these  tails  requires  (i) 
they  behave  approximately  as  y'3(*r)  as  seen  from  a 
neighboring  atom,  as  pointed  out  previously  by 
Butler/  and  !ii)  they  are  expandable  in  /= 3  angu¬ 
lar  functions.  These  requirements  might  be  expect¬ 
ed  to  point  the  way  to  crystal  structures  and  chem¬ 
ical  configurations  with  larger  values  of  /3,  and 
perhaps  17.  The  NaCl-structure  compounds,  it 
should  be  noted,  have  Nb  / 3  ratios  twice  at  large 
as  for  the  .415  compounds,  much  of  which  prob¬ 
ably  derives  from  C  or  N  p  states. 

However,  the  contribution  to  77  involving  /3 
arises  in  real  space  from  near  the  sphere  boundary, 
where  the  RMTA  model  of  the  screened  potential 
is  most  approximate.  Indeed  it  has  been  noted3 
that  Tc's  of  elemental  metals  seem  to  show  better 
agreement  with  measured  values  if  77;, 3  is  (some¬ 
what  arbitrarily)  divided  by  two/9  It  is  at  first 
glance  surprising  that  this  df  contribution,  from 
near  Rs,  is  not  sensitive  to  the  value  of  Rs,  al¬ 
though  the  present  treatment  of  RMTA  clarifies 


this  point  (neither  /3  nor  l\ }  is  sensitive  to  Rs). 
We  want  merely  to  caution  here  that  requirements 
(i)  and  (ii)  above  are  necessary  to  maximize /3 
within  RMTA,  but  that  the  RMTA  is  least  certain 
here.  It  seems  likely  that  the  inclusion  of  the  non- 
rigid  potential30  and  the  concomitant  relaxing  of 
the  /  •*=*/  +  1  selection  rule  in  a  more  rigorous 
theory  may  be  necessary  to  lead  to  a  more  funda¬ 
mental  understanding  of  the  best  mechanism  by 
which  to  increase  tj. 

For  /-band  metals  with  atomic  contributions  to 
the  /  =  3  quantities,  that  is,  rare  earths  and  ac¬ 
tinides,  Butler4  has  suggested  that  the  /  bands  may 
be  useful  in  reaching  large  17  values.  Calcula¬ 
tions17  3I">:  for  La,  Ce,  and  Th  do  not  confirm 
these  expectations,  however,  and  the  discussion  of 
Sec.  II  shows  why  narrow  /  bands  will  not  lead  to 
large  rj's,  in  spite  of  extremely  large  densities  of 
states. 

Note  added  in  proof.  Calculations  of  Tc  for 
NbN  and  NbC  have  been  presented  by  W.E.  Pick¬ 
ett,  B.M.  Klein,  and  D.  A  Papaconstantopoulos 
rPhvsica  107B.  667  (1981)1. 
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A  microscopic  treatment  of  the  consequences  for  superconductivity  of  a  nonconstant  electron¬ 
ic  density  of  states  is  presented.  Generalized  Eliashberg  gap  equations  valid  for  a  varying  densi¬ 
ty  of  states  are  presented,  from  which  the  change  of  Te  with  static  or  thermal  disorder  can  be 
calculated.  The  temperature  dependence  of  the  effective  mass  is  shown  to  be  altered  by  disor¬ 
der.  Use  of  these  results  provides  a  possible  experimental  approach  for  deducing  the  energy 
variation  of  the  density  of  states  of  superconductors. 


I.  INTRODUCTION 

There  is  currently  an  active  interest  in  the  effects 
upon  the  electronic  properties  of  metals  when  the 
density  of  states  .V(£)  cannot  be  assumed  to  be  con¬ 
stant  near  the  Fermi  level  £>.  Recent  band-structure 
calculations,  especially  on  the  A 15  compounds,  con¬ 
firm  that  in  transition- metal  compounds  with  several 
atoms  per  cell,  sharp  structure  in  ,V(£),  at  least  on 
the  order  of  the  Debye  frequency  k9o~~2S  meV,  oc¬ 
curs  in  a  region  around  £>.  Since  there  is  no  general 
agreement  (nor  will  there  be  soon)  among  band- 
structure  calculations  on  this  scale,  the  general  trend 
has  been  to  attempt  to  fit  experimental  data  with 
models  of  ,V(£),  usually  with  ad  hoc  phenomenologi¬ 
cal  expressions  or  oversimplified  theories.  Since  the 
phenomenological  approach  appears  to  be  the  only 
way  at  present  to  deduce  the  structure  in  X  ( £)  on 
this  scale,  it  is  evident  that  correct  theoretical  expres¬ 
sions  are  needed  to  allow  a  proper  unfolding  of  the 
data. 

The  effect  upon  the  critical  temperature  Tc,  due  to 
static  disorder  has  been  studied  extensively  in  A 1 5 
compounds,  where  Tc  may  either  decrease  or  increase 
with  disorder.1  Explanations  of  this  phenomenon  has 
postulated  a  smearing  of  structure2  in  X(.  £) ,  gap  an¬ 
isotropy.3  changes  in  the  phonon  spectrum,4  or  the 
special  type  of  defect3  which  is  present,  although  not 
all  of  these  could  raise  Te.  Unfortunately,  there  is  a 
great  deal  of  uncertainty  in  the  literature  about  the 
form  the  energy  smearing  will  take,  with  the  only  mi¬ 
croscopic  theory  (Labbe  and  van  Reuth,  Ref.  2)  be¬ 
ing  based  on  a  special  type  of  disorder  and  an  over¬ 
simplified  model  of  an  A  15  compound. 

In  this  paper  a  more  rigorous  theoretical  treatment 
than  previously  available  is  given  for  the  supercon¬ 
ducting  properties  of  a  metal  with  a  nonconstant 
iV(£).  The  results  will  be  confined  primarily  to  the 
effect  on  Tc  and  the  effective  mass.  Theoretically, 
the  effect  upon  Te  has  been  considered  in  different 
approximations  by  Horsch  and  Rietschel6  (HR)  and 
Ho,  Cohen,  Pickett7  (HCP),  with  apparently  conflict- 
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ing  conclusions:  HR  suggested  that  the  structure  in 
•V  ( £)  is  responsible  for  the  high  Tc  found  in  many 
A  15  compounds,  while  HCP  contended  that  structure 
in  X(E)  is  ultimately  of  no  importance  in  raising  Tc. 
Here  it  is  shown  that  both  situations  are  consistent 
with  a  full  theory.  In  addition  it  is  found  that,  for 
metals  in  which  the  spectral  function  a2F  is  accurate¬ 
ly  known,  the  dependence  of  Tc  on  (nonmagnetic) 
impurity  concentration  n,  and  the  altered  temperature 
dependence  of  the  mass  enhancement  [1  J-\(  D) 
due  to  disorder  may  provide  information  on  the  vari¬ 
ation  of  X(E)  in  the  region  of  Ef.  The  change  in 
\(  T)  has  apparently  not  been  anticipated. 

II.  ELIASHBERG  GAP  EQUATIONS 

The  Eliashberg  gap  equations  assume  a  particularly 
simple  form  when  expressed  on  the  imaginary- 
frequency  axis.  The  complex  frequency -dependent 
normal  [Z(<u)l  and  anomalous  [<6(&»)1  self-energies 
collapse  to  functions  defined  at  discrete  imaginary 
frequencies  /&»„- i2trnkaT.  Unfortunately,  in  this 
representation  there  is  no  straightforward  intuitive  in¬ 
terpretation,  and  numerical  results  which  are  ob¬ 
tained  with  relative  ease  are  not  easily  analytically 
continuable  to  the  real  axis.  Since  one  aim  here  is  a 
clear  understanding  of  the  effects  of  a  varying  X(E), 
the  real-frequency  expressions  will  be  developed  and 
interpreted. 

The  irreducible  self-energy  I  -  ■+■  Ic  I,  is  to 

be  determined  self-consistently  from  contributions 
from  phonon,  electron,  and  impurity  scattering.  The 
latter  usually  is  not  treated  explicitly,  as  it  has  only 
the  effect  of  changing  the  reference-band  structure, 
i.e.,  changing  ;V(0).  The  self-energy  is  expressed  as 

£(  k,u>)«a>[l  —  Z(k,(u))ro 

+  x(  k ,  «u) tj  +  <6(  k .  o»)t]  ,  (1) 

where  T|  and  r}  are  the  Pauli  matrices,  and  r0  is  the 
2x2  identity.  The  phonon  contribution  is  given  (on 
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the  imaginary-frequency  axis)  by1 

l*(it .(«„)  --0"'  XI*rr^TjC(  * “v* 

T«* 

*  T}D(  k  —  k  —  no^)  (2) 

where  G,  Dan  the  electron.  ,-*ionon  Green's  func¬ 
tions,  respectively,  and  g  is  (he  screened  electron- 
phonon  matrix  element.  In  this  approximation  for 
I,*  it  is  assumed  that  Migdal's  theorem  is  valid,  i.e., 
that  other  diagrams  give  a  contribution  to  £,k  which 
is  smaller  than  Eq.  (2)  by  the  factor  (k9D/Ef), 
which  is  of  the  order  of  (electron  mass/ion  mass) 1/2 
in  ordinary  systems.  Although  Migdal's  analysis  no 
longer  guarantees  that  other  diagrams  always  give  a 
negligible  result  when  the  system  contains  very  nar¬ 
row  bands,  there  should  remain  an  important  regime 


where  the  following  treatment  of  Eq.  (2)  wilt  give  the 
dominant  contribution  to  I*,. 

To  convert  the  wave-vector  <  k)  dependence  to  en¬ 
ergy  ( E)  dependence,’  define  for  any  function  C(k) 

f(£)N(£)«JC(k)Slf-£r)  .  (3) 

r 

It  will  be  assumed  that  the  replacement 
2C1(k)C2(k)5(£-£r)-C,(£)C2(£).V(£)  (4) 

k 

is  valid.  This  amounts  to  ignoring  anisotropy  and 
contributions  which  are  off-diagonal  in  the 
(suppressed)  band  index,10  but  allows  us  to  retain  the 
energy  dependence,  which  is  the  point  of  interest  in 
this  paper.  This  allows  Eq.  (2)  to  be  written 


..  , „ .  ,  r  „-#(£“)  2fla!f(£,r,n)  „/rw.  , 

!*(£.,«,)  ~  e  >2J  dE*— w  J  dClj— 


where 


<V(£)W(£')a*f(£,£';ft)  —  <V(0)  £  |*r_.j2a(  k-k';n)8(£-£r)6(£'-£r) 

r  r 


(6) 


The  phonon  spectral  density  is  denoted  by  B,  and 
a2f(0,0;O)  is  the  usual  electron-phonon  spectral 
function  a2£. 

From  Eq.  (5)  it  is  evident  that  the  £  dependence  of 
I  arises  from  the  £  dependence  of  o2F  and  the 
analogous  Coulomb  and  impurity  functions,  a  fact 
noted  by  HR.  In  systems  where  N(  E)  varies  appre¬ 
ciably,  very  little  is  known  about  this  energy  depen¬ 
dence.  For  clarity  we  make  the  approximation 

a2£(£.F;n)-o2£(0.0;n)-a2£(n)  ; 

it  should  be  kept  in  mind  that  this  approximation  is 
untested.  The  full  energy  dependence  can  be  re¬ 
tained  in  the  development  which  follows,  but  doing 
so  complicates  the  expressions  considerably. 

What  must  be  evaluated  then  is 

n(o>)-~ L  f"  dE  ,V(  E)G(  E,  at)  .  (7) 

This  can  be  done  formally  by  introducing  the  con¬ 
tinuation  jV ( £ )  of  <V(  £)  off  the  real  axis 

A'(£)-J_P  dE  jV(£)/(£-£)  .  (8) 

ITT 


and  £  is  substituted  for  £  in  G.  The  contour  can  be 
removed  to  infinity,  leaving  only  the  contributions 
from  the  poles  of  G  at 

£-±Q («•).  QU)«U«Z)2-* Ml/I  . 

X  will  be  ignored,  as  its  effect  can  be  shown  to  be 
negligible  in  all  except  the  most  pathological  cases. 
The  result  is 

ij(w)-/|[(MZr0+ir,)/0]iV.((?)  +r3tf *((?))  , 

(10) 

with 

tf±(<?)-l(V(<?)  ±N(~Q))/2  .  (11) 

In  the  limit  of  a  constant  density  of  states, 

\-(Q)  —  .V(0)  and  N+(Q)  — ‘0.  In  general  ReN _  is 
the  dominant  term,  as  discussed  below.  The  other 
contributions  dm#.,  Re,V*,  and  ImjV+]  are  much 
smaller,  since  cancellation  is  present  in  their  defini¬ 
tion  and,  being  oscillatory  about  zero  as  a  function  of 
to*  through  Q(to'),  contributions  to  the  a*’  integral 
will  tend  to  cancel.  These  factors  will  be  neglected 
relative  to  N- 


.Vis  analytic  in  both  the  upper  and  lower  half-planes 
of  £  but  has  a  discontinuity  2  S(E)  along  the  real 
axis.  This  property  allows  n  to  be  written 

n(«)--~/c«f£ #(£)<?(£.«)  .  (9) 

where  the  contour  C  encircles  the  real  axis  in  the 
negative  sense  at  a  distance  5  — 0*  above  and  below. 


III.  CONSEQUENCES  OF  ENERGY  DEPENDENCE 

It  is  straightforward  to  apply  the  same  approxima¬ 
tions  to  the  Coulomb  and  impurity  (or  more  general¬ 
ly  any  static-disorder)  contributions  and  continue  the 
resulting  expressions  to  the  real  axis.1  The  equation 
for  the  gap  A  -  d>tZ  —  Apt,  A,  is  most  easily 
derived  using  Iph  -  +  ( 1  -  Z)pnA,  etc.,  with  the 
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result  at  T  -  7V  becoming 


A  ( tu )  ”  f  dot  S(at)  : —  [«_(  at.  at)f(  —at)  *_(  at,  — <o  ')f(at') ! 

•'O  f.i 


_  [  *.,(  at.  at  ) /( —  <ii  1  —  k.(  (u.  —  tu  )/( <u  )  ]  —  I/c  tanh-^- 

(u  2 

V  [  A1"*  ^  I  \u  .  \ 


+  /«■  I*"  <ifl  a2£(  (!)/»( fl)  2  .via.,)  . 

Here  at,  -at  ■*■(— )*fl  Uc  is  the  Coulomb  pseudopocential.  u»c  is  the  cutoff  frequency, 

k±(oi.  at'  >  - J* daa:F(n)in«'-r«  +  n)  ±ru  -  <u-n)]  , 

;V(0)iV(u»)  “Re  .V-[Q(<«)  ]  —  Xliili-  f“  -  ^---.(,V[£-a>+.t/(a>))  f  ,V(  £ -at  -  ,V/(  at)  \  i  . 

2  it  £•  -r-(at) 


/(n)  denotes  the  fermion  (boson)  thermal  occupa¬ 
tion  function,  T(x)-(x+i&)-\  and 
Xf(at) -tT(<u)  -  oil  1  —  Z (at)]  denotes  the  normal- 
state  self-energy  at  T  Tc.  Within  the  approxima¬ 
tions  outlined  above,  the  sole  change  in  the  gap 
equation  (12)  is  the  replacement  of  .V(0)  by 
.V(0)/V(w),  given  by  Eq.  (14).  Physically  Si  at)  is  a 
dimensionless  function  (as  defined  here)  which  re¬ 
flects  the  average  of  the  density  of  states  at  the 
quasiparticle  energies  ±[a»  — A/(<u)l,  broadened  by 
the  quasiparticle  damping  T(<u).  It  should  be  noted 
that,  for  strong  electron-phonon  coupling  and  a 
varying  ,V(£)  near  Ef,  neglecting  the  temperature 
dependence  of  the  energy  shift  Xf  may  not  be  justi¬ 
fied. 

The  first  effect  [averaging  over  ±(at  -XI)]  was 
studied  in  some  detail  by  HR,  who  neglected  the 
smearing  and  made  the  replacement 
,V(0)  —  jliVU-.l/)  +,V(-«  +  A/l.  Using  a  model 
density  of  states  with  an  E~m  singularity,  numerical 
solutions  indicated  an  enhancement  of  7V,  over  the 
value  for  constant  N( 0),  for  Fermi-level  placements 
Ef  =*  j(l  +\)<u o  relative  to  the  singularity,  while 
other  placements  gave  a  reduction  of  Tc(aio  is  the 
Debye  frequency).  This  behavior  was  interpreted  as 
a  reduction  of  the  repulsive  part  of  the  interaction 
resulting  from  the  interplay  of  the  energy  dependence 
of  S(  E)  and  the  frequency  dependence  of  the  in¬ 
teraction.  Part  of  the  effect  no  doubt  is  due  to  the 


fact  that  }[.V(oj-.Vf)  -r.V<  -at  +Xf)  ]  >  ,V(0>  for 
the  important  frequencies  for  this  placement  of  £>. 

In  a  sense  the  effective  value  of  the  density  of  states 
is  increased.  Since  the  other  parameters  are  fixed, 
this  results  in  an  increase  in  Tc.  The  treatment  of 
HR  is  valid  if  the  damping  r(u>)  is  negligible:  this 
approximation  is  dose  to  breaking  down  at  T  =  Te  in 
high-  Tc  A  15  compounds.' 

The  second  effect,  the  smearing  of  .V(£)  resulting 
from  quasiparticle  damping,  has  been  considered  by 
HCP,  who  argue  that  for  strong-coupling  (and  high 
Tc)  superconductors,  peaks  in  ,V(£)  become  ineffec¬ 
tive  in  raising  Tc.  A  large  A  leads  to  a  large  damping 
T(at)  and  thus  a  greater  smearing  of  .V(  £).  It  was 
argued  that  for  large  enough  A,  the  peaks  in  ,V(  £) 
will  be  smeared  out  at  T  -  Tc ,  implying  that  large 
values  of  .V(0)  (resulting  from  sharp  structure  in  the 
T-0  band  structure)  are  ultimately  useless  in  reach¬ 
ing  very  high  Tc.  Roughly  speaking,  the  broadening 
effect,  with  peaks  heights  decreasing  in  proportion  to 
A,  will  dominate  over  the  increase  in  Tc,  (ultimately) 
proportional11  to  Al/:.  By  assuming  that  Ef  lies  near 
the  center  of  a  very  narrow  peak  in  S(E),  it  was  es¬ 
timated  that  Tc  could  be  reduced  by  —2-4  K  by  this 
effect  in  high- 7V  (—25  K)  .415  compounds.  The 
averaging  of  S(E)  over  a  range  ±atD  was  neglected 
by  HCP,  but  in  fact  tends  to  strengthen  their  conclu¬ 
sions  if  Er  lies  precisely  at  a  peak  in  S(E). 

The  normal  self-energy  at  T^TC  is 


XT  (at)  —ir(at)  —  dot'  S(at')  (#c-(u,  «’)/(—  u>’)  +  *-(«.  -ot')f(ai')  1 

-in  d{la1F((l)n(nHS(ut  +  (l)  +JV{«- n)]  +(1  -Z  ),«-/*, iV(0)|  Kf.VU)  .  (15) 


The  Coulomb  contribution  ( 1  -  Z)e«  arises  from 
high-frequency  processes  and  is  unchanged  by  struc¬ 
ture  in  .V(£);  the  impurity  contribution  is  obtained 
by  the  replacement  S(0)  —S(0)S(ut)  and  a  small 
real  part  has  been  ignored.  The  impurities,  with 
(constant)  matrix  elements  I  F'|J,  have  been  treated 


in  the  second  Born  approximation11  and  in  the  small 
impurity  concentration  limit  n,  «  1.  It  is  this  self¬ 
energy  which  must  be  inserted  in  Eqs.  (14)  and  (12) 
when  solving  for  Tc.  The  impurities  do  not  enter  ex¬ 
plicitly  inio  the  gap  Eq.  (12),  a  manifestation  of 
Anderson’s  theorem.13  [This  appears  not  to  be  true  if 
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the  energy  dependence  of  a2F(  £.£*;  0 )  is  important, 
in  which  case  the  gap  is  renormalized  by  the  impuri¬ 
ties!]  The  impurities,  however,  do  add  to  the  damp¬ 
ing,  and  hence  contribute  to  the  smearing  of  .V(  E)  in 
Eq.  (12).  This  implies  a  variation  of  Te  with  impurity 
content  which  in  principle  can  be  used  to  provide  in¬ 
formation  on  the  variation  of  /V(  £) . 

The  temperature  dependence  of  the  effective  mass 
( m*fm )  -  1  +  X(  T)  is  also  altered  by  a  variation  in 
N(E).  The  phonon  contribution  is  given  by 

X(D-2  do,N(<*) 


X 


£  dCla:F(n) 


A-  a.)  fim) 

(ft+cu)2  (fl-eu)2 


a  less  than  obvious  generalization  of  the  usual  ex¬ 
pression.14  This  relation  generalizes  the  known 
result11  that  strong  impurity  scattering  does  not  des¬ 
troy  the  electron-phonon  mass  enhancement.  Impur¬ 
ity  scattering  does  enter  into  \  and  may  lead  to  either 
an  increase  or  decrease  in  X. 


IV.  DISCUSSION 

In  previous  treatments2-7  a  dc  relaxation  time  r, 
usually  taken  from  the  resistivity,  together  with  the 
uncertainty  relation  A£  -»/r,  has  been  used  to 
smear  an  assumed  N(E).  The  resulting  effective 
density  of  states  N^iEr),  together  with  the  assump¬ 
tion  Ax  «  \N(  £f),  was  used  to  estimate  changes  in 
Tc.  The  form  of  Eq.  (12),  as  well  as  the  analysis  of 
HR,  suggests  that  this  procedure  will  give  a  reason¬ 
able  approximation  only  if  the  scale  of  the  important 
structure  in  N(E)  j$  much  larger  than  <#D ,  as  in  the 
case  of  large  disorder  (concentrated  alloys,  large  im¬ 
purity  concentrations,  or  radiation  damage).  The 
present  treatment  of  course  is  only  strictly  valid  for 
small  amounts  of  disorder. 

A  quantitative  study  of  the  usefulness  of  these 
results  in  empirically  inferring  structure  in  N(E)  will 
be  discussed  elsewhere.  A  few  observations  should 
be  made  here.  Superficially  Eq.  (16)  appears  to  be 
the  more  promising  approach,  since  the  relation 
between  .V  ( £)  and  the  measured  quantity  is  more 
direct,  and  X(  D  can  be  extended  (in  principle)  to 
temperatures  T  >  Te  where  the  energy  smearing  is 
larger.  Although  light  effective  masses  have  been 
seen  up  to  T  — 100  K  in  de  Haas-van  Alphen  oscilla¬ 
tions.'4  the  heavy  masses  which  should  be  associated 
with  structure  in  N(E)  most  likely  will  be  observable 


only  over  a  much  more  limited  range  of  T.  Another 
point  to  consider  is  whether  a  given  experiment  will 
actually  measure  the  T-dependent  mass  enhancement 
as  given  by  Eq.  (16).  For  example,  the  usual  T 
dependence14  of  X  [i.e.,  when  ,V(£)  can  be  con¬ 
sidered  constant]  is  measured  by  cyclotron  reso¬ 
nance17  but  apparently  not  by  the  de  Haas-van  Al¬ 
phen  effect.17-  '*  More  theoretical  work  is  necessary 
to  determine  the  appropriate  experiments  to  measure 
X(  T)  of  Eq.  (16). 

In  the  case  of  the  dependence  of  Tc  on  disorder, 
the  difficulty  is  less  with  the  observation  than  with 
the  interpretation  of  the  effect.  The  variation  of  Tc 
with  large  amounts  of  static  disorder  has  been  known 
for  some  time.”  However,  for  large  static  disorder 
the  variation  of  a:F(  ft  ),M  as  well  as  the  variation  of 
N(to),  is  important  in  Eq.  (12)  which  determines  7V 
Of  course  a2F  should  be  determined  by  inverting  Eq. 
(12)  rather  than  the  usual,  simpler  N(E)  —N( 0) 
equation;  there  is  a  difference  even  at  T  -0. 

A  different  approach  to  these  equations  has  been 
taken  by  Lie  and  Carbotte,21  who  have  evaluated 
8  Te/bN(,E)  for  a  number  of  superconductors.  Their 
results  allow  an  estimation  of  the  change  in  Te  due  to 
a  change  in  N(  £)  more  directly  than  the  procedure 
outlined  in  this  paper.  However,  since  their  results 
are  numerical,  so  far  they  have  not  provided  much 
insight  into  the  physical  mechanisms  which  affect  Te. 

In  this  paper  a  microscopic  formulation  of  the 
manner  in  which  a  varying  density  of  states  affects 
the  superconducting  properties  of  a  metal  has  been 
developed  and  discussed.  It  is  shown  that  the  tem¬ 
perature  dependence  of  the  effective  mass  reflects 
the  energy  variation,  as  does  the  Tc  dependence  on 
impurity  concentration.  Finally,  it  is  suggested  that 
the  energy  spectrum  near  the  Fermi  level  may  be  in¬ 
ferred  from  a  careful  inversion  of  experimental  data 
on  rt(n,)  and  X(  T). 

Note  added  in  proof.  Although  the  chemical  potential 
has  not  been  exhibited  explicitly  in  this  paper,  its 
variation  with  temperature  can  be  important  and 
must  be  inducted  in  Eqs.  (12) -(IS).  This  point 
often  has  been  neglected  in  previous  work. 
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INFLUENCE  OF  ELECTRONIC  STRUCTURE  ON  SUPERCONDUCTING  PROPERTIES 
OF  COMPLEX  CRYSTALS:  THEORY  AND  APPLICATION  TO  NbjSn 

W.  E.  Pickett  and  B.  M.  Klein 

Naval  Research  Laboratory. 

Washington.  D.C.  20375 


A  thermodynamic  theory  valid  for  complex  crystalline  superconductors  is  applied  to  NbjSn.  It  is  shown  how  the 
structure  of  the  Eliashberg  equation  and  the  solution  is  altered  by  density  of  slates  fine  structure.  One  important  impli¬ 
cation  is  that  the  current  method  of  inversion  of  tunneling  data  can  be  only  approximate. 


1.  Introduction 

The  anomalous  temperature  (T)  dependence  of  normal 
stale  properties  of  A15  compounds1  have  often  been  supposed 
to  arise  from  a  variation  in  the  electronic  state  density  N(E) 
near  the  Fermi  energy  Ef  on  the  scale  of  the  maximum  phonon 
frequency  fl .  That  such  structure  should  occur  in  crystals  con¬ 
taining  several  transition-metal  atoms  per  cell  can  be  argued 
quite  generally/  and  recently  several  band  structure  studies  1  of 
A15  compounds  have  verified  structure  on  this  scale.  Unfor¬ 
tunately  the  usual  theory  of  strongly  coupled  electron-phonon 
systems  is  based  on  a  constant  density  of  states  (CDOS)  as¬ 
sumption  and  is  not  directly  applicable  in  the  more  general 
case. 

Since  the  early  weak-coupling  models  of  the  critical  tem¬ 
perature  Tr  due  to  Barisic  and  collaborators4,  only  a  few  at¬ 
tempts5  at  a  general  strong  coupling  theory  have  been  made. 
However,  each  of  these  makes  unnecessary  approximations 
and/or  uses  pathological  density  of  states  models  which  obscure 
the  underlying  principles.  Recently  one  of  the  authors6  has 
generalized  Eliashberg  strong-coupling  theory,  subject  only  to 
the  isotropic  approximation.  In  this  paper  we  report  the  applica¬ 
tion  of  this  theory  to  NbjSn.  We  point  out  several  novel  impli¬ 
cations  of  the  theory,  emphasizing  that  superconducting  proper¬ 
ties,  and  in  particular  the  gap  function  A(u).  will  be  affected 
significantly  by  even  non-pathological  structure  in  N(E). 

2.  Formalism 

The  electron  self-energy  Z„  -  Z  (*,/<»,)  in  the  Nambu- 
Matsubara  formalism  is  given  in  standard  notation1  by 

Z.  “  i»,(l  —  Z„)  +  x.r (1) 

In  Eliashberg’s  approach  Z  is  determined  self-consisienlly  in 
terms  of  the  renormalized  Green’s  function,  given  by 

G;'  -  i«,  -  (£,  -Or,-  Z„  (2) 

where  {  is  the  chemical  potential  and  «,  «  (2«f  l)irk*r.  In 
the  isotropic  approximation  the  wavevector  (k  )  dependence  of 
Z  is  averaged  over6,  leaving  only  a  weak  dependence  of  Z( E) 
on  the  energy  E  which  can  be  disregarded.  For  a  general  state 
density  the  generalized  gap  equation  at  7*r  can  be  written 

v'  I- 

£  |A  (n  —  m)-t-A(»+-m)  (3) 

-2ii*-»»,(2n  +  D-^-j -  0, 


lv  - 

Z,  “  1  +  - - -  T  |A  In— m)  —  A  (a  +  m)|A/„,  (4) 

2n  +  1 

X.m  — irkg T£  (A  (n  -  m)  +  A  (n  +  m )  —  hi  *!/*,.  (5) 


The  state  density  enters  through  the  even  and  odd  averages 


ZJ-.I 
-c  + 


SiE)dE 

(£  -  {  +  x.)2  +  <*;Z.2 


(6) 


It  is  notable  that  the  normally  ubiquitous  factor  .V  ( £ri  nowhere 
appears  in  the  theory.  Accordingly  the  COOS  quantities  A  and 
/»'.  which  are  formally  proportional  to  S(Ef).  are  replaced  by 
the  corresponding  intensive  quantities  X  and  ii  "  which  depend  on 
the  character  of  the  wavef unctions,  but  not  on  the  density  of 
states,  in  the  vicinity  of  £>.  In  the  limit  of  a  slowly  varying 
.V(£),  AM£>)A  — A  and  > •— p.  •  in  the  equations  above. 
a  0)  is  given  by 

A  0)  -  2/ toa!F(u)da/(u3  +  (2ir;T]J)  (7) 

in  terms  of  an  intensive  electron-phonon  spectral  function  oJF. 

Equations  (3)-(6)  and  the  relation*  determining  {( T)  form 
a  system  of  coupled  non-linear  matrix  equations  which  must  be 
solved  to  find  Tc  and  the  self-energies.  In  the  CDOS  limit  V. 
approaches  NIEr ).  P,  vanishes  and  the  gap  equation  reduces  to 
the  well-known  form.*  We  show  below,  however,  that  for  a 
realistic  N(E),  .V,  and  P,  differ  greatly  from  these  limits  for 
NbjSn  (and  no  doubt  for  other  compounds).  Before  discussing 
the  numerical  results  it  is  helpful  to  understand  the  implications 
of  this  new  behavior. 

The  eigenvector  components  of  the  gap  Eq.  (3)  are  given 
by  .V,  A„,  rather  than  by  A„  as  in  the  CDOS  theory.  The  ker¬ 
nel  is  altered  from  the  usual  theory  only  by  the  changes  in  the 
diagonal  term  Z„/.V,,  with  changes  in  Z,  tending  to  be  can¬ 
celled  by  the  denominator.  If  for  the  moment  the  correction  to 
the  kernel  is  ignored.  T  and  the  corresponding  eigenvector  is 
unchanged  but  the  gap  function  A,  is  altered  by  a  fraction 
\/N,.  This  leads  to  the  conclusion  that  the  gap  function  Alai), 
which  is  given  by  the  analytic  continuation  of  A,  to  the  real 
axis,  is  fundamentally  altered  by  a  non-constant  .V(£)  It  then 
becomes  evident  that  A(w)  and  hence  a:F(w)  will  he  different 
from  that  inferred  from  tunneling  data  using  the  CDOS  inversion 
procedure.  It  will  be  shown  below  that  Tr  in  fact  is  altered 
appreciably  by  the  corrections  to  the  kernel  diagonal. 


3.  Application  to  Nb,  Sn 


where  A.  -  A„  ~  d„/Z„  and  Nr  is  the  cutoff  integer.  We  have  solved  the  coupled  equations’0  for  Nb,$n  using 

The  normal  self-energies  «„Z„  and  x.  are  given  by  a 2F  and  p  *  deduced  by  Wolf  et  at.  "  and  the  calculated  S(E) 
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Fig.  1  —  (a)  The  experimentally  deduced  specirai  function  a1  Flat)  of  Wolf  et  at.  for  NbjSn  (Ref.  II)  for  NbjSn 
and  the  calculated  function  ,V(£)  of  Klein  ei  at.  (Ref.  3)  on  the  same  scale  as  a:F.  (b)  .V(£)  shown  on  a  larger 
scale  to  display  its  fine  structure. 


of  Klein  ef  at. 3  We  make  the  somewhat  arbitrary  identification 
d!F(u)  -  a:F  ( u)/N(Ef ),  and  similarly  for*!*,  as  the  most 
reasonable  choice  available.  The  functions  a:F  and  .V(£)  (the 
latter  interpolated  linearly  between  values  on  a  1  mRy  grid)  are 
shown  in  Figure  1  on  the  same  scale  to  emphasize  that, 
although  N(E)  has  fine  structure,  it  is  utterly  non-pathalogical, 
in  contrast  to  those  used  previously3  which  had  singularities 
and/or  discontinuities.  The  general  behavior  we  discuss  is  not 
critically  dependent  on  the  chosen  functions a2FM  and  . V ( £ ) . 
[It  should  be  noted  that  we  find  that  accurate  numerical  results 
cannot  always  be  obtained  if  C  and  x  are  disregarded.  U  is  not 
appreciated  generally  that  { ( D  is  strongly  influenced  by  the 
electron-phonon  interaction.] 

In  Figure  2a  we  show  N,  and  P,  for  NbjSn  at  Tc  «  IS  K. 
Clearly  N differs  greatly  from  N(Ef)  and  /*,  is  strongly  non¬ 
zero.  Both  vary  rapidly  for  small  n.  reflecting  the  fine  structure 
near  Ef.  before  tending  to  saturate  at  larger  values  of  n.  As  dis¬ 
cussed  above,  the  variation  of  .V,  leads  to  corrections  to  A, 
relative  to  the  constant  V,  approximation  A,’.  In  Figure  2b  this 
difference  A„  -  A”  is  displayed.  The  normalization 
A,  -  l  —  A*  is  used  since  the  magnitude  of  A  is  arbitrary 
(infinitesmal)  at  Tc,  but  below  T,  the  magnitude  as  well  as  the 
frequency  dependence  of  A  will  be  altered.  To  learn  the  effect  of 
ME)  on  the  complex  gap  function  A(«)  the  solution  |A«]  of 
the  (appropriately  generalized3)  equations  at  T  <  Tt  can  be 
analytically  continued  to  the  real  axis.  Unfortunately  we  have 
found  the  Pad<  approximant  approach11  to  be  inadequate  at 
T,~  IS  AT  due  to  having  only  3-4  Matsu bara  frequencies  in  the 
region  <*,<(!  where  A(u)  has  sharp  structure.  However, 
Vidberg  and  Serene11  have  shown  that  this  continuation  to  the 
real  axis  can  be  performed  reliably  from  the  first  100-200 
values  of  A,  at  low  temperature.  Since  the  shape  of  Afw.D  is 
nearly  T-indepcndent,  the  sensitivity  of  A(«)  upon  the  values 


!A„)  found  by  Vidberg  and  Serene  indicates  important  correc¬ 
tions  to  A(i»)  will  result.  Numerical  solutions  to  Eqs.  (3)-(6) 
at  T  «  Tc  will  be  needed  to  study  this  effect  more  quantita¬ 
tively. 

Figure  2  indicates  that  corrections  even  from  non- 
pathological  DOS  functions  can  be  —  5*H>  of  A,  %  Alw  -  0). 
However,  to  give  a  reliable  picture  of  the  structure  in  o  :F  it  is 
necessary"  to  invert  tunneling  data  to  a  relative  accuracy  of 
—  1(TJ.  Here  it  is  noteworthy  that  this  has  not  been  achieved 
in  A1S  compounds,  although  heretofore  this  deficiency  has 
been  ascribed  solely  to  difficulties  in  preparing  ideal  junctions. 

To  illustrate  the  effect  of  DOS  structure  upon  Tr  we  allow 
the  assumed  Fermi  level  Ef  to  vary  over  a  range  of  5  mRy  near 
the  calculated  value  £f  „  of  Klein  et  al.  In  Figure  3  it  can  be 
seen  that  Tc  resulting  from  the  CDOS  theory,  denoted  77. 
varies  from  10  to  26  K  within  3  mRy,  whereas  Tr  from  Eq.  (3) 
varies  only  from  13  to  18  K.  Not  least  among  the  implications 
of  this  drastic  smoothing  of  Te  is  that  band  theory  need  not  be 
absolutely  accurate  on  a  1  mRy  scale  to  provide  reasonable 
values  of  Tr;  conversely,  Tc  contains  less  straightforward  infor¬ 
mation  about  .V(£)  than  usually  assumed. 

Currently  there  is  considerable  interest  in  obuining 
approximate  Tc  values  by  using  an  appropriate  average  over 
NiE)  in  the  CDOS  theory  in  the  spirit  of  Nettel  and  Thomas.5 
Equation  (3)  suggests  an  average  St/f  given  by 

(») 

*  1 

The  resulting  critical  temperature  TJrf  also  is  shown  in  Figure 
3.  Although  77//  is  too  large  by  a  nearly  constant  amount 
(—  l.SK),  it  reflects  the  trends  in  Tc  quite  well.  Improvements 
upon  this  result  will  be  investigated  in  a  subsequent  paper. 
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Fig.  2  —  The  even  and  odd  density  of  states  averages  .V,  and  P,  respectively  ffor  the  calculated  value 
Tc  —  16.6  K  for  Nb}Sn)  plotted  versus  Matsubara  frequency  and  the  correction  to  the  gap  function  A,  due  to 
the  density  of  states  shown  in  Figure  1.  Note  the  change  of  scale  at  n  -  20.  above  which  only  every  tenth  value 
is  shown. 
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Fig.  3  —  The  behavior  of  the  critical  temperature  versus  the  assumed  position  of  Er.  Results  of  the  general 
theory  are  denoted  by  T-.  those  of  the  CDOS  theory  by  77  and  those  of  the  average  value  of  ,V(£)  suggested  in 
the  text  by  TJ,f. 
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4.  Discussion 

We  have  shown  that  when  fine  structure  in  the  electronic 
spectrum  exists  near  Ef  it  is  necessary  to  take  it  into  account  in 
the  inversion  of  the  tunneling  data  to  obtain  accurate  values  of 
A  and  a2F as  well  as  Tc.  Our  use  in  the  calculations  of  a ‘Flu) 
<*  a2F(u),  with  the  latter  deduced  from  the  CDOS  theory,  was 
an  expediency  which  will  not  be  satisfied  by  a  self-consistently 
determined  <> 2  F. 

It  should  be  emphasized  that  the  present  calculations  have 
disregarded  the  disorder  which  is  present  in  all  NbjSn  samples. 
It  is  well  known  that  disorder  broadening  of  the  DOS  will  tend 
to  smooth  fine  structure  and  thereby  reduce  the  effects  dis¬ 
cussed  in  this  paper.  The  theory  presented  here  recently  has 
been  extended  to  include  the  effects  of  disorder.  Preliminary 


work  indicates  that  for  a  resistivity  ratio  pl300  KJ/p(20  K)  ™ 
10  (a  'good"  NbjSn  sampler,  ihe  correction  (A,  -  A„')/A.,'  in 
Figure  2  is  reduced  by  roughly  a  factor  of  2.  This  should  leave 
observable  consequences  of  fine  structure,  which  for  more 
highly  disordered  material  will  become  progressively  weaker 
As  an  experimental  lest  of  this  effect,  high  quality  NbjSn  tun¬ 
nel  junctions  could  be  lightly  irradiated  thereby  increasing  ihe 
disorder  in  the  NbjSn,  and  the  observed  tunneling  DOS  moni¬ 
tored  for  alterations.  A  full  discussion  of  the  theory  and 
further  numerical  results  will  be  published  elsewhere.* 
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PARAMETER-FREE  CALCULATION  OF  THE  ENHANCED  SPIN 
SUSCEPTIBILITY  OF  Nb3Sn  INCLUDING  ELECTRON -PHONON  EFFECTS 

Warren  E.  Pickett 


Naval  Research  Laboratory 
Washington,  D.C.  2073S 

A  parameter-free  calculation  of  the  spin  susceptibility  x_p(T)  of  NbjSn, 
which  includes  static  and  dynamic  disorder  as  well  as  exchange-correlation 
enhancement,  is  presented.  It  is  found  that  xsp  is  only  157.  of  the 
measured  susceptibility,  and  its  small  temperature  dependence  cannot 
account  for  che  experimental  findings. 


It  is  a  well  established  fact  that  Che 
anomalous  variation^-  of  the  susceptibil¬ 
ity  x(T)  below  room  temperature  in  V- 
and  Nb-based  A15  compounds  is  directly 
correlated  with  the  superconducting 
transition  temperature  Tc:  d  log  x/dT 
is  largest  for  the  high  Tc  compounds. 
Attempts  to  explain  the  behavior  of  x(T) 
have  focussed  on  the  spin  susceptibility 
Xsp,  which  is  presumed  to  vary  due  to 
thermal  repopulation  of  a  narrow  peak  in 
Che  density  of  states  N(E)  near  the  Fermi 
level  Ep.  In  this  paper  I  present  Che 
results  of  a  parameter -free  calculation 
of  x9p  for  Nb^Sn  which  includes  (l)  the 
effects  of  both  static  and  dynamic  dis¬ 
order  and  (2)  electron-electron  inter¬ 
actions  through  an  ab  initio  Stoner  en¬ 
hancement.  It  is  found  that  X«p  is  only 
15X  of  the  measured  susceptibility, ^  and 
that  dXsp/dT  accounts  for  only  157.  of  the 
measured  temperature  variation. 

A  theory  for  the  electron  self -energy  Z 
has  recently  been  developed^  which  ac¬ 
counts  for  cases  where  N(E)  varies  on  the 
scale  of  a  typical  phonon  frequency,  such 
as  often  occurs  in  A15  compounds .  Both 
defects  and  phonons  are  included  in  a 
self-consistent  way.  Calculation  of  Z 
requires  solving  several  nonlinear  cou¬ 
pled  matrix  equations  involving  N(E)  and 
the  Green's  function  G.  The  input  to  the 
calculation  of  I  is  (1)  che  a2F  function 
derived  from  proximity  ettect  tunneling 
on  NbjSn  by  Wolf  ec  al.^  which  describes 
the  electron-phonon  interaction  (EPI), 

(2)  the  N(E)  function  calculated  for 
NbjSn  by  Klein  et  al.*  and  (3)  a  Lorentz- 
lan  defect  broadening  corresponding  to 
10  uDcm  of  residual  resistivity,  typical 
of  the  best  samples  of  Nb3Sn. 
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The  quasiparticle  density  of  states  is 
evaluated  from  the  relation 

nC^.T)--^"1!®  TrG((i)) ,  (1) 

and  therefore  includes  the  effects  of 
static  and  dynamic  disorder.  z  is  actu¬ 
ally  evaluated  on  the  imaginary  frequency 
axis,  leading  to  the  matrix  equations 
mentioned  above,  and  TrG  is  analytically 
continued  to  the  real  axis .  From  n  the 
unenhanced  spin  susceptibility  Xsp  is 
given  by 

*sJ(T)"uB2  /d“|liI±^if1(1J.T)-uB2n(c)  (2) 

where  f  denotes  the  Fermi  distribution. 

A  crucial  feature  of  the  present  calcu¬ 
lation  is  that  the  chemical  potential C 
is  retained  explicitly. 

Liu  et  al.^  have  shown  that,  with  local 
density  (LD)  theory,  exchange  and  corre¬ 
lation  enhancement  of  Xg2  may  be  included 
at  finite  T  through  the  Stoner-like  re¬ 
lation 

Xsp(T)"UB2lT(C)/a'r<c)^<;)}*  (3> 

where  I ( c)  is  the  thermal  average  around 
C  of  the  LD  Stoner  parameter  I,  which  can 
be  evaluated  from  the  band  calculations 
of  Klein  et  al.  Here  it  will  be  assumed 
that  Eq .  (3)  holds  when  n(c)  includes  the 
effects  of  disorder,  a  case  which  was  not 
considered  by  Liu  et  al. 

The  calculated  values  of  X  °  and  Xjp  are 
compared  in  Fig.  1  with  the  measured 
total  susceptibility  of  Reywald  et  al.^ 
Although  d  log  Xsp/dT*-lxlO“^K-^  is 
similar  to  the  measured  value  of 
d  log  X/dT,  Xjp  accounts  for  only  about 
15X  of  x.  This  last  result  is  surprising, 
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since  Klein  et  al.^  were  able  to  fit  X(T) 
with  only  two  parameters,  T  and  a  T-inde- 
pendent  ("orbital")  contribution  Y,  by 
ignoring  disorder  effects  and  including 
only  thermal  redistribution.  The  differ¬ 
ences  between  this  fit  and  the  present 
calculation  arise  in  two  ways.  First, 
the  calculated  value  (only  weakly  T-de- 
pendent)  of  T(C)"  2.3  mRy  is  only  about 
half  the  fitted  value  of  4.4  mRy.  Both 
the  Stoner  enhancement  (l-fn)’1  and  its 
temperature  variation  are  reduced  (rela¬ 
tive  to  the  fit  value)  by  more  than  a 
factor  of  two._  Secondly,  including  dis¬ 
order  reduces  n(c),  and  thereby  reduces 
X9°  as  well  as  further  reducing  the  en¬ 
hancement  factor. 

With  regard  to  this  second  point.  It  is 
notable  that  n(c)  is  reduced  more  by  the 
variation  of  C(T)  than  by  the  explicit 
defect  broadening.  The  effect  of  includ¬ 
ing  broadening  corresponding  to  10  uft  cm 
of  residual  resistivity  is  to  displace 
;(0)  upward  by  2  mRy  toward  a  region  of 
lower  N(E),  and  this  effect  alone  leads 
to  a  value  of  n(c)  which  is  85Y  of  N(Ep). 
Moreover,  merely  turning  on  the  EPI  fur¬ 
ther  shifts  electronic  soectral  weight 
our  of  the  peak,  shifting  C  upward  by 
another  3  mRy  (at  25K,  the  lowest  tempera¬ 
ture  considered  here)^  and  resulting  in  a 
further  decrease  of  77(;)  to  70Y  of  N(Ep). 
At  200K,  the  highest  temperature  consider¬ 
ed  thus  far^  I  find  c-10.5  mRy  (relative 
to  Eg)  and  n(C)-0.62  N(Ep)  .  The  temper¬ 
ature  variation  of  Xg^  shown  in  Fig.  1 
results  from  this  nearly  linear  decrease 
in  n(c)  with  temperature. 

The  results  presented  in  Fig.  1  indicate 
that  Xsp  cannot  account  for  experimental 
data  for  tft^Sn.  It  should  be  noted  that 
these  results  are  not  sensitive  to  de¬ 
tails  of  the  density  of  states  which  was 
used  (the  very  fine  structure  is  smoothed 
out  by  the  disorder)  as  is  the  case  when 
disorder  is  not  included,  and  only  the 
strength  of  a2F  (X«1.8)  is  important  in 
determining  C(T)  and  n(c).  It  now  ap¬ 
pears  that  the  orbital  susceptibility  may 
be  larger  than  Xsp  and  may  even  show  con¬ 
siderable  temperature  variation.  A  more 
detailed  discussion  of  all  aspects  of  this 
study  will  be  published  elsewhere. 

The  author  wishes  to  acknowledge  the  use 
of  unpublished  results  of  the  Ift^Sn  cal¬ 
culations  (Ref.  S)  as  well  as  the 


discussions  with  D.A.  Papaconstantopoulos 
and  C.M.  Soukoulis  about  disorder  broad¬ 
ening  and  with  B.M.  Klein  about  the  sus¬ 
ceptibility  fit  described  in  Ref.  7. 
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perature  as  discussed  in  the  text. 
Xexp  is  from  Ref.  2. 
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A  thermodynamic  formulation  for  the  electron  self-energy  is  given  which  is  applicable 
when  the  electronic  spectrum  possesses  structure  on  the  scale  of  phonon  frequencies,  pro¬ 
vided  only  that  the  ratio  of  phonon  phase  velocity  to  electron  Fermi  velocity  is  small. 
Electron-phonon,  Coulomb,  and  electron-defect  interactions  are  included  on  an  equal  foot¬ 
ing  and  it  is  shown  that  their  different  frequency  dependencies  lead  to  specific  effects 
on  the  Eliashberg  self-energy:  (a)  The  Coulomb  interaction  contributes  nothing  of  essence 
to  the  normal-state  self-energy  (in  this  isotropic  approximation)  but  retains  its  usual  de¬ 
pairing  effect  upon  the  superconducting  gap  function,  (b)  defects  affect  superconducting 
properties  primarily  through  a  broadening  of  the  electronic  spectrum,  and  (c)  phonons 
contribute  a  thermal  shift  and  broadening  as  well  as  the  mass  enhancement.  A  generali¬ 
zation  to  intensive  electron-phonon,  electron-electron,  and  electron-defect  interaction  con¬ 
stants  is  necessary  to  redevelop  an  intuition  into  the  effects  of  these  interactions.  The 
change  in  the  structure  of  the  Eliashberg  equation  due  to  a  nonconstant  density  of  states 
(DOS)  and  the  consequent  interplay  between  static  and  thermal  disorder  is  analyzed  in 
detail,  with  a  central  feature  being  the  change  in  frequency  dependence  of  the  self-energy 
compared  to  a  constant  DOS  solution.  The  effect  of  DOS  structure  on  the  superconduct¬ 
ing  transition  temperature  Tc,  which  is  manifested  in  the  defect  dependence  of  Tc,  is 
analyzed  in  detail.  Further  it  is  proposed  that  an  extension  of  the  self-consistent  Eliash¬ 
berg  approach  be  extended  above  Tc  to  determine  the  normal-state  self-energy  and  there¬ 
by  the  electronic  contribution  to  thermodynamic  quantities.  Phonon  broadening  is  shown 
to  affect  the  spin  susceptibility  at  finite  temperature.  Reinterpretation  of  several  of  the 
anomalous  properties  of  ,415  compounds  in  terms  of  the  present  theory  is  suggested. 

Several  aspects  of  the  theory  are  compared  to  experimental  data  for  NbjSn. 
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I.  INTRODUCTION 

Deeply  ingrained  in  the  formal  theory  of  the  in¬ 
teracting  electron-phonon  (e-ph)  system  in  metals 
are  two  simplifying  approximations.  The  first  is 
an  extension  of  the  adiabatic,  or  Bom-Oppen- 
heimer  approximation1  in  which  the  light  electrons 
are  considered  to  respond  instantaneously  to  the 
heavy  ions  (of  mass  A/).  Central  to  the  theory  of 
e-ph  systems  is  Migdal’s  theorem,2  which  demon¬ 
strates  that  nonadiabatic  effects  can  be  obtained 
accurately  by  low-order  Feynman- Dyson  perturba¬ 
tion  theory,  to  lowest  order  in  an  expansion  param¬ 
eter  of  the  order  of  (m/M)xn  « 1.  The  second 
simplification  is  the  assumption  of  a  constant  den¬ 
sity  of  states  (CDOS)  over  a  region  -Q  around  the 
Fermi  energy  EF ,  where  12  is  a  few  times  of  the 
mean  phonon  frequency.  This  approximation  al¬ 
lows  the  DOS  function  N(E)  to  be  approximated 
by  1 V(£f)  in  certain  energy  integrals.  The  two  ap¬ 
proximations  in  fact  are  related,  and  it  often  seems 
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to  be  assumed  that  Midgal’s  theorem  is  inapplica¬ 
ble  if  NIE)  is  not  constant  [to  within  (m/M)l/7] 
over  a  range  ±12  around  EF.  As  will  be  shown  in 
this  paper,  however,  there  exists  an  important  re¬ 
gime  within  which  the  CDOS  approximation  may 
be  relaxed  in  a  straightforward  manner  while  re¬ 
taining  Migdal’s  simplification.  The  resulting  gen¬ 
eralizations  of  the  CDOS  expressions  often  are  not 
intuitively  obvious,  and  the  consequences  involve  a 
reinterpretation  of  many  of  the  properties  of  this 
class  of  materials. 

That  structure  in  the  DOS  on  the  scale  of  12 
should  be  expected  in  crystals  containing  several 
transition-metal  atoms  per  unit  cell  can  be  deduced 
from  general  considerations.3  Elemental  transition 
metals  are  known  to  have  peak  structure  in  their 
DOS  which  may  be  only  a  few  tenths  of  an  eV 
wide.  A  compound  with  (for  example)  ten  atoms 
per  unit  cell  will  have  10  times  the  number  of 
bands  in  the  same  overall  bandwidth,  leading  to 
structure  on  the  order  of  hundredths  of  an  eV. 
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Recent  band-structure  calculations4  5  on  .-115  com¬ 
pounds  (eight  atoms  per  unit  cell)  indeed  have  veri¬ 
fied  structure  on  this  scale.  DOS  structure  can  be 
important  even  in  elemental  metals  if  the  Fermi 
level  EF  lies  near  a  van  Hove  singularity,  such  as 
occurs  in  Pd  and  Pt. 

The  most  intriguing  consequence  of  the  e-ph  in¬ 
teraction.  both  theoretically  and  technologically,  is 
superconductivity,  and  it  has  not  been  overlooked 
that  the  superconducting  transition  temperature  Te 
might  be  affected  by  DOS  fine  structure.  An  early 
study  of  A 15  compounds  was  carried  out  by  Bar- 
isic  and  co-workers0  in  the  weak-coupling  limit, 
where  it  was  noted  that  the  DOS  within  a  region 
-r(l  around  EF  contributes  in  the  determination  of 
Te.  Cohen.  Cody,  and  Vieland7  applied  a  strong- 
coupling  version  of  the  Koonce-Cohen8  formalism 
to  investigate  the  effect  of  model  A 15  DOS  func¬ 
tions  upon  Te. 

In  a  model  solution  of  the  Eliashberg  equation 
with  non-CDOS  effects  included  approximately 
Nettel  and  Thomas’  suggested  an  average  DOS 
given  by 


•Vnt  = 


f°'  del 
rt  eoZial) 


X  [XloiZl  +  .V(  —<oZ)]  /ln(  ilt/Te) 


(1.1) 


as  being  most  relevant  in  determining  Te.  Here 
Z(e>)  is  the  strong-coupling  renormalization  (real 
part),  nc»n  is  a  cutoff  frequency,  and  the  DOS 
function  will  be  taken  with  origin  at  EF  through¬ 
out  this  paper.  It  was  evident  to  Nettel  and  Tho¬ 
mas  that  there  is  great  utility  in  identifying,  when¬ 
ever  possible,  an  effective  value  of  XlEF)  for  use 
in  the  CDOS  version  of  Eliashberg  theory  rather 
than  having  to  deal  always  with  the  full  energy 
dependence  explicitly.  This  simplication  will  be 
explored  further  in  this  paper. 

Horsch  and  Rietschel10  obtained  numerical  solu¬ 
tions  to  the  Eliashberg  equation  in  the  small 
kTc/fi  limit  where  the  imaginary  part  of  the  self¬ 
energy  can  be  neglected.  They  found  an  enhance¬ 
ment  of  Te  relative  to  the  CDOS  value  T°  for 
placements  of  a  square  root  DOS  peak  near 
Ef  —  nZ(0).  This  was  interpreted  as  a  reduction 
of  the  repulsive  part  of  the  e-ph  interaction  due  to 
DOS  variation.  It  could  equally  well  be  considered 
as  an  enhancement  of  the  attractive  part  of  the  in¬ 
teraction  (which  occurs  at  w  <  fi>,  with  an  approxi¬ 
mate  value  of  the  enhancement  given  by  ;VNT. 

A  more  general  approach  was  taken  by  Lie  and 


Carbotte,"  who  calculated  a  functional  derivative 
&Te/&X(E)  for  several  superconductors.  Its  shape 
was  found  to  be  quite  insensitive  to  the  metal  con¬ 
sidered.  having  a  maximum  at  E  =  EF.  decreasing 
to  half  maximum  at 

\E-Ey,  =5-107; 

and  becoming  negative  (but  remaining  very  small) 
above 


£-£>  =50—707;  . 


This  (linear)  approach  leads  to  an  average  DOS 
*T'  dTt 


X 


LC 


=  /  dE 


5  JV<£) 


.V  (£) 


/ 


dXlEy) 


(1.2) 


which  is  useful  for  estimating  Tc  from  a  CDOS 
theory.  In  Eq.  (1.2)  dTe/dX(EF )  is  the  integral 
over  bTe/&N(E),  and  Tc  is  given  by 

r‘=r^c-'v'£^]^’  tlJ) 


Written  in  this  form  .VLC  clearly  has  an  interpreta¬ 
tion  as  an  effective  value  of  .V  (£f ).  Since  the  im¬ 
portant  contributions  to  the  integral  in  Eq.  (1.2) 
come  from  the  region  E  —  EF  <  10TC  ~ffe  (for 
strong  coupling  superconductorsl.Jor  nonpatholog- 
ical  DOS  functions  the  averages  .VNT  and  XLC 
should  be  similar. 

Whereas  an  average  over  the  DOS  may  be  suffi¬ 
cient  for  understanding  Tc,  for  many  properties  of 
interest  (e.g.,  the  T  dependence  of  the  spin  suscep¬ 
tibility  Y„p)  no  such  simplification  will  be  possible. 
In  this  paper  a  general  approach  for  obtaining  the 
electronic  self-energy,  and  thereby  the  thermo¬ 
dynamic  properties,  is  described.  Both  static  and 
dynamic  disorder  are  included  in  a  straightforward 
generalization12  of  the  usual  Eliashberg  approach, 
and  it  is  shown  that  each  type  of  disorder  contri¬ 
butes  to  the  shift  in  chemical  potential  as  well  as 
to  a  shift  and  renormalization  of  the  electronic 
spectrum.  Although  the  emphasis  here  will  be 
focused  on  the  superconducting  state,  implications 
for  normal-state  properties  also  will  be  discussed. 

A  preliminary  report1'  of  the  application  of  this 
theory  to  NbjSn  has  been  published  elsewhere. 

The  plan  of  the  paper  is  as  follows.  Section  II  is 
devoted  to  the  description  of  the  system  of  in¬ 
teracting  electrons,  phonons,  and  defects  which 
will  be  studied,  and  the  approximations  leading  to 
our  expression  for  the  electronic  self-energy  are 
discussed  and  justified.  The  treatment  of  the  ener¬ 
gy  dependence  in  the  Eliashberg  equation  is  given 
in  Sec.  III.  where  it  is  argued  that  the  band-energy 
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dependence  of  the  self-energy  should  be  negligible 
in  most  cases.  This  leads  in  turn  to  a  less  formid¬ 
able  numerical  procedure  for  solving  for  the  self¬ 
energy.  In  Sec.  IV  the  generalized  Eliashberg  sys¬ 
tem  of  equations  is  presented  and  the  alterations 
arising  from  DOS  structure  are  discussed.  As  a 
byproduct  of  retaining  the  band-energy  dependence 
of  the  e-ph  interaction  in  the  equations,  a  better 
understanding  of  the  Coulomb  pseudopotential  and 
of  the  role  of  impurity  scattering  is  obtained.  A 
discussion  of  several  aspects  of  the  theory  at  Te  is 
presented  in  Sec.  V.  In  Sec.  VI  the  spin  suscepti¬ 
bility  is  discussed  in  terms  of  a  self-energy  deter¬ 
mined  from  an  Eliashberg-type  equation  at  T>TC. 
Unlike  previous  uses  of  Eliashberg-type  equa¬ 
tions,14  which  of  necessity  require  an  infinite  sum¬ 
mation  of  diagrams  to  describe  appropriately  the 
appearance  of  a  gap  (superconducting  or  spin- 
density-wave)  in  the  spectrum,  this  novel  applica¬ 
tion  accomplishes  more  easily  what  could  be  ac¬ 
complished  (at  least  approximately)  in  some  finite 
order  of  perturbation  theory.  It  is  suggested  that 
this  approach  may  be  fruitful  for  studying  other 
thermodynamic  properties.  Indeed,  there  has  been 
little  concern  over  non-CDOS  corrections  above  Te 
with  the  exception  of  the  T  dependence  of  and 
the  elastic  constants,  both  of  which  show  anoma¬ 
lous  behavior  in  high-7"e  A 15  compounds. 


H.  GENERALIZATION  OF  THE 
ELIASHBERG  SELF-ENERGY 

The  usual  treatment  of  Eliashberg  theory  on  the 
imaginary  frequency  axis  is  clearly  set  forward  in 
the  literature.'*  We  will  provide  only  the  back¬ 
ground  necessary  to  clarify  the  generalizations  we 
propose  and  the  approximations  which  remain. 

The  system  we  consider  is  described  by  the  Hamil¬ 
tonian 

H  +ffph  +^f«.ph  +  f^imp  >  (2.1a) 


where 


k 

(2.1b) 

Hpji*  2nc*-Mc*  • 

(2.1c) 

X'Ft^>,r3'FJk. , 

(2.  Id) 

tf«-ph=  2  £*Av(bfc_ie\v  +  hfc-_i.v) 
k,k\v 

(2.10 


H, mp=  2  2  , 

where  CT  denotes  a  second-order  term.  In  the 
Nambu  scheme  the  spin  index  is  eliminated  in 
favor  of  the  two-component  field  operator 


Ckt 

t 

c-kl 


(2.2) 


in  terms  of  the  electron  annihilation  operator  cka. 
The  Pauli  matrices  are  denoted  by  t\,  r2,  and  r3. 
Electron-band  and  phonon-mode  indices  will  not 
be  exhibited  explicitly  except  where  necessary. 

In  this  form  of  the  Hamiltonian  'V  describes 
band  electrons,  for  which  the  electron-static  lattice 
and  electron-electron  interactions  have  been  includ¬ 
ed  in  a  mean-field  sense.  For  the  electron-lattice 
interaction  the  remaining  coupling  is  given,  to 
second  order  in  the  ion  displacement,  by  the 
electron-phonon  Hamiltonian  The  second- 

order  term,  which  has  not  been  displayed  explicit¬ 
ly,  is  required  to  keep  the  theory  translationally  in¬ 
variant.16  The  residual  Coulomb  interaction  be¬ 
tween  band  electrons  is  assumed  to  be  representable 
in  the  usual  four-body  form  and  is  denoted  by  V  in 
the  electron-electron  Hamiltonian  H,^.  The  effect 
of  this  term  on  superconducting  behavior  is  not 
understood  in  detail  but  fortunately  an  approxi¬ 
mate  treatment  in  terms  of  an  empirical  "Coulomb 
pseudopotentiai”  (see  below)  seems  sufficient  for 
most  purposes. 

The  “impurity”  term  ff,mp  represents17  n,  identi¬ 
cal,  randomly  distributed  imperfections  centered  at 
positions  Rj.  The  scattering  potential  K,mp 
represents  the  difference  between  the  potential  in 
the  region  of  the  imperfection  and  the  perfect- 
crystal  potential.  The  phonon  Hamiltonian  is 
expressed  in  terms  of  bare  phonons  (with  creation 
operator  bqj  of  frequency  which  are  dressed 
to  the  observed  frequency  jug,  by  the  band-electron 
polarization  as  described  for  example  by  Migdal. 

In  this  paper  we  will  not  consider  the  effect  of  lat¬ 
tice  imperfections  on  either  the  phonon  spectrum 
or  the  electron-phonon  coupling  g,  although  the  ef¬ 
fect  may  become  large  for  highly  disordered  ma¬ 
terials. 

The  electron  thermodynamic  Green’s  funotion  is 
given  by 
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—  (£*—  f)rj— Hk,iw9)  , 

(2.3) 

where  =(2 n  + 1) irT  and  the  self-energy  2  con¬ 

ventionally  is  given  by15 

J.{k,ita9)^i(aH[  \  —  Z(k,iam)\ 

+*(£,/«,  )T3+^(fc,i<ull)r,  .  (2.4) 

For  bulk  superconducting  properties  the  r2  term  in 
2  is  proportional  to  <t>  and  is  assumed  to  be  elim¬ 
inated  by  the  choice  of  phase.  Two  points  should 
be  noted  here.  First,  the  chemical  potential  f, 
which  is  determined  by  the  relation 

N,=*T'2tTrG(k,ioH)ei**‘  (2.5) 

k,n 

in  terms  of  the  number  of  electrons  N,  (e  is  a  posi- 


D  V  it,  i1 


ti  c  c 


PHONON  SCREENED  DEFECT 
COULOMB 

FIG.  1.  Proposed  phonon.  Coulomb,  and  defect  con¬ 
tributions  to  the  electron  self-energy  2.  The  double 
lines  indicate  the  renormalized  Green's  function  G<2). 
giving  a  self-consistent,  infinite  order  relation  for  2. 

tive  infinitesimal  quantity),  will  be  retained  expli¬ 
citly.  Second,  we  also  retain  the  energy  shift  X 
(not  to  be  confused  with  the  susceptibility)  and 
find  that  although  it  does  not  reduce  to  a  trivial 
shift  in  the  chemical  potential,  its  effects  are  usual 
ly  of  a  secondary  nature. 

The  Eliashberg  equation  for  the  self-energy, 
represented  diagramatically  in  Fig.  1  is15 


2(/c, /<!>„)=*  —  T  2  r3<7( k )r}  2  lS**vl20v(*  —k',i(o9—i<0 „•) 

*•.«•  v 


+  V(k,k’\itoH  — iav)+[iif  |  t{k,k',ia„)\  l/T)  5. 


The  first  term  results  from  coupling  to  phonon 
branch  v  with  renormalized  phonon  Green's  func¬ 
tion  Dv.  The  second  term  arises  from  the  dynami¬ 
cally  screened  and  Coulomb-vertex-corrected 
electron-electron  interaction  P,  which  is  discussed 
further  below.  The  third  term,  which  is  usually 
not  included,  results  from  neglecting  correlated 
multiple  scattering  between  distinct  static  defects, 
but  includes  multiple  scattering  from  a  single  de¬ 
fect  by  means  of  the  t  matrix17  f.  This  treatment 
strictly  applies  to  the  dilute  limit;  however,  the 
transport  theory  which  results  usually  is  found  to 
be  valid  to  rather  large  defect  concentrations. 

In  writing  the  phonon  contribution  to  2  as  in 
Eq.  (6),  Migdal's  theorem2  has  been  invoked.  In  its 
usual  form  this  theorem  states  that,  to  within  a 
quantity  of  order  ( m/M)xn ,  the  electron-phonon 
vertex  function  r,.Ph  can  be  taken  as  unity  (for  the 
purpose  of  calculating  the  electron  self-energy). 
Within  simple  metal  language  an  equivalent  expan¬ 
sion  parameter  is  U/E/.  This  had  led  to  specula¬ 
tion  that  Migdal’s  theorem  may  not  hold  for  sys¬ 
tems  which  have  DOS  structure,  and  in  some  sense 
an  effective  degeneracy  temperature  (Ef)  on  the 
scale  of  H.  As  Scalapinol5(cl  has  emphasized, 
however,  the  validity  of  the  approximation  rests  on 
a  small  value  of  the  ratio  of  phonon  phase  velocity 


cjq/Q  to  electron  group  velocity  vk.  Even  for  A  15 
compounds  with  large  N(Ef)  and  low  electron  ve¬ 
locities  the  ratio  of  these  quantities  averaged  over 
the  Brillouin  zone  is  (2— 3)X  10-2  (e.g.,u  { o>q/Q  ) 
—  5 X  105  cm/sec  and5  <ut>~<u*),/2si>j r 
=2x  1G7  cm/sec  for  Nb3Sn).  The  approximation 
rf.ph=s  I  can  break  down  in  the  immediate  vicinity 
of  a  van  Hove  singularity  vk  =0,  but  such_a  small 
number  of  electrons  within  j  Ek  —Ef  |  <  12  have 
<^k  <uq/Q  ^at  the  contribution  to  the  Brillouin- 
zone  sum  in  Eq.  (2.6)  will  be  negligible. 

Recently  Drozhov19  has  provided  more  insight 
into  this  problem.  The  situation  Drozhov  studied 
is  that  in  which  an  M  |  van  Hove  critical  point 
falls  at  Ef,  which  is  one  situation  for  which 
vk  <,wq/Q .  It  was  found  that  so  much  dispersion 
is  introduced  into  the  electronic  states  by  the  e-ph 
interaction  that  the  corresponding  quasiparticles 
are  for  the  most  part  removed  from  the  region  of 
the  critical  point  (and,  of  course,  are  highly 
damped),  The  quasiparticle  velocities  tend  to 
diverge  a:  the  critical  point.  Thus,  if  the  fully  re¬ 
normalized  Green's  function,  rather  than  its  bare 
counterpart,  is  used  in  Eq.  (2.6)  (as  is  done  in 
Eliashberg  theory  and  as  it  is  shown  in  later  sec¬ 
tions  must  also  be  done  in  the  normal  state),  it 
may  well  be  the  case  that  the  criterion  Vk  ><jq/Q 
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in  terms  of  renormalized  velocities  Vk  (if  they  can 
be  defined)  is  satisfied  even  near  critical  points 
vk  —0  where  Migdal’s  theorem  as  stated  is  unprov¬ 
en.  Assuming  this  to  be  the  case,  the  only  remain¬ 
ing  difficulty  is  with  optic  modes,  for  which  near 
zone-center  phase  velocities  (<jq  /Q  diverge.  It  ap¬ 
pears  that  a  small  enough  number  of  modes  are  in¬ 
volved  that  the  Brillouin-zone  sum  for  I  is  unaf¬ 
fected. 

Migdal2  also  demonstrated  (for  the  normal  state 
at  T  — 0 )  that  the  renormalized  Green’s  function  G 
appearing  on  the  right-hand  side  of  Eq.  (2.6)  can 
be  replaced,  in  the  CDOS  system  he  considered,  by 
the  unrenormalized  Green’s  function,  thus  allowing 
a  formal  solution  for  low  frequencies.  No  corre¬ 
sponding  simplification  is  possible  when  variation 
in  V  ( £)  is  not  negligible,  as  observed  by  Eliash- 
bergl5(hl  for  the  superconducting  state  in  which  the 
opening  of  the  gap  gives  a  sharply  varying  density 
of  quasiparticle  states.  A  central  point  of  this  pa¬ 
per  is  that,  when  the  energy  variation  of  ,V(£)  oc¬ 
curs  on  the  scale  of  phonon  energies,  2  must  be 
determined  self-consistently  from  Eqs.  (2.S)  and 
(2.6)  for  the  normal  state  as  well. 

The  three  contributions  to  2  in  Eq.  (2.6)  differ 
fundamentally  in  their  frequency  dependence.  The 
phonon  Green’s  function  D¥(iton)  has  its  important 
frequency  variation  for  |  »„  |  <  fi  and  decreases  as 
a »,  for  |  <u„  |  »fi.  The  Coulomb  interaction 
PUtoJ  varies  only  on  the  scale  of  the  electronic 
plasma  frequency  (peak  in  the  energy-loss  function) 
to pj~  10  eV.  The  impurity  contribution  is  energy 
conserving,  with  (1  /T)bn%H-  in  Eq.  (2.6)  represent¬ 
ing  the  6  function  6 (o>— &>')  in  the  thermodynamic 
formulation.  For  the  phonon  and  impurity  contri¬ 
butions  the  sum  over  |  a>„*|  in  Eq.  (2.6)  can  be 
truncated  at  a  cutoff  frequency  <ue  —  S  —  lOfi  with 
negligible  loss  of  accuracy.  The  practical  necessity 
of  using  a  minimal  frequency  range  for  solving  for 
2  has  prompted  the  folding  down  of  the  Coulomb 
potential  P  into  a  pseudopotential  U  as  described 
in  detail  elsewhere. 15,20  The  treatment  is  general¬ 


ized  somewhat  in  Appendix  A,  with  the  result 

l/U^l/P+.VtcjptnniWpt/V')  .  (2.7) 

This  relation,  and  the  effective  DOS  .V(o»pi),  will  be 
discussed  further  in  Sec.  IV  A.  Using  this  pseudo- 
potential  and  utilizing  the  evenness15"*1  in  frequen¬ 
cy  cj„  of  the  functions  Z,  X,  and  6,  Eq.  (2.6)  can 
be  reduced  to  an  .Vc  x.Vc  matrix  equation  with 
(2-Vc  —  1  )irT  =coc  (i.e.,  n  =0, 1,  ....  .Vf  —  1). 

III.  TREATMENT  OF  THE 
ENERGY  DEPENDENCE 

In  its  most  general  form  the  Eliashberg  equation 
is  extremely  difficult  to  solve  even  if  the  kernel  in 
large  parentheses  in  Eq.  (2.6)  is  known.  Besides 
the  four-dimensional  frequency-momentum  vari¬ 
ables  which  are  summed  over,  G  and  2  also  are 
matrices  in  band  index,21  although  this  fact  is 
nearly  always  ignored.  To  date  little  of  a  quantita¬ 
tive  nature  is  known  about  the  importance  of  off- 
diagonal  (in  band  index)  contributions  to  G;  how¬ 
ever,  cases  where  these  corrections  seem  to  be 
necessary  have  been  extremely  rare  (but  see 
Chakra  bo  rty  and  Allen21).  In  the  present  paper  all 
such  “band-mixing”  effects  will  be  neglected.  We 
will  concentrate  on  including  the  energy  depen¬ 
dence  of  the  band  density  of  states  within  an  iso¬ 
tropic,  band-diagonal  approximation. 

The  isotropic  average  A  (E)  of  a  wave- vector- 
dependent  quantity  A  (k)  is  defined  by 

A(E)  =  ^A(k)8(E-Ek)  ,  (3.1) 

t 

where  to  simplify  notation  we  have  introduced  a 
dimensionless,  normalized  averaging  function  given 
by 

8(E-e)=8(E-e)AV(E) . 

A  further  approximation  necessary  to  reduce  Eq. 

(2.6)  is 


- 1 

'2lG(k’,iuH  )B(k,k'-,ia)n  —ian)%lE  —  £*)5(£'  — Ek')=*G(E',ioH)B{E,E'\ian  ,  (3.2) 

kk' 

where22 

G (£,/«„ )_l£i<u,  —  (£—  f)rj— 2(£,/o>,) .  (3.3) 

B(k,k'\ian—i(on)  is  any  of  the  three  kernels  in  Eq.  (2.6)  and  £(£,£';/<u„  —ieo„  )  is  defined  by  Eq.  (3.2). 

The  self-energy  becomesl5u,,,5,el 

2(£,to»,)*r  2  f  dE,TiG(E’,ioH')TiN(E')(k(E,E’;ian—iu>n  )—li*(E,E')—\T‘(E,E')/wT\8m-\  , 

(3.4) 
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.  • 
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where 

k{E,E‘;io)n)- 


!  —  2  I  ?**'- 


lDv(k  -k'JwJbiE  -Ek)&{E'-Ek  )  , 


kk'y 


fi*(E,E')=2i/(k,k‘;0)6(E-Ek)S(E'-Ek  >  =  U(£,E’;0>  , 

**• 

f(£,£')  =  irn,  2  I  t(k,k", 0)  1 2S(£  -£*  )&£'-£*•)  . 

kk‘ 


(3.5) 

(3.6) 

(3.7) 


In  the  CDOS  limit  .V(£)-*V(£f ),  the  usual 
electron-phonon  coupling  constant  A.,  pseudopoten¬ 
tial  fim,  and  impurity  width  I"  are  given  by 
N  (EF)k(EF,EF‘, 0),  N{Ef)(A*(Ef,Ef),  and 
N(EF)f‘(EF,EF),  respectively.  It  will  become  ap¬ 
parent,  however,  that  the  normally  ubiquitous 
quantity  N(EF)  nowhere  appears  explicitly  in  this 
more  general  theory,  although  we  will  identify  in 
Sec.  V  an  effective  average  density  of  states  which 
is  useful  in  an  approximate  determination  of  Te. 
The  Coulomb  pseudopotential  /I*  is  discussed  fur¬ 
ther  in  Sec.  IV  A  and  in  Appendix  A. 

The  £  dependence  of  2  arises  solely  from  the  £ 
dependence  of  the  kernel 

kdf^E,E‘)shE,E‘)-^im^E,E^) 

—( l/nT)f  (£,£')  , 

which  in  turn  is  due  to  variation  with  energy  of 
the  character  of  the  wave  functions  and  scattering 
properties  but  not  the  density  of  states.  In  general 
we  expect  this  variation  to  be  small  compared  to 
DOS  effects,  although  this  question  deserves  fur¬ 
ther  study.  Information  on  the  £  dependence  of 
k(E,E)  can  be  gained  by  studying  the  £  depen¬ 
dence  of  the  mean-square  electron-ion  scattering 
matrix  element  /z(£)  which  enters  k.  Using  the 
calculated  electronic  structure  and  wave  functions 
of  Klein  et  a/.1  (to  be  used  in  the  calculations 
described  in  Sec.  V),  it  is  found  that  /*(£)  varies 
by  only  2%  in  a  4-mRy  region  centered  at  EF 
where  N(E)  changes  by  a  factor  of  3.  Neglecting 
this  small  energy  dependence  of  the  kernel  gives 


=  T  2  K^it»H-iaH)Fm 


(3.8) 


where  2.  etc,  and  the  density-of-states 

effects  are  confined  to  the  factor 

£.»  /  dE  rfj(E,io)n  )rtN  (£)  . 

Inverting  G-1  in  Eq.  (3)  leads  to 

-i<onZn -(£-£+*•„  )rj 


(3.9) 


TjG(£,J<i»,  )Tj  ’ 


C(E,iaH ) 


(3.10) 


with 


C(E,iwh  )  =  (£  —  £ 4-  )* -r  +  Ajj )  . 


(3.11) 


The  gap  function  is  given  by  A„  —d>„  /Z„.  Even 
and  odd  averages  over  .V(£)  result,  given  by 


=  —  f  dE 

if  * 


,V(£) 

C[E,ia>.) 


2„(o>;  +  A;  )1/* 


E-;+xn 


(3.12) 


with  each  average  incorporating  in  the  denomina¬ 
tor  C(E,it )  an  energy  shift  f  —  T,  and  broaden¬ 
ing  half- width  Z^itol  )1/z.  Both  functions  are 
even  in  <on.  The  following  sections  describe  the  ef¬ 
fects  of  a  nonconstant  .V  ( £)  that  are  included  in 
Nn  and  P„,  which  for  a  constant  density  of  states 
reduce  to 


IV. 


msf) 

o 


(3.13) 


IV.  GENERALIZED  ELIASHBERG  THEORY 

The  generalized  equations  for  the  Eliashberg 
self-energies  are 

•v* 

canZl,=cjl,+irT  2  ^'Nmem  +  r.\ne„(oH  , 


m  *0 


(4.1) 


X„  =  -rrT  2  a(£-2f?)Pmem-rPneHo>H 

m  *0 


(4.2) 


(2n  +  I)4.-  2  (A£  -2#Z’)-=- Nmemim/Zn 


m  m0 


0>n 


+  f  N.e.i.co./Z. 


(4.3) 


where  k1 1 '  and  the  DOS  “enhancement"  *?„  are  de¬ 
fined  by 


43 


1192 


warren  e.  pickett 


26 


•  (4.4)  R.ee(<y'i  =  Re{a>/[ti>2  —  A2(t<>)]l/2|  .  14.5’) 


A  nonzero  solution  to  Eq.  (4.3)  exists  only  for 
T  <TC.  However,  Eqs.  (4.1)  and  (4.2)  may  be  use¬ 
ful  at  higher  temperatures  as  well,  a  fact  which  is 
discussed  in  Sec.  VI. 

Converting  Eq.  (4.3)  into  an  equation  for  the  gap 
function  leads  to 

\  '  „ 

A„=trr  2  (A^;’-2/Z*hN 'mem 

m  *0 

v,  _ 

-5.m 

/-  0 

(4.6) 

Explicitly  K  is  given  by 
k(io)n  —i(om ) 

=2/dnna:Fm)/tn2+<u„-<o2]  i*7) 

in  terms  of  an  intensive  electron-phonon  spectra 
function  (i.e.,  coupling  per  electronic  state) 

a2F(E,E'-,Cl)=  2 

k.k\v 

xS  <£'-£*•> 


X5(fl— <Ufc_*',»)  ,  (4.8) 

and 

3zF(n)*S2F(£,£':ft)  |  £a£-=£(r  - 

Equations  (2.5),  (3.12),  and  (4.1)  — (4.3)  [or  (4.3) 
replaced  by  (4.6)]  form  a  system  of  coupled  non¬ 
linear  matrix  equations  which  must  be  solved  itera¬ 
tively  for  the  self-energies  at  each  temperature  of 
interest.  Their  solution  allows  (at  least  in  princi¬ 
ple)  the  direct  calculation21  of  electronic  thermo¬ 
dynamic  properties,1 51  fl  such  as  Te,  specific  heat, 
etc.  Numerical  solutions  of  these  equations  have 
been  presented  elsewhere13  (and  see  below)  but  a 
number  of  effects,  and  their  interpretation,  result¬ 
ing  from  a  complex  electronic  structure  will  be  dis¬ 
cussed  in  this  section.  The  care  which  must  be  ex¬ 
ercised  in  any  interpretation  of  the  imaginary  fre¬ 
quency  equations  is  exemplified  by  the  “enhance¬ 
ment"  defined  by  Eq.  (4.5).  Evidently  e„  lies  in 
the  range  0<  ;  e,  I  $  1.  The  designation  "enhance¬ 
ment"  is  appropriate  since,  when  continued  to  the 
real  axis,  the  real  part  of  this  function  gives  the 
well-known  DOS  enhancement  on  the  real  axis 


A.  The  Coulomb  pseudopotential 


From  Eq.  (4.2)  it  follows  that  the  contribution  to 
X  from  the  Coulomb  interaction  is  independent  of 
a>„.  This  results  in  a  constant  shift  Xc  in  all  quasi¬ 
particle  energies  and  consequently  a  compensating 
shift  £■ —  £-*-Xc  in  the  chemical  potential.  There¬ 
fore  the  Coulomb  contribution  to  X  in  Eq.  (4.2) 
can  be  disregarded  without  loss  of  generality.  On 
the  other  hand,  phonon  and  impurity  scattering 
give  rise  to  nonconstant  energy  shifts  X„  which 
must  be  retained  in  an  accurate  theory. 

The  pseudopotential  U  in  Eq.  (3.6)  is  given  as 
usual20  as  a  solution  to  an  integral  equation  de¬ 
scribing  the  folding-down  of  high-frequency 
scattering  processes  included  in  V.  The  zone- 
averaged  pseudopotential  jl*  satisfies  the  approxi¬ 
mate  relation  (see  Appendix  A) 

-L  =  ±+ivWln  ^  ,  (4.9) 

M  F  [ 

where  N(o>)  denotes  an  average  over  N(E)  on  the 
scale  of  to,  and  jl  is  given  in  terms  of  V  by  an 
equation  analogous  to  (3.6).  However,  according  to 
the  gap  equation  (4.6) _it  is  approximately  the  di¬ 
mensionless  quantity  V(<u4)jZ*  which  enters  into 
the  determination  of  the  gap  function  (and  hence 
rc).  Here  denotes  the  range  of  the  average  over 
.V(£)  appropriate  to  the  gap  equation,  which  will 
be  discussed  further  below.  In  the  limit  of  a  slow¬ 
ly  varying  DOS  near  Ep,  iV(<uA)—  .V(£f )  and 


J _ l_  ff(a)pl)  0>pt 

ji*  n  +  N(Ef)  °  oe 


(4.10) 


This  leads  to  a  dependence  of  n *  upon  .V(£f )  simi¬ 
lar  in  form  to  that  of  Bennemann  and  Garland,24 
who  used  p=0.26  and  .V(<up|)In(tOp|/&>c) 

=  3.85  (eV  atom)-1  for  transition  metals.  How¬ 
ever,  since  fi  in  Eq.  (4.10)  is  formally  proportional 
to  ,V(£f ),  which  can  vary  widely  within  a  given 
class  of  metals,  Eq.  (4.9)  is  the  appropriate  relation 
from  which  to  determine  Ji* .  The  quantities  /Z  and 
<V(w P|)  vary  slowly  within  a  class  of  metals,  with 
the  dominant  variation  arising  from  scaling  with 
bandwidth. 


B.  Impurity  scattering  and  “Anderson’s  theorem” 

In  the  CDOS  limit  .V(£f),  F„  —  0,  X„ 
vanishes,  and  the  chemical  potential  is  constant. 
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The  gap  equation  (4.6)  reduces  to  the  usual  form 
and,  for  T sN(EF )T— *0,  the  equation  (4.1)  for  Z 
does  likewise.  For  nonzero  T  the  static  impurities 
give  a  contribution  to  Z  related  to  disorder 
broadening  of  the  energy  spectrum.  Obviously  this 
has  no  effect  on  the  normal  state,  since  (in  the 
CDOS  limit)  there  is  no  structure  to  broaden.  The 
gm-induced  structure  in  the  superconducting  state, 
where  the  quasiparticle  DOS  is 

.V(£f)Re|tu/[a):-AJ(6))]-1/2j  , 

also  is  unaffected  by  the  broadening  since  the 
equations  for  Z  and  A  are  uncoupled  in  this  limit. 
This  is  a  manifestation  (indeed,  a  proof)  of 
Anderson's  theorem,23  which  states  that  supercon¬ 
ducting  properties  (of  CDOS  systems)  are  unaffect¬ 
ed  by  dilute  nonmagnetic  impurities. 

The  generalised  Eliashberg  equations  above  show 
that  this  statement  of  Anderson's  theorem  breaks 
down  when  structure  is  present  in  the  DOS,  due  to 
the  disorder  broadening  of  the  underlying  electron¬ 
ic  structure.  Although  T  does  not  enter  explicitly 
into  the  gap  equation  (4.6)  (due  to  the  energy- 
conserving  nature  of  this  interaction),  disorder 
scattering  can  lead  to  important  spectra]  broaden¬ 
ing  via  Z,  and  conceivably  to  non-negligible  spec¬ 
tral  shifts  X„  —  £.  The  effect  on  the  gap  and  on 
Tf,  is  transmitted  to  the  gap  equation  entirely 
through  the  resulting  set  f.VH  j. 


c.  r=rt 


At  T=TC  (A„-*0)  the  alteration  of  the  quasi- 
particle  density  of  states  due  to  the  gap  vanishes, 

w„/(<ui+Ai)1/2— sgrnu,  , 


and  the  expressions  (4.1),  (4.2),  and  (4.6)  simplify 
somewhat.  The  equation  for  the  (infinitesimal)  gap 
becomes  linear  and  decouples  from  the  remaining 
equations,  which  still  must  be  solved  iteratively  for 
.V„.  The  equation  determining  Tc  can  be  written 
[Eq.  (4.6)  with  e„— *  1] 


1 


-2A* 


-6. 


In  -r  1 

•V, 


Vc 

+  2  Zi 

ImO 


,ft 

•V, 


(4.11) 


where  Am  =  Am  /cjm .  Written  in  this  form  the  Te 
equation  is  a  straightforward  generalization  of  the 


CDOS  expression,  as  given  for  example  by  Allen 
and  Dynes.26  Ignoring  for  the  moment  corrections 
to  the  kernel  diagonal,  the  primary  change  is  the 
replacement  of  the  CDOS  eigenvector  A,  by  .V„  A„ . 
This  replacement  has  no  effect  upon  either  Te  or 
the  eigenvector,  whose  components  are  now  .V„  A„ . 
This  means  that  (at  this  level  of  approximation)  A„ 
will  be  altered  by  DOS  structure  in  proportion  to 
.V"1.  A  similar  correction  to  the  gap  occurs  at 
T  <  Tc,  and  this  indicates  that  the  analytic  con¬ 
tinuation  A(cu)  can  be  altered  substantially13  27  by 
energy  dependence  of  the  DOS. 

The  correction  due  to  DOS  structure  to  the  pho¬ 
non  contribution  to  Z„  —  1,  proportional  to  .V,  in 
Eq.  (4.1 1),  tends  to  be  canceled  by  Nn,  while  the 
factor  .V"1  multiplying  (2n  + 1 )  in  the  diagonal 
can  be  regarded  as  the  self-consistent  response  of 
A,  to  the  DOS  structure.  It  is  this  latter  effect 
which  gives  rise  to  much  of  the  correction  to  Te 
compared  to  the  CDOS  limit.  The  structure  of 
Eq.  (4.11)  will  be  clarified  further  in  Sec.  V  where 
an  approximate  Te  equation  is  discussed. 

V.  DISCUSSION  AND  CALCULATION AL 
RESULTS  AT  LOW  TEMPERATURE 

A.  How  much  DOS  variation? 

The  system  of  generalized  Eliashberg  equations 
has  been  applied  to  NbjSn  and  Nb  at  T~TC. 
Technical  aspects  of  the  numerical  solution  are  dis¬ 
cussed  in  Appendix  B.  For  Nb3Sn  the  corrections 
to  the  CDOS  results  are  substantial,  both  for  Tc 
and  for  the  self-energies  A,,  o>„Z„,  and  Xn.  Crys¬ 
talline  Nb  was  checked  as  a  possible  fringe  case, 
where  corrections  might  be  noticeable  if  not  really 
important.  The  DOS  of  Nb  is  characterized  by  a 
large  value  of  .WiEF)s  10  (Ryspin)-1  with  EF  ly¬ 
ing  20  mRy  above  a  peak,  and  the  energy  variation 
being  described  sufficiently  well  by  a  slope  of 
—  200  (Ry2spin)-1  over  a  range  of  ±20  mRy 
=  ±3\  I03  K  around  EF.  The  correction  to  Tc 
was  found  to  be  0.4%  (downward),  with  corre¬ 
spondingly  small  corrections  to  the  self-energies  at 
Tc.  Although  this  correction  is  well  below  the  ab¬ 
solute  accuracy  of  the  Eliashberg  equations  as  de¬ 
rived  (keeping  only  certain  classes  of  diagrams)  and 
applied  (especially  the  treatment  of  the  Coulomb 
interaction),  it  nevertheless  represents  faithfully  the 
accuracy  of  the  CDOS  approximation  in  Nb  and  in 
metals  with  similar  or  less  DOS  structure  near  EF . 
A  preliminary  report  of  the  NbiSn  calculations  has 
been  reported13  and  further  results  are  discussed 
below. 
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B.  Approximate  relations  for  Te 

In  the  CDOS  limit  Te  is  determined  by  azF  and 
/i*.  or  equivalently, 

rc  =  rc[^(^),M*]  ,  (5.1) 

where  arF  is  described  by  a  strength  A.  and  shape 
function 

g(o>)=(2/A*j)a2F(aj)  . 

Generalizing  the  pioneering  work  of  McMillan.28 
Allen  and  Dynes26  showed  that,  for  widely  varying 
strengths  and  shape  functions,  this  functional  of  g 
can  be  replaced  to  high  accuracy  by  a  function  in¬ 
volving  only  two  moments  <u)0(  and  on  of  g. 

Te  as  1  •  (5.2) 

One  of  the  most  useful  applications  of  such  an  ap¬ 
proximate  Te  equation  is  the  extraction  of  A.  from 
experimental  data  on  Tc  and  reasonable  estimates29 
of  a,,  o»;,  and  /z*. 

In  the  previous  sections  it  has  been  shown  that 
the  general  functional  form  of  Te  is 

Te  =  Tc[k.H*>),p;N(E),r]  .  (5.3) 

where  g  is  defined  analogously  to  g.  The  form  of 
the  gap  equation  (4.6),  together  with  the  work  of 
Allen  and  Dynes,  suggests  an  approximate  Te 
equation  of  the  form 

Tc  ~  Te(X,a}ot,co2,iim f ) ,  (5.4) 

where  iWK,  jVM,  and  iV4  are  averages  over  1V(£).  A 
few  general  features  of  these  averages  will  be  noted 
here. 

The  average  denoted  N is  that  which  should 
multiply  A.'  ~ 1  in  the  kernel  diagonal  in  the  approx¬ 
imation 

.V  .V 

(5.5) 

ImO  lm  0 

Since  X  decreases  as  a>P2  for  <om  »H,  iVx  will  cor¬ 
respond  approximately  to  a  Lorentzian  average  of 
iV(£)  over  a  frequency  half- width  H  around  EF 
(actually  £). 

In  the  off-diagonal  terms  the  corresponding  re¬ 
placement  is 

2  .  (5.6) 

*•0  mm  0 

However,  since  Am  decreases  rapidly  with  increas¬ 
ing  um,  even  changing  sign  for  com  <H,  <V4  aver¬ 
ages  over  a  smaller  region  than  does  ,V\,  for  exam- 


•v« 

2 


pie,  on  the  order  of  fl/2. 

The  average  .VM  multiplying  /I*  is  less  straight¬ 
forward,  since  /I*  has  no  frequency  cutoff  to  make 
a  partial  sum  coverage.  Numerical  solutions20- 31 
show  that  A„  approaches  a  negative  constant  value 
A,  for  large  n  (for  nonvanishing  pm).  Obviously 
N„  also  approaches  a  constant  average  DOS  value 
-V  „ .  If  both  i V,  and  A,  have  approached  their 
asymptotic  value  at  n  =k  <SC,  it  follows  that  ,VM, 
given  by 


2'^A,  + jiV.AJn^/fc) 

o _ 

*-!  _ 

2  A„  +  TA.lndVj/k ) 

o 


(5.7) 


is  weakly  dependent  on  the  cutoff  Se.  This  cutoff 
dependence  is  artificial  and  arises  from  the  incon¬ 
sistent  treatment  given  to  p*  in  folding  jl  down  to 
p*.  The  energy  integral  leading  to  Eq.  (4.9)  does 
not  take  into  account  the  full  variation  of  iY(£) 
whereas  the  “unfolding”  to  the  cutoff  Ne  in  Eq. 
(5.7)  is  taking  this  energy  dependence  into  account. 
It  is  straightforward  to  correct  this  discrepancy  by 
generalizing  Eq.  (4.9),  that  is,  by  using  the  relation 
(A5).  However,  the  following  approximation  will 
suffice  to  illustrate  some  effects  of  DOS  structure. 

Owing  to  the  rapid  decrease  of  A„  with  n  and 
the  (2/t  + 1)-1  factor  in  A„,  the  DOS  average. 


M  _  /  M  _ 

2  ^»  /  2  * 

n -0  •  n.O 


(5.7') 


for  small  Af  “converges”  (i.e.,  becomes  stationary 
with  respect  to  At)  rapidly  before  the  cutoff  depen¬ 
dence  mentioned  above  becomes  a  consideration 
and  leads  to  the  unwanted  limit  JVM— » Nm.  In  a 
preliminary  report13  of  this  work  applied  to  NbjSn 
(also  see  below)  the  first  term  alone 


iYM(Af=0)  =  N{  irTeZ0 )  s  lVefr 

was  used  in  the  CDOS  theory,  giving  an  approxi¬ 
mate  value  T which  could  be  compared  to  nu¬ 
merical  solutions  of  the  full  system  of  matrix  equa¬ 
tions  for  Te.  In  essence  T^a  is  given  by  (for  f =0) 

T?=  T^iN'ifX.u^^N'M 2* )  (5.8) 

(although  numerical  solutions  to  the  CDOS  equa¬ 
tions  were  actually  used  rather  than  the  Alien- 
Dynes  equation).  Note  that 

itTcZq  a  0. 85 — 0. 90<y ;  a  fl  /2 


■  ^  ^ 
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FIG.  2.  Superconducting  transition  temperature  Tc  of 
NbjSn  vs  Fermi-level  position,  calculated  within  dif¬ 
ferent  approximations.  Full  line,  solution  from  full 
energy-dependent  equations;  dashed  line,  solution  assum¬ 
ing  CDOS  iV(£)=iV(£f);  dot-dash  line,  solution  assum¬ 
ing  CDOS  E)—NtS(  (see  text). 

for  NbjSn  <7V  =  18  K,  Z0=2.7,  <u2=175  K).  The 
results  of  this  approximation  are  shown  in  Fig.  2, 
where  it  is  evident  that  Tf1  reflects  the  trends  of 
the  numerical  solution  quite  well  as  the  assumed 
value  of  EF  is  varied,  although  being  — 1.5  K  too 
large.  The  crude  approximation  7?  resulting  from 
the  assumption  N(E)=N(EF)  is  seen  in  Fig.  2  to 
give  a  much  stronger  N(EF)  dependence  than  actu¬ 
ally  occurs. 

For  nonzero  impurity  scattering  each  DOS  aver¬ 
age  N„  has  its  broadening  half-width  u>nZn  in¬ 
creased  by  TNI  [see  Eq.  (4.1)].  A  low-order  ap¬ 
proximation  can  include  this  by  increasing  by 
each  of  the  widths  which  determine 
Nx,  N&,  and  NM,  of  which  the  former  must  be  cal¬ 
culated  self-consistently  (at  least  in  principle).  The 
calculations  described  below  result  from  numerical 
solutions  to  the  full  energy-dependent  equations, 
however. 

The  quantities  JV4,  and  ,VM  offer  a  reason¬ 
able  possibility  of  obtaining  a  realistic  value  of  Te 
without  solving  nonlinear  coupled  matrix  equa¬ 
tions.  Unfortunately  it  is  a  monumental  task  to 
ascertain  even  an  appropriate  form  of  the  approxi¬ 
mate  Te  equation  envisioned  in  Eq.  (5.4),  as  this 
would  involve  full  solutions  for  a  wide  variety  of 
shapes  of  N(E)  and  £(o>)  as  well  as  coupling 
strengths  X,  ji *,  and  T.  Worse,  it  begins  to  appear 
hopeless  to  extract  unambiguous  information  from 
an  experimental  value  of  Te  in  the  face  of  so  many 
unknown  parameters.  The  only  compensation  is 
that  the  behavior  of  Te  with  the  level  of  imperfec¬ 


tion  [7*e(f )]  can  be  incorporated  into  the  analysis. 
A  more  realistic  approach  for  interpreting  experi¬ 
mental  data  is  to  use  a  single  average  A^ff.  with 
N'n=Wi.  being  a  good  choice,  in  the  CDOS  ex¬ 
pression  of  Allen  and  Dynes.  This  in  effect  is  the 
usual  practice,  except  that  the  value  of  A.tfr=  ,Veff X 
extracted  in  this  manner  may  differ  from  the  mass 
enhancement  value  derived  from  heat  capacity  or 
other  measurements. 

Although  the  DOS  average  appropriate  for  an 
approximate  Te  equation  given  above  differs  in  de¬ 
tail  from  that  of  Nettel  and  Thomas,9  the  spirit  of 
the  approximation  is  similar  to  theirs.  The  more 
rigorous,  but  linear,  approach  of  Lie  and  Car- 
botte11  for  calculating  8TC /SNiE)  follows  from  the 
present  formulation  by  calculating  the  correction 
A  7V  arising  from  an  infinitesimal  variation  of  en¬ 
ergy  E  from  a  constant  DOS  given  by 

AJV(E')=e6(E'-E)  . 

It  could  be  of  interest  to  study  the  effects  of  the 
supralinear  terms  for  a  particular  DOS  function 
using  the  present  theory.  However,  the  Lie- 
Carbotte  results  provide  an  important  intuitive 
grasp  of  the  general  behavior  of  3Te/8N(E)  as 
well  as  a  first  approximation  for  numerical  calcula¬ 
tions. 

C.  Defect  dependence  of  Tc  and  N(EF) 

It  was  observed  in  Sec.  V  B  that  experimental 
data  on  7V[  f  ]  [i.e.,  Tc  versus  defect  concentration 
h,,  since  f =irn,  1 12 1  from  Eq.  (3.7)]  provide 


-80  -SO  -«0  -Z0  0  20  «0  60  80 


E-Er  (mfty) 

FIG.  3.  Comparison  of  electron  and  phonon  energy 
scales  for  NbjSn.  (a)  a:£(<u)  from  tunneling  data  by 
Wolf  «  al.  (Ref.  32)  and  .V(£)  function  calculated  by 
Klein  et  al.  (Ref.  5),  shown  on  the  same  mRy  scale,  (b) 
Same  iV(£)  function  on  a  scale  of  —  1  eV. 
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direct  information  on  the  energy  variation  of 
.V (£).  Calculation  of  Tf(  F)  has  been  carried  out 
for  Nb3Sn  using  the  DOS  function  calculated  by 
Klein  et  al.s  For  alF  we  have  assumed 

a2F{ai)=N(ioph)a1F{(j)  ,  (5.9) 

<wph  =  300  K,  with  a:Fi w)  derived  from  tunneling 
data  by  Wolf  et  al.12  Both  .V(£)  and  a2F  are 
shown  in  Fig.  3  on  the  same  scale  to  allow  an  easy 
comparison  of  their  energy  variations.  For  jl *  the 
value  p*  =0. 16  determined  by  Wolf  et  al.  has  been 
used  in  the  form  X  ( Ef  )jZ =0. 16. 

To  compare  with  experimental  data  the  defect¬ 
broadening  half-width  F  was  related  to  the  residual 
resistivity  po  by 

Po=4irr/il2.  (5.10) 

(l;={4ire1/3)N(EF)v'F  ,  (5.11) 

with33  Mlp=4.0  eV.  The  experimental  data  on 
p(Te)=p{ 0)  has  been  convened  to  the  ‘‘true"  resi¬ 
dual  resitivity  po  by  the  use  of  the  parallel  resistor 
formula 


I 

piO) 


- r -  . 

Po  Ptnu 


(5.12) 


pmmi  =  1 50  pH  cm,  since  Weismann  et  al.1*  have 
found  that  po  given  by  this  expression,  and  not 
p< 0),  is  approximately  linear  with  damage-inducing 
radiation  dose.  The  theoretical  prediction  is  com¬ 
pared  to  four  sets35  of  experimental  data  in  Fig.  4. 
The  calculation  is  consistent  with  experimental 
data  for  residual  resistivities  below  75  pH  cm.  For 
larger  disorder  the  theoretical  prediction  ap- 
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FIG.  4.  Calculated  Tt  of  NbjSn  (full  curve)  vj  "true" 
residual  resistivity  po  (i.e.,  from  parallel  resistor  formula) 
compared  with  experimental  data  (Ref.  35).  The  theory 
is  not  expected  to  be  accurate  for  large  defect  concentra¬ 
tion  (po>75  pflcm)  since  changes  in  5  Fare  not  includ¬ 
ed  in  the  theory. 


proaches  the  asymptotic  limit  7V  =9  K  whereas  the 
experimental  values  saturate  at  Tc  ~3  K.  This  in¬ 
dicates  that,  for  po>  75  pflcm,  the  disorder  is  de¬ 
creasing  the  electron-phonon  coupling  52(o>) 
and/or  hardening  the  phonon  spectrum  F( at),  nei¬ 
ther  of  which  is  included  in  the  present  theory 
which  treats  a2F  as  independent  of  disorder. 

The  correspondence  between  defect  concentra¬ 
tion  and  residual  resistivity  (T— po)  is  not  as  sim¬ 
ple  as  Eq.  (5.10)  implies.  It  should  be  observed, 
for  example,  that  T  in  Eq.  (5.10)  must  be  deter¬ 
mined  self-consistently  from  the  relation 

r  =  iV(Df,  (5.13) 

where,  as  will  be  shown  below,  the  F  dependence 
of  X(D  is  strong.  Therefore,  although  f  increases 
linearly  with  disorder,  F  is  distinctly  sublinear. 
However,  the  factor  fl2  in  a  suitably  generalized 
version  of  Eq.  (5.10)  also  contains  the  factor  N(D 
so  the  DOS  dependence  of  po  cancels.  As  a  result 
the  defect  concentration  dependence  of  p0  goes  in¬ 
versely  with  that  of  vF,  a  result  analogous  to  an 
equivalent  observation  for  the  phonon-limited  resis¬ 
tivity  made  previously  by  Allen.36  This  result  can 
be  written 

po»3f/e2u/[r(f)J  .  (5.14) 

We  find  by  calculation,  however,  that  uj?(D<x  1/ 
iV( D  [since  Q2(D  is  constant  to  within  2%], 
which  makes  Eq.  (5.14)  numerically  equivalent  to 

p(r-[4jr/n2(£/E-)],V(r)f  .  (5.14') 

Conceptually  Eqs.  (5.14)  and  (5.14')  are  not  at  all 
equivalent  and  the  former  relation  gives  the  correct 
physical  picture. 

In  Fig.  5  the  effective  defect  broadening  (or 


r„  (mRy) 


FIG.  5.  Self-consistently  determined  defect  broaden¬ 
ing  width  T  (equivalently  residual  resistivity  p0) — full 
line — vs  disorder  r<>,  determined  from  .V(£)  shown  in 
Fig.  3.  The  dashed  line  gives  the  CDOS  approximation 

r«r0. 
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FIG.  6.  Effective  density  of  states  .V(  D  vs  residual 
resistivity  p0  for  Nb,Sn.  Solid  line:  theory,  including  e- 
ph  interaction  with  As  1.8.  Dashed  line:  theory,  no  e- 
ph  interaction.  Symbols:  experimentally  inferred  values 
(Ref.  38),  assumed  two  different  constraints  on  S(EF)vF 
with  varying  disorder. 

equivalently,  p0)  is  plotted  versus 

r0=JV(£f)f  =  jrn1JV(£f.)|r2|  . 

Owing  to  broadening  of  the  peak  in  N(E)  near  EF, 
r  increases  more  slowly  than  that  for  a  CDOS 
,V(£)=jV(£f),  shown  by  the  dashed  line.  Howev¬ 
er,  for  po>50  pil  cm,  T  is  approximately  linear 
with  defect  concentration,  just  as  the  data  of 
Weismann  et  al ,34  for  po  approximately  linear 

with  radiation  fluence  in  this  range. 

The  behavior  of  the  factor  of  N(EF)  in  the 
linear  specific-heat  coefficient  y  versus  residual 
resistivity  may  provide  more  direct  information  on 
N(E)  than  does  Tc(p'a),  since  only  a  straightfor¬ 
ward  Lorentzian  broadening37  of  AT(£)  is  involved. 
In  the  few  cases  which  have  been  studied  systemat¬ 
ically,  however,  y  has  been  inferred  instead  from 
the  temperature  dependence  of  the  upper  critical 
magnetic  field,3*  which  involves  independent  infor¬ 
mation  (or  assumptions)  about  the  values  of  ma¬ 
terial  constants  such  as  k,  vF,  mean  free  path  /, 
etc.,  as  well  as  an  assumption  about  the  behavior 
of  band-structure-related  quantities  with  defect 
concentration.  The  resulting  values  of  N(  D  for 
M>jSn  obtained  by  Ghosh  et  al.,3*  resulting  from 
the  assumption  that  N[EF)vF,  or  secondly,  vF  it¬ 
self,  remains  constant  with  increasing  disorder,  are 
compared  in  Fig.  6  with  that  calculated  from  the 
DOS  function  in  Fig.  3.  Calculated  values  of 
;V(D  are  shown  both  for  the  (usually  assumed) 
case  where  only  defect  broadening  is  taken  into  ac¬ 
count  (dashed  line),  and  the  more  general  case  in 
which  e-ph  broadening  is  included  (full  line). 

The  calculated  values  including  e-ph  broadening 
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FIG.  7.  Effect  of  defect  broadening  on  the  .V(£) 
function  shown  in  Fig.  3.  Vertical  lines  indicate  the 
respective  Fermi  levels,  determined  from  conservation  of 
electrons.  Circles  denote  .V(  D  =«  .V  ( £f ).  The  best  sam¬ 
ples  of  NbjSn  have  residual  resistivities  p0=  10  pH  cm. 
which  precludes  extreme  energy  variation  of  the  density 
of  states. 

are  in  slightly  better  agreement  with  the  empirical 
values  obtained  assuming  \(EF)vF=- const  [howev¬ 
er,  the  argument  given  above  in  the  discussion  of 
7V<po>  indicates  N(EF)vF)= const  would  be  the 
proper  assumption].  The  analysis  of  Ghosh  et  al. 
of  course  also  needs  generalizing  along  the  lines  of 
the  present  theory.  It  is  clear  from  Fig.  6  that  ex¬ 
trapolation  from  the  existing  data  to  obtain  a  “per¬ 
fect  crystal  value”  of  A’(£f )  may  not  be  valid, 
since  a  kink  may  occur  at  or  below  10  pH  cm. 

The  degree  of  broadening  of  the  full  .V  (£)  curve  is 
illustrated  more  clearly  in  Fig.  7.  A  peak  of  width 
3—4  mRy  or  less  is  virtually  lost  for  a  residual 
resistivity  of  10  p£l  cm,  which  corresponds  to  the 
highest  quality  samples  of  NbjSn.  Such  a  large 
perfect-crystal  value  of  NiEF)  implies  a  tendency 
toward  defect  formation,  which  serves  to  lower 
NiEF )  and  thereby  decreases  the  number  of  elec¬ 
trons  at  high  energy,  resulting  in  a  more  stable  ma¬ 
terial. 


D.  Defect  dependence  of  Tc  : 

Previous  studies 

There  have  been  several  previous  applications  of 
broadening  to  account  for  properties,  especially  Tc, 
of  A 1 5  compounds.  Most  of  these  have  been 
phenomenological,  as  typified  particularly  by  the 
work  of  Mattheiss  and  Testardi.39  40  The  studies 
of  Aleksandrov,  Elesin,  and  Kazeko,41  and  Huang, 
Chu,  and  Ting,41  however,  have  used  a  more  fun¬ 
damental  approach.  Both  groups  studied  Tc  using 
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a  weak -coupling  formalism  and  model  DOS  func¬ 
tions  to  illustrate  the  effect  of  defect  broadening  of 
V(£)  upon  Tt. 

Aleksandrov  et  al.*[  purported  to  explain  the  sa¬ 
turation  in  NbjSn  of  Te  at  3  K  for  high  defect 
concentration  solely  in  terms  of  DOS  broadening. 
The  present  generalized  theory  and  realistic  DOS 
function  shows  saturation  at  —  8  — 9  K  if  only  de¬ 
fect  broadening  is  taken  into  account.  In  addition, 
this  limit  depends  on  the  DOS  function  on  only  a 
large  (  —  1  eV)  scale,  for  which  all  band-structure 
calculations*4 5  give  results  similar  ( modulo 
5  —  10  "fc  bandwidth-type  differences  arisings  from 
different  exchange-correlation  approximations,  etc.) 
to  those  of  Klein  et  al. 5  which  were  used  here.  A 
similar  limit  ( ~  7  K)  follows  from  the  DOS 
broadening  study  of  Soukoulis  and  Papaconstanto- 
poulos,*2  if  Tc(  10  p.1 1cm)  is  normalized  to  the  ex¬ 
perimental  value  of  18.6  K.  To  understand  the  3- 
K  saturation  it  is  necessary  to  invoke  a  weakening 
in  the  strength  of  a2F  by  approximately  a  factor  of 
2.  This  weakening  may  result  from  weaker  cou¬ 
pling  or  harder  phonons.  In  fact,  both  of  these  ef¬ 
fects  have  recently  been  observed  in  tunneling  stud¬ 
ies  of  NbjAl  (Ref.  43)  and  Nb3Ge  (Ref.  44)  with 
varying  degrees  of  disorder. 

Huang  et  al.*1  assumed  a  singular  one-dimen¬ 
sional  DOS  and  obtained  results  similar  to  those  of 
Aleksandrov  et  al.  and  the  present  results.  They 
concluded  that,  due  to  DOS  variation  near  £> ,  a 
high  Te  compound  can  have  a  small  value  of 
.V(£f ).  A  look  at  Fig.  7  shows  that  this  cannot  be 
the  case  for  A  15  compounds  with  appreciable 
amounts  of  disorder,  as  in  all  samples  of  Nb3Ge, 
Nb3Al,  or  for  that  matter,  Nb3Sn.  For  p0=  10 
fid  cm,  i.e.,  a  clean  Nb3Sn  sample,  the  broadening 
half-width  T  is  1.6  mRy =250  K,  which  is  roughly 
equal  to  the  maximum  phonon  frequency  12.  Thus 
,V(  12 ),  which  determines  Te  (see  Sec.  V  B),  cannot 
bemuch  different  from  the  value  of  .V(£)  at  EF 
[,V(D  in  the  present  notation].  A  very  low  p0, 
high-7^  crystal  may  have  a  relatively  low  value  of 
S(EF )  if  a  large  peak  lies  within  -12  of  EF,  how¬ 
ever. 

Other  mechanisms  involving  defect  broadening 
have  been  proposed  to  account  for  the  degradation 
of  Te  by  defects.  Meisel  and  Cote*5  suggested  that 
an  assumed  inability  of  phonons  with  wavelength 
longer  than  the  electron  mean  fre*  nath  to  scatter 
electrons  effectively  (and  hence  bind  Cooper  pairs) 
could  account  for  the  defect  dependence  of  high  Te 
compounds.  This  cannot  explain  the  increase  in  Tc 
with  po  in  low-  Te  materials,*4  however.  Moreover, 
if  this  “phonon-ineffectiveness”  concept  is  accepted 


it  is  very  difficult  to  understand  the  rather  large 
Tc  and  measured  a~F  of  highly  disordered  Mo 
(Kimhi  and  Geballe,  Ref.  46).  Ruvalds  and 
Soukoulis*7  have  attributed  most  of  the  high  Tc  'at 
least  in  Nb3Ge  and  V3Si)  to  an  acoustic  plasmon 
mechanism  which  decreases  in  strength  as  the  elec¬ 
tronic  spectrum  is  broadened.  However,  there  is 
strong  evidence  from  tunneling  measurements  that 
the  electron-phonon  interaction  itself  is  strong 
enough  to  account  for  the  high  Tc  in  V3Ga  (Ref. 
48),  Nb3Sn  (Ref.  32),  NbjAl  (Ref.  43),  and  NbjGe 
(Ref.  44),  so  there  is  now  little  reason  to  except 
the  acoustic  plasmon  mechanism  to  apply  to  any 
A  15  compounds. 

VI.  NORMAL-STATE  SPIN  SUSCEPTIBILITY 

In  Sec.  IV  it  was  shown  that  in  general  Tc  does 
not  provide  a  measure  of  N  ( EF ),  but  rather  of  an 
average  of  the  DOS  over  a  region  nTcZ0  or  larger. 
In  this  section  we  investigate  what  information 
about  ,V(£)  is  contained  in  the  spin  susceptibility 
Y,p.  For  the  sake  of  generality  the  full  k  depen¬ 
dence  of  the  expressions  will  be  retained  as  far  as 
possible. 

A.  General  considerations 

The  magnetization  M  is  given  by  the  difference 
in  number  of  spin-up  and  spin-do  .vn  electrons 
times  the  moment  per  electron, 

M(S,T)=nB(Nt'X-N'A) 

-  (6.1) 

a 

with 

(6.2) 

k,n 

and 

GZ'(k,io)n)~  ia)H-(Ek-$-p.BoH) 

-2 „(.k,io)„) ,  (6.3) 

where  H  is  the  magnetic  field.  Luttinger*4  has 
shown  how  this  "self-evident”  exact  expression  for 
.Vf  can  be  derived  diagrammatically. 

The  measured  zero-field  susceptibility  is  given 
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»-2nlT'2GHk,i(oH)[^  +  ^k,ian)]  , 

k,H 

where  we  have  used  Eq.  (6.3)  and  the  expansion 

Za[k,iQ)„)—Z{ktici)H)—ng<7H'S.'(k,ial,)+0{.H1)  . 


(6.5 


(6.6 


Equation  (6.5)  is  an  exact  expression  for  in  terms  of  the  renormalized  Green's  function  and  the  field 
derivative  of  the  self-energy.  As  can  be  seen  from  Wolffs50  diagrammatic  theory  for  Xip,  the  G:  term  in 
Eq.  (6.5)  arises  from  the  simple,  but  renormalized,  bubble  diagram  in  Fig.  8(a),  whereas  the  G:T  term  arises 
from  all  vertex  corrections  (i.e.,  diagrams  with  other  than  self-energy  insertions  to  G).  With  I  given  by  the 
diagrams  of  Fig.  1,  we  can  write,  introducing  J„  =  -dGp  1  /dlfigcrH)  |  H  _0=./_a  —J, 


J(k,io*u)s* 


1- 


d 

d  (flgaH) 


T  2  kfgik.k’Ua), ,  - iton  )Ga(k',ia>„ •) 

k'n’ 


(6.7 


=  1-f-T  ’2tXdi(k,k';icjn—itan‘)G1(k',i(i)n)J(k',i(i)m) 

k*a' 


difigoH ) 


k^k.k'Uoin-ia^) 


G{k‘,iu>n) 


ff-o 


i6.7r; 


The  second  term  in  Eq.  (6.7*)  describes  “ladder- 
type”  vertex  corrections  (see  below),  with  a  typical 
graph  shown  in  Fig.  8(b).  The  third  term  repre¬ 
sents  field  corrections  to  the  effective  electron- 
electron  interaction,  with  the  lowest-order  phonon 
contribution  shown  in  Fig.  8(c).  There  is  of  course 


FIG.  8.  Susceptibility  diagrams  discussed  in  the  text. 
For  simplicity  the  renormalized  Green’s  function  is 
denoted  by  single  full  lines  rather  than  double  lines  as  in 
Fig.  1.  Other  notation  is  as  in  Fig.  I. 


a  large  class  of  diagrams  for  Y,p  which  are  not  in¬ 
cluded  in  Eq.  (6.7)  by  virtue  of  the  approximation 
for  2  from  which  this  equation  is  derived. 

For  T— *0  Luttinger49  has  shown  that,  as  long  as 
yia)  [2(«)sm(w)— iyia)]  vanishes  as  fast  as  a r 
as  <y— *0,  M  is  given  by 

M=fiB'2)<7Gi;-fkt,)  (6.8 

ka 

where  6  is  the  unit  step  function  and  € ka  is  the 
renormalized  energy  given  by 

•  <6-9 


The  susceptibility  is  given  by 


Y>p(r=0)  =  2^a2 

k  o 


d*ka 

dH 


5 i;-6k9)  . 


(6.101 


The  derivatives 

^k&ko~  ^  k^ko  ~  B  kmv(k,d  ko  —  b  ) 

—dum(k,o>)  |  UmAk0-{V k  6 ko  • 


(6.11 
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dH  ~  dH  3  H 


m,{k,6  ha-*' 


d6k. 


at  the  Fermi  surface  6  k(r  =  £  become 


7*^4<t  =  [  v*  -3tm„!fe.0)]/(  1  -A.* 

^*<7  dm„{k,0)  1  / 

d//  ~~Ms  3<Mstf)  j/vl 


l/ll-X.*)  , 


,  /  3m,(k,0)  i  / 

Wr=o».to.v,f.(!,,^rj/ 

Here  u*  =  v4  '  and  vk=vkvk,  and  the  angle 
brackets  (  )  denote  a  Fermi*surface  average. 

Equation  (6.10)  expresses  Vsp  in  terms  of  a  renor¬ 
malized  density  of  states,  including  mass  enhance¬ 
ment  from  Eq.  (6.9),  and  a  susceptibility  per  state 
—dSka/dH  which  also  includes  a  mass  enhance¬ 
ment  factor  [Eq.  (6.14)].  Equation  (6.15)  expresses 
Y,p  in  terms  of  the  bare  factor  ,V(£)  in  which  the 
total  mass  enhancement  cancels  out.  That  the 
phonon  mass  enhancement  cancels  out  was  first 
recognized  by  Quinn  and  Ferrell,51  who  interpreted 
the  cancellation  as  the  result  of  each  of  the  spin- 
split  bands  carrying  its  own  mass  renormalization 
with  it  rigidly. 

Reverting  to  the  isotropic  approximation  for  G. 
the  relation  for  the  thermal  distribution  function 

/(£-f)=r2C(£,(ai,)e'V  (6.16) 

It 

can  be  differentiated  to  yield 


where  \k  =  —3 um  1  k,cj) 1  is  the  usual  mass  re¬ 
normalization.  Equation  '6.10)  becomes 


V  km  ( k, 0)/vk 


stein  phonon  spectrum  of  strength  A.,  where  it  was 
found  that,  with  increasing  A.,  —df/'dE  becomes 
ever  more  long  ranged,  resulting  in  considerably  in 
creased  averaging  over  .V(£)  as  well  as  a  greater 
variation  of  £(  T)  given  by  the  isotropic  form  of 
Eq.  (6.2): 

.V,  =2  f  dE  f(E  -fhV(£)  .  (6.2 

For  even  moderate  values  of  A.,  of  the  order  of  0.5, 
/differs  significantly  from  /0  and  results  in  an 
enhanced  temperature  dependence  of  Y,p  in  Eq. 
(6.18).  Bhatt53  has  found  previously  that  a  low- 
temperature  expression  leads  to  an  enhancement 
!  1  +  A.)  of  the  7*  dependence  of  T,p. 

Returning  to  the  expression  (6.7)  for  the  en¬ 
hancement  /,  neglect  of  the  field  dependence  of  the 
interaction  leads  to 

Jii(o„)=  1 +  T  2  -‘a*-) 


dflE-C ) 


=  —T 2  GHE,i<an ) , 


the  convergence  factor  being  unnecessary  here. 
Equation  (6.5)  for  T,p  can  therefore  be  written 


f  dE  -■ 


rl£-L> 

dE 


X.Y(£)(  1  +  £'(£,£)]  .  (6.18) 

with  S'(£,D  being  defined  from  Eq.  (6.5)  by  this 
equation. 

The  crucial  feature  of  this  expression  is  that  in¬ 
teractions  can  alter  /(£)  drastically  from  its  free- 
particle  counterpart 

/„(£!  =  [  exp(£/D+ l]*1  . 

The  behavior  of  /  has  been  presented  elsewhere53 
for  the  case  of  electrons  interacting  with  an  Ein- 


x/d£  N  (E)G2[Eri(i)„')J 


in  the  isotropic  approximation.  As  was  the  case 
for  the  self-energy  in  Sec.  III.  the  (assumed)  in¬ 
dependence  of  A.cl-f  on  Ek  and  £*■  leads  to  an  £- 
independent  J.  Also  as  was  found  for  2,  the  dis¬ 
tinct  frequency  dependences  of  the  Coulomb,  pho¬ 
non,  and  defect  interactions  lead  to  different 
characteristic  behaviors  for  J. 


B.  Coulomb  interaction 

In  the  presence  of  only  Coulomb  interactions  'in 
the  approximation  of  Secs.  II  and  III)  f  is  essen¬ 
tially  unchanged  from  /0.  In  the  range  of  interest 
in  Eq.  (6.18).  of  the  order  of  rf.  the  frequency 
dependence  of  the  screened  exchange  interaction. 
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i 


n 


GENERALIZATION  OF  THE  THEORY  OF  THE  ELECTRON- 


which  will  be  denoted  as  usual  by  /,  can  be 
neglected.  Equation  (6.19)  becomes  (J,  =Ji  =/) 


/=!+/ fdE  - 


f(£-£) 


\N(E\J 


«l/[l 


where 


(6.20),  Eq.  (6.18)  becomes  the  usual  Stoner-en¬ 
hanced  expression  for 

C.  Interaction  with  defects 

Defects  drastically  alter  the  occupation  of  crys¬ 
talline  eigenstates,  even  at  T  =0.  Assuming  a  con¬ 
stant  density  of  states,  it  is  straightforward  to 
show  that  /  is  given  (in  the  dilute  limit)  by  the 
convolution  of  /o  with  a  Lorentzian  of  width 

r=N(EF)T : 


*«*<£)=  fdE  - 


df(E -C) 


LV(£)  (6.21) 


/(£)=  fdE'UE') 


(T/ir) 

(£-£')2+r2 


is  the  effective  density  of  states  at  £  at  temperature 
T.  With  the  enhancement  /  =  l  +  2'  given  in  Eq. 


If  .V(£)  varies  on  the  scale  of  N(EF) f ,  T  must  be 
determined  self-consistently  as  discussed  in  Sec. 
VC. 


Defects  give  rise  to  a  susceptibility  enhancement  given  by 

l+(f  /ir)f  dEMB)Gi(EJw.irU»m)-l/[l-r,U*M)]  , 

where  It  is  given,  after  partial  integration,  by 


7j(/o).)=-(f/ir) / dE^j~G(E,ian) . 


Clearly  no  enhancement  occurs  for  a  constant  NiE).  To  provide  an  estimate  of  the  magnitude  of  this 
correction  we  evaluate  the  lowest  order  correction  SY<)11,  pictured  in  Fig.  Rid), 

8Ti"=2Mi(f/i rlfj  ^GUcJo,,  )2G(k‘,ia>n)2 

h  kk' 

=2^(f  /ir)  2  )G(k',iaH ) 

kk-  dEk  dEk-  n 

a^lif/i r)g  Ek  l_Ek-T2lG{k,io>n)-G(k',ia>H))eim't 

.  2,*,  d  d  /<£*)-/(£*•) 


The  k  dependence  of  the  self-energy  has  been 
neglected  in  writing  the  energy  denominator  as 
£*■—£*•  Using  the  isotropic  approximation  and 
performing  partial  integration  in  each  variable 
gives 


5T|n*Va  —  /  dEN’(E) 


X  /  dE’N'(E')‘ 


r(£)— /(£') 
£’-£ 


We  consider  first  the  case  where  iV'(£) 
»dN(E)/dE  can  be  taken  to  be  constant  in  the 
important  range  around  EF.  Setting 


N‘(E)=AN(Ef )/W  ,  (6.2' 

where  A  is  a  constant  of  order  unity  and  W  is  the 
bandwidth,  leads  to  the  result 

8*:"=2 ^V(£fl42^£.  (62| 

Thus  defects  may  lead  to  a  nonvanishing  but  usu¬ 
ally  small  contribution  to  the  Stoner  7.  For  com¬ 
pounds  with  several  transition-metal  atoms  per 
cell,  however,  W  may  represent  a  subband  width 
such  that  T/IFis  not  small,  in  which  case  it  is 
necessary  to  do  the  integral  in  Eq.  (6.26)  more 
carefully.  Nevertheless,  we  expect  this  contribu¬ 
tion  to  X  generally  to  be  secondary  to  the  change 
in  .V(£f )  in  the  simple  bubble  term  due  to  defect 
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broadening.  Finally,  it  should  be  noted  that  this 
defect  contribution  to  Y,p  need  not  always  be  posi¬ 
tive  as  in  Eq.  >6.23).  If  EF  occurs  at  a  peak,  or  a 
dip,  in  .V(£).  in  which  case  .V'(£)  and  .V’(£')  in 
Eq.  (6.26)  have  opposite  signs  in  the  important  re¬ 
gion  of  integration,  the  correction  to  ¥sp  may  be 
negative. 

D.  Interaction  with  phonons 

It  is  not  generally  recognized  that  phonons  can 
seriously  alter  the  thermal  occupation  of  crystalline 
eigenstates  if  the  e-ph  interaction  is  even  moderate¬ 
ly  strong.  It  has  been  shown  elsewhere54  that  e-ph 
interactions  shift  the  spectral  weight  of  low-energy 
electronic  excitations  drastically,  and  as  a  result 
the  thermal  occupation  /( £)  broadens,  correspond¬ 
ing  to  occupation  of  higher-energy  bare  electrons 
and  holes.  The  resulting  modification  of  the  first 
term  in  Eq.  (6.18)  can  be  substantial.  Since  this  ef¬ 
fect  of  the  e-ph  interaction  on  /(E),  as  well  as 
many  of  its  implications,  is  described  elsewhere,52 
only  phonon  effects  on  the  susceptibility  will  be 


addressed  here. 

Previous  treatments53'56  of  the  effect  of  e-ph  in¬ 
teractions  on  the  susceptibility  have  approached 
the  problem  from  the  viewpoint  of  a  phonon  con¬ 
tribution  to  the  Stoner  interaction  parameter  I. 
These  treatments  have  all  assumed  a  CDOS  sys¬ 
tem,  however,  and  in  the  general  case  it  is  not  clear 
that  this  “Stoner  r  viewpoint  provides  a  useful 
approach.  To  see  this,  consider  the  Eq.  (6.7') 
neglecting  the  field  dependence  of  the  interaction: 
The  frequency  dependence  of  X  from  the  e-ph  in¬ 
teraction  is  essential  to  the  correct  evaluation  of 
this  term  and  results  in  a  frequency-dependent 
enhancement  J(iion ).  As  a  result  J  cannot  be  writ¬ 
ten  in  the  simple  (constant)  form  (6.20)  as  is  the 
case  for  Coulomb  interactions,  nor  even  in  the  still 
simple  (but  frequency-dependent)  form  (6.23)  as  for 
defects. 

We  content  ourselves  here  with  studying  briefly 
the  first  correction  to  ¥,p  from  e-ph  interaction 
(beyond  self-energy  insertions  to  G).  This  correc¬ 
tion  is  shown  in  Fig.  8(e)  and  is  given  by 


f*~  ^  G'(k,i<gm  )A Ak  —  k'\icjn  —  i(on  )G‘‘(k\i(i)n-)  . 

kn 

k’n‘ 

In  the  isotropic  approximation  this  becomes 


(6.29) 


=  >  f  dE  f  dE’ X’(E)G(E,ia)JN'(E')G(E',ia>'  )  . 


(6.30) 


If  :V(£)  varies  on  the  scale  of  fl,  no  further  sim¬ 
plification  of  this  expression  is  possible.  In  partic¬ 
ular,  is  not  bounded  in  magnitude  by  any 
small  parameter  and  it  may  be  of  either  sign. 

Since  this  contribution  to  ¥,p  is  the  first  [Fig.  8(e)] 
in  the  e-ph  ladder  series  leading  to  a  generalized 
Stoner  enhancement,  the  e-ph  contribution  to  a 
“Stoner  /“  may  be  large  and  of  either  sign. 

If  iV'(£)  can  be  approximated  by  a  constant 
N(EF)/W  over  the  range  EF  —  W  to  EF  W,  and 
is  negligible  otherwise,  Eq.  (6.30)  can  be  evaluated 
for  an  Einstein  phonon  spectrum  a1F(co) 
=(Xfl£/2)5(<u  —  SlE).  The  energy  integral  gives 

/ dE  S'(E)G  (E,iton ) 

— -isr[.V(£/1)/W']sgna>(t  .  (6.31) 


For  W  the  frequency  sums  lead  to  the  result 


h\s(ef) 


A .SlE  w 

TiF  '"n7 


(6.32) 


suggesting  that  a  ladder  summation  of  such  contri¬ 
butions  would  lead  to  a  contribution  to  the  Stoner 
I  given  by 

kClE  w 

/'  ph==  4  w  lnn7  •  (6-33) 

Oscillation  in  .V'(£)  in  transition-metal  compounds 
is  likely  to  severely  reduce  the  dependence  of 
on  the  cutoff  W,  so  the  net  result  is  likely  to  be 
dominated  by  DOS  structure  near  EF. 

It  is  notable  that  Eq.  (6.33)  is  of  exactly  the 
same  form  deduced  by  Fay  and  Appel55  from  the 
same  diagram  without  invoking  the  isotropic  ap¬ 
proximation  used  here.  Fay  and  Appel  also  em¬ 
phasized  that  other  diagrams  may  give  phonon 
contributions  to  ¥,p  of  the  same  order  as  that  re¬ 
sulting  from  Eq.  (6.33).  Although  the  contribution 
(6.30)  to  V,p  is  not  bounded  by  any  small  quantity 
like  Cl/W,  there  is  still  no  assurance  that  other  di- 
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a  grams  will  not  give  important  contributions. 
However,  Eq.  (6.30)  has  the  advantage  that  the 
electron  Green’s  function  is  dressed  by  phonons, 
implicitly  accounting  for  an  infinite  subset  of  dia¬ 
grams  for  X^. 

Kim,5*  taking  another  point  of  view,  has  con¬ 
sidered  the  magnetic  field  dependence  of  the  pho¬ 
non  frequencies  arising  from  changes  in  electronic 
screening  of  the  ion-ion  interaction.  The  change 
A (Oq  depends  on  the  spin  polarization,  which  may 
be  Stoner-enhanced  by  the  Coulomb  interaction. 
Kim  finds  the  result  that  the  phonon  contribution 
to  7  may  itself  be  enhanced  by  the  same  factor, 
thereby  greatly  increasing  the  importance  of  1,.^ 
for  nearly  magnetic  metals.  Figure  8(c)  gives  the 
lowest-order  diagram  which  arises  from  field 
dependence  of  the  phonon  spectrum.  Recently 
MacDonald  and  Taylor57  have  suggested  that  re¬ 
finements  of  Kim’s  theory  will  not  lead  to  correc¬ 
tions  as  large  as  were  envisioned  by  Kim. 


E.  Combined  effects  of  interactions 

A  realistic  calculation  of  Xv  requires  that  all 
three  of  the  interactions  discussed  above  be  taken 
into  account  simultaneously.  The  ladder  diagrams, 
e.g..  Fig.  7(b),  can  only  be  fully  included  by  a  nu¬ 
merical  solution  of  Eq.  (6.19),  and  corrections 
beyond  the  ladder  approximation  present  additional 
difficulties.  It  might  be  expected  that  a  reasonable 
first  approximation  would  be  to  retain  only  the 
Coulomb  part  of  X^t  in  Eq.  (6.19)  but  include  de¬ 
fect  and  phonon  contributions  to  G.  The  set  of  di¬ 
agrams  included  in  this  approach  is  typified  by  the 
Coulomb  ladder  diagram  of  Fig.  8(f),  and  the  cor¬ 
responding  contribution  to  Xv  is  given  by 


(6.34) 


with  iVrff  given  by  Eq.  (6.21). 

The  results  of  the  application  of  this  expression, 
with  7  calculated  from  local-density-functional 
theory,  have  been  described  elsewhere.5*  It  was 
found  that,  although  the  temperature  dependence 
was  qualitatively  similar  to  the  data  of  Reywald 
et  al.,i9  both  the  calculated  absolute  magnitude 
and  T  dependence  were  only  — 15%  of  the  experi¬ 
mental  values.  The  orbital  susceptibility  can  cer¬ 
tainly  account  for  much  of  the  discrepancy  in 
magnitude,  as  well  as  for  some  of  the  T  depen¬ 
dence  which  arises  from  the  variation  of  £  with 
temperature.  However,  it  seems  likely  that  disor¬ 
der  contributions  to  ¥,p  other  than  those  entering 
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Eq.  (6.30)  through  will  be  necessary  for  a 
more  quantitative  theory. 

VII.  CONCLUSION 

A  generalized  formulation  of  the  Eliashberg  ap¬ 
proach  to  the  electron  self-energy  has  been  devel¬ 
oped  which  is  valid  for  crystals  for  which  the  elec¬ 
tronic  spectrum  varies  on  the  scale  of  phonon  fre¬ 
quencies.  The  formulation  has  been  kept  on  the 
imaginary  frequency  axis  where  it  is  numerically 
tractable,  and  further  numerical  application  of  the 
theory  to  ^15  compounds  will  be  presented  else¬ 
where.  Although  there  is  no  difficulty  in  extend¬ 
ing  the  formalism  to  the  real  axis,12,60  the  resulting 
calculations!  difficulties  make  this  approach  unat¬ 
tractive  when  detailed  results  are  wanted. 

The  present  generalization  of  Eliashberg  theory 
proposes  a  system  of  equations  for  describing  the 
superconducting  onset  at  Te  as  well  as  the  gap  and 
renormalization  functions  below  Tc.  Calculations 
at  and  below  Te  have  shown  that  tunneling  spec¬ 
tra13,27  as  well  as  thermodynamic  behavior40  can  be 
significantly  altered  by  DOS  structure.  In  addi¬ 
tion,  it  is  proposed  that  the  self-consistent  ap¬ 
proach  of  Eliashberg  be  extended  above  Tc  to 
determine  the  normal-state  electronic  self-energy. 
This  theory  lays  the  foundation  for  a  more  unified 
picture  of  the  relationship  between  high- 
temperature  superconductivity  and  the  anomalous 
normal-state  behavior  of  A\5  compounds,  as  in  the 
conclusion  that  the  temperature  dependence  of  the 
susceptibility  is  increased  substantially  by  the 
strong  electron-phonon  interaction. 

The  satisfying  agreement  between  the  calcula¬ 
tions]  results  discussed  in  Sec.  V  and  the  experi¬ 
mental  data  indicates  that  (1)  the  present  proposal 
for  the  self-energy  (Fig.  1 )  is  adequate  at  low  tem¬ 
perature,  and  (2)  the  DOS  for  NbjSn  calculated  by 
Klein  et  al.  is  realistic.  The  disagreement  between 
the  theoretical  and  experimental  susceptibility  is 
certainly  too  large  to  be  ascribed  to  an  incorrect 
theoretical  DOS  function.  A  conceivable  source  of 
this  discrepancy  is  that  this  ansatz  for  the  self¬ 
energy  is  inadequate  for  temperatures  approaching 
the  Debye  frequency;  it  is  known  (see  Allen,  Ref. 
21)  that  even  in  CDOS  systems  there  are  further 
adiabatic  corrections  to  I  which  cannot  be  shown 
to  be  small.  However,  I  consider  it  more  likely 
that  the  expression  (6.34)  for  Xv  is  insufficient. 

An  important  area  which  has  not  been  addressed 
in  this  paper  is  the  normal-state  specific  heat  C„, 
which  has  been  used  often  to  extract  empirical 
values  of  NiEp).  Let  us  assume,  as  the  simplest 
situation,  that  the  specific  heat  can  be  divided  into 
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electron  and  phonon  contributions,  as  was  shown 
to  be  the  case  for  CDOS  systems  by  Grimvall.01 
Generalizing  Grim vall's  general  relation  Eq.  '  A7>  * 
along  the  lines  of  the  present  theory  leads  to  addi¬ 
tional  contributions  to  Cv  t  due  to  DOS  variation, 
arising  from  the  temperature  variation  of  C  and  of 
yl<u).  Grimvall01  has  already  concluded  that  Cv  t 
is  “far  from  linear  at  T ^sTe"  for  certain  strong 
coupling  superconductors.  In  addition,  many  .-i  IS 
compounds  are  known  to  have  highly  anharmonic 
phonons,  which  complicates  the  identification  of 
C0.ph.  Thus  much  theoretical  work  remains  before 


the  measured  specific  heat  can  be  used  to  deduce 
detailed  information  about  variation  in  ,V(£). 
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APPENDIX  A 

The  usual  expression15  20  for  the  Coulomb  pseudopotential  U  generalizes  in  the  present  case  to 
-  -  _  VUc,k"\ia>n)U(k",k'\0) 

(/(k,k';0)=Rfc,k';0)-i-r  T  T - ~  - - - 5—  .  (A1 

The  treatment  to  this  point  ignores  the  possibility  of  spin  fluctuation  contributions  and  also  assumes15101  that 
V  is  essentially  positive,  eliminating  the  consideration  of  enhancement  of  superconductivity  by  plasmon  or 
exciton  mechanisms.  In  principle,  the  dynamically  screened  Coulomb  interaction  V  includes  Coulomb  renor¬ 
malization  and  Coulomb  vertex  corrections. 

For  the  “frequencies"  on  and  energies  £*■•  of  interest  in  Eq.  (Al)  it  suffices  to  ignore  the  self-energy 
corrections  to  the  denominator.  Performing  the  energy  surface  averages  leads  to 

jS  *(£,£')=*/!(  £,£')+£  2  /  dE"}HE,E")  ^[E"\  .  1A21 

M.  >m.  E 


To  make  contact  with  later  notation  £*(£,£’) 
st7(£,£';0)  and  /!(£,£' Is  P(E,E',Q)  have  been 
introduced,  and  also  the  frequency  dependence  of 
V  has  been  ignored. 

Assuming  further  that 


fi,  I  £  j ,  |  £'  |  <o»p| 

/Z(  £,£')=«  0  otherwise  U 

and  that  N(E)  vanishes  for  I  £  |  >  wpl,  one  finds 
l//x*=»l//I-l-L(<ue,<Up)) ,  ( A 


where 


L{(ve,(0pi)sT 


T  f"p'  jr  -V(£) 


Hae,apt)ssN((opi  )H<oe,a)vi) 

*  -Vlwpi  )ln(<Wp ,/<uc ) ,  (A t 

where  iV(o>pj)  is  the  mean  DOS  over  a  region 
Ef±a p)  and  Ha,a')  (introduced  by  Allen13'*1)  is 
the  function  L  with  N(E )  replaced  by  unity. 
In  the  (strict)  CDOS  limit  ^(0^)— iV(Ef ),  the 
usual  expression 

l/M*  =  l/M4-ln(cup|/<tfe)  (A  7 

is  recovered,  where  n  =  !V(Ef)ji  and  n * 
3sN(EF)(i9.  However,  the  CDOS  limit  of  Eliash- 
berg  theory  applies  if  .V(£)  is  approximately  con¬ 
stant  over  the  range  EFzue>  in  which  case  a  more 
precise  expression  for  n*  is  given  by 


l /fi •  -  1  /n  +  (iV(Wp,)AV(£f  )]ln(a)p,/«c )  . 


For  a  specific  model  of  iV(£),  L  can  be  evaluated 
explicitly.  Given  the  level  of  approximation  lead¬ 
ing  to  Eq.  (A4)  however,  it  suffices  to  note  that 


This  relation  is  discussed  further  in  the  text. 
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APPENDIX  B 

The  calculations  described  in  this  paper  require 
the  self-consistent  solution  of  the  Ne  X.Vf  matrix 
equations  (2.5),  (3.12),  and  (4.1)— (4.3).  Two  tech¬ 
niques  used  in  the  present  computations  will  be 
described  briefly  here. 

These  equations  have  the  feature  that  the  low- 
frequency  (small  |  o)n  | )  rows  and  columns  are  in 
some  sense  more  critical  for  obtaining  convergence 
than  those  of  larger  frequency.  It  has  been  found 
that  the  technique  used  by  Allen62  of  using  the 
solution  of  the  m  x  m  subsystem  as  an  approxi¬ 
mate  solution  of  a  larger  subsystem  provides  an  ef¬ 
ficient  method  of  providing  an  iterative  solution  of 
the  system  of  equations.  The  sequence  m  — 1,2, 

4,8, . . . ,  Ne  used  by  Allen  has  been  used  here, 
with  Ne  =  128  for  all  calculations  described  in  this 
paper. 

The  only  arbitrary  assignments  in  this  procedure 
occur  in  setting  initial  values  of  N„  and  P„.  The 
iteration  is  begun  by  setting  N0=X{EF),  P0= 0, 
then  calculating  Z0,  To,  Aq,  and  £°\T),  the  first 
approximation  to  £(  D  (see  below).  Then  N0  tnd 
Pq  are  recalculated  until  self-consistency  is  ob¬ 
tained.  In  preceding  from  the  solutions  of  m  X  rn 

- 1 


subsystem  to  the  2m  x  2m  subsystem,  the  assign¬ 
ment 


m <j <2m 


(Bl) 


was  made.  A  search  for  more  sophisticated  extra¬ 
polations  of  these  functions  did  not  result  in  signi¬ 
ficantly  more  efficient  solutions. 

A  crucial  step  in  attaining  reasonably  accurate 
solutions  to  this  system  is  finding  a  procedure  for 
determining  £(  T)  from  Eq.  (2.5).  The  free  Green’s 
function  G0,  for  instance,  gives  rise  to  the  Fermi- 
Dirac  distribution  function  /„ 


/o<£)=  T  2  G0lE,ieo„  )e‘a't  (B2) 

H 

only  after  an  infinite  summation;  truncation  of  the 
frequency  sum  in  Eq.  (B2)  at  any  point  leads  to  a 
function  which  is  nonexponential.  In  the  present 
calculation  £  has  been  evaluated  as  follows. 

For  the  m  xm  subsystem,  m  =  1,2,4, .... Ne, 
we  make  the  definition 


Zj  =Aj  =Xj  =0,  j>m  .  (B3> 

Then  Eq.  (2.5)  can  be  written 


JV.*2r2  / <f£fV<E)ReG(£,/a^;£)*'*,,, 

«*o 

-27-Re  / dEN(E)  2  jG0(£,/®.;f)+[(?(£,ia>;£)-Go(£,i«;£)]| 

■  20 

-2/d£/0(£-£W(£)+2TRe  2  /  dE NiE)[G(E,ion\S)-G0(E,io>n,$)]  •  (B4> 

■  -o 

The  relation  G(£,— /«„)— G  (£,/«,)•  has  been  used  in  writing  the  sums  over  at„>0  only.  The  dependence 
on  £,  the  interacting  chemical  potential  has  been  displayed  explicitly  to  emphasize  that  introducing 
Go(£,i6>„;£)  and  /0(£  —  £)  is  merely  a  mathematical  method  for  evaluating  the  infinite  summation.  In  par¬ 
ticular, 

2  /  d£/0<£  -£)AT(£)*Ar,  .  (BS) 

Finally,  the  noninteracting  chemical  potential  £0(  T),  determined  by 

N,m  /</£/„(£ -£0W(£),  (B6> 

was  used  to  write  the  equation  determining  £  as 

27-Re  2  f  dE N{E)[G{E,io„-,{;)-Go{E,ioni£)]~2  /  dE[f0(E  -£0)-/0(£  — £)]AT(£)  .  (B7> 

»-o 

This  equation  is  solved  iteratively  for  £  at  each  step  of  the  iteration  of  the  system  of  equations. 
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Renormalized  Thermal  Distribution  Function  in  an  Interacting  Electron* Phonon  System 

Warren  E.  Pickett 
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(Received  1  December  1981) 

The  electron-phonon  Interaction  la  used  to  demonstrate  the  Important  effect  of  Inter¬ 
actions  on  the  electronic  distribution  Auction  at  finite  temperature.  It  la  shown  that  the 
usual  picture  of  "thermal  (Fermi)  smearing’’  la  a  greatly  oversimplified  one.  The  dls- 
trltation  function  resulting  from  an  Einstein  spectrum  with  various  couplh*  strwgths  Is 
presented  and  interpreted,  and  an  exact  expression  for  the  spin  susceptibility  la  used  to 
Illustrate  the  utility  of  this  novel  viewpoint  for  thermodynamics. 

PACS  numbers:  65.50. +m,  63.20.Kr,  71.38. +1 


One  of  the  primary  characteristics  of  a  metal 
at  a  finite  temperature  T  is  that  the  crystalline 
electronic  eigenstates  (labeled  by  index  h)  of  en¬ 
ergy  E,  are  occupied  according  to  a  thermal  dis¬ 
tribution  Almost  universally  in  the  theory 

of  metals  this  distribution  is  taken  to  be  that  of 
noninteracting  fermions,  i.e.,  the  Fermi-Dirac 
distribution /„(£,).  Typically  this  assumption  ap¬ 
pears  as  a  thermal  broadening  (“Fermi  smear¬ 
ing"),  given  by  -dy0(E)/dE,  of  some  quantity  over 
a  region  around  the  Fermi  energy  £  F.  The  effect 
of  interactions,  if  included  at  all,  is  not  intro¬ 
duced  into  the  occupation  function.  For  systems 
where  the  density  of  states  (DOS)  function  N(E) 
varies  sufficiently  slowly  near  £  F  thermal  averag¬ 
ing  is  expected  to  be  insensitive  to  the  actual  form 
of/.  For  varying  DOS  systems,  however,  of 
which  the  A 15  class  of  compounds  provides  the 
canonical  example,1  the  understanding  of  their 
anomalous  thermal  behavior  may  depend  critical¬ 
ly  on  the  proper  description  of  /  and  df/dE . 

Two  questions  arise:  What  is  the  effect  of  in¬ 
teractions  on  the  thermal  distribution /,  and  Is 
our  understanding  of  thermodynamic  properties 
clarified  by  a  viewpoint  which  includes  interac¬ 
tions  In/?  In  this  paper  we  use  the  example  of 
the  electron-phonon  (EP)  interaction  to  show  that 
/,  and  thereby  the  interpretation  of  thermodynam¬ 
ic  quantities,  can  be  radically  altered  by  interac¬ 
tions. 

The  distribution  function  /(E  J  is  defined  as  the , 


thermal  expectation  of  the  number  operator  nt . 
For  T  >0,  /  is  given  in  general  in  terms  of  the 
thermodynamic  Green’s  function  G  by  the  rela¬ 
tion2 

/(£,)*T  Tj  G(fe,iw„)expiu)„n,  (1) 

ft  *  •• 

where  *  (2 n  + 1  )jtT  and  n  is  a  positive  infinites¬ 

imal.  (In  our  units  ft  =*fcB  » 1,  and  for  simplicity 
an  isotropic  approximation  for  electrons  will  be 
used.)  It  is  easily  verified  that  (i)  in  the  absence 
of  interactions,  /-/„,  and  (11)  by  converting  the 
sum  to  a  contour  integral  /  can  be  written 

/(E,)-X-dw/o(«)A(*,w)  (2) 

in  terms  of  the  spectral  density  A.  Evidently  / 

*/„  if  and  only  if  A  is  a  6  function  at  <*»■£,.  When 
A  is  broadened  by  interactions,  /  can  differ  con¬ 
siderably  from /0  (a  result  not  often  stated  in 
quantum  statistical  theory  texts),  as  I  now  ex¬ 
plicitly  demonstrate. 

For  simplicity  let  us  initially  consider  a  con¬ 
stant- DOS  electronic  system  interacting  with  an 
Einstein  phonon  spectrum  with  EP  spectral  func¬ 
tion  a*F  given  by 

o2F(w)-(kn/2)6(w-n),  (3) 

where  fl  Is  the  Einstein  frequency  and  X  Is  the  EP 
coupling  constant.  A  straightforward  calculation 
of  the  electronic  self-energy  Z  *Af  -»T  on  the 
real  axis  gives,  with  energies  measured  relative 
to  £  F, 


r(w)  ■  (irAfl/2){/0(Q  -  w) +/0(fl  +  w)+ 2n0(O)}  sgn(w).  (5) 

Here  *  and  n0  denote  the  dlgamma  and  Bose-Elnsteln  functions,  respectively.  The  spectral  density  Is 
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given  simply  by  (£  a  chemical  potential) 


A  (fc,<u)=- 

7 J 


Im 


i 


«-<&, -£)-Z(<*>)| 


(6) 


In  Fig.  1 1  display  M(w)  and  T(u>)  for  several 
values  of  T.  Spectral  functions  at  T  =  0  have  been 
studied  previously  by  Engelsberg  and  Schrieffer3 
and  by  Shimojlma  and  Ichimura,4  who  find  that, 
even  for  |£t  -£  F|  s  0,  spectral  weight  is  spread 
over  a  region  of  several  times  ft  around  E  f .  At 
T  >  0  this  spread  is  further  increased  because  of 
the  increasing  width  T(w)  (Fig.  1).  Neither  of  the 
previous  studies  have  noted  the  affect  of  EP  in¬ 
teractions  on  the  distribution  function. 

The  change  in /  due  to  EP  Interaction  is  repre¬ 
sented  most  easily  in  unitless  differential  form 
-Tdf/dE  ■  The  results,  for  several  values 
of  T  and  X,  are  presented  in  Fig.  2  for  E  >£  *  0. 
[Note  that /'(E)  */'(-£).]  At  T  =  0  Fig.  2(a)  in¬ 
dicates  that  increasing  X  leads  to  the  displace¬ 
ment  of  weight  In  -/'(E)  (i.e.,  occupation  of  bare 
electron  and  hole  states)  from  the  region  \E  -£| 
s  2R  to  the  higher  excitation  energy  tails.  At  T 
-0/2  [Fig.  2(b)]  the  behavior  is  similar.  How¬ 
ever,  as  the  frequency  dependence  of  E  becomes 
sharper  at  T«a  (Fig.  1),  qualitatively  new  be¬ 
havior — negative  weighting  near  £ — can  occur  at 
low  energy,  as  shown  in  Fig.  2(c)  for  T 
This  unusual  behavior  is  exaggerated  by  the  Ein¬ 
stein  spectrum  used  here,  although  similar  be- 


I=M-iT.  Einstein  Spectrum 


FIG.  1.  Electron  self-energy  S  L>) »  Min)  -i  rfa)  for 
an  Einstein  spectrum  at  energy  (I,  for  several  temper¬ 
atures  r.  The  self-energies  are  proportional  to  X  and 
are  shown  for  X*  1. 


havior  may  occur  even  for  a  realistic  spectrum 
at  lower  temperatures.  General  features  of  -/' 
include  these:  (a)  at  high  energy  ,E  -  £|  n  it 
decreases  as  a  Lorentzian  of  width  T  =  (irxn/2)[l 
+  2n0(O)]  rather  than  exponentially  as  does  -//, 
(b)  at  high  T  -/'  is  essentially  Lorentzian 
[T  =  rrAGn0(fi)]  everywhere,  and  (c)  the  behavior 
at  T  =  0  is  as  given  previously  by  Shimojima  and 
Ichimura4:  /  possesses  a  discontinuity  of  (l  +  x)‘l 
at  £ ,  which  indicates  that  the  Fermi  surface  re¬ 
mains  sharp,  and  -/'  has  a  6-function  contribu¬ 
tion  of  corresponding  amplitude,  with  the  re¬ 
maining  weight  X/(l  +  X)  displaced  over  the  range 


E/T 


I  .  I  ^ 

0  2  4  6  8  10 

E/T 


E/T 


FIG.  2.  The  derivative-/'  (£)  ”  -Tdf{E)UE  of  the  ' 
thermal  distribution  function  /  calculated  from  Eqs. 

(2)  and  (6)  and  the  self-energies  of  Fig.  X,  for  several 
values  of  X.  The  dashed  line  shows  -/'  (■  -/#’>  in  the 
absence  of  electron-phonon  interaction  (X»  0).  Note 
that  Interactions  broaden  the  thermal  distribution  con¬ 
siderably  even  for  modest  values  of  X  at  all  tempera¬ 
tures. 
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;£  -£r|  -AH. 

For  systems  in  which  the  electronic  spectrum 
varies  on  the  scale  ol  n,  Z  becomes  a  functional3 
ot  Nt£)  and  the  computation  of  /  becomes  corre¬ 
spondingly  more  involved.  However,  because  of 
the  substantial  widths  in  the  spectral  density 
peaks  for  finite  w  and/or  T,  the  qualitative  be¬ 
havior  of  -/'  will  generally  reflect  that  of  the 
constant- DOS  model.  The  proper  interpretation 
of  thermodynamic  properties ,  however,  can  be 
substantially  altered  in  nonconstant -DOS  systems, 
as  will  now  be  shown. 

The  case  of  the  spin  susceptibility  x,p  will  be 
used  to  demonstrate  that  the  renormalized  ther¬ 
mal  smearing  function  -/' ,  rather  than  -/0' , 
often  arises  naturally  in  the  formalism.  That  this 
should  be  the  case  is  suggested  by  the  relation¬ 
ship  [Eq.  (1)]  between/  and  the  (renormalized) 
Green's  function  G  on  the  one  hand,  and  the  ex- 


99*  Ms  =  UB  LcrfdENGYTLG 

O*  *  3*  *  n 

in  terms  ol  the  number  N „  ot  spin  a  electrons, 
ment)  as 


pansion  of  the  thermodynamic  potential,  which  is 
accomplished  most  concisely*  in  terms  of  G  (rath¬ 
er  than  G0),  on  the  other.  The  conservation  of 
electron  number  Ntl ,  lor  example,  is  given  nat¬ 
urally  as  the  thermal  redistribution  of  crystalline 
states: 


ATel=2/<iE/(E-W<£), 

(7) 

although  it  can  also  be  written 

Mel  =2jdw/0(w-£)R(w) 

(V) 

in  terms  of  a  noninteracting  electron  distribution 
over  a  broadened  and  renormalized  DOS7  3t  given 
by 

&(«  +  £)*  JdE  .V(E  )A  (£  -f ,  u/).  (8) 

This  is  simply  a  question  of  whether  the  w  or  £ 
integral  is  carried  out  first.  However,  from  the 
magnetization  3R,  given  by* 

(9) 

Xip  can  be  written  exactly  (within  this  isotropic  treat- 


„(£,;'«„  )expiw„T7, 


L  c fdEN<£)T£G„2(£  ,iuj  exp m. r,  j  ,.B 

*2u,*fdEN<t)[-  ^~-t)]H>£(£,r)],  (10) 

where  -df/dE  is  identified  by  differentiating  Eq.  (1)  wlthC „*'(£, - (£  -t  -<748tf) -  £„(»«„), 
and  the  thermal  average  X  ot  the  field  derivative  of  the  self-energy  is  defined  by 

K(£  ,T)»- TEC1®  ,i | ITEG*(E  ,iw.  )]*‘ .  (11) 


Equation  (10)  gives  an  interpretation  of  x,p  as  arising  from  the  bare  DOS,  appropriately  enhanced  by 
1 +K  and  averaged  around  t  according  to  the  interacting  thermal  smearing  function.  Equation  (10) 
gives  directly  an  enhancement  due  to  the  EP  interaction  of  the  T  dependence  of  x.p  arising  from  a  peak 
in  N(JE),  as  surmised  by  Bhatt*  for  low  temperature. 

Essentially  all  other  interpretations  of  Xio(T)  have  assumed  a  form  like  Eq.  (10)  with /-/„.  By 
means  of  the  standard  analytic  continuation1  to  express  the  frequency  sum  in  Eq.  (9)  in  terms  ot  an 
integral  over  real  frequencies,  Xip  can  be  written  in  terms  of  /„  as  two  contributions  Xip<0>  +x,p(  °  : 

X,p""-2p,*Jifw(-j^)3l(w>£),  (10') 

which  is  reminiscent  of  but  not  identical  to  the  first  (unenhanced)  term  in  Eq.  (10),  and  the  “enhance¬ 
ment” 


W“  —• y  tojT  dw JdE  jV(£)G*(E  ,w)^|/0(w)  Zf  '  j-  ,  (10”) 

where  Sa'*dZg*/<f(aMi#)U.0  and  ZJ*»Z*/dw  (superscript  ft  denotes  retarded  functions). 
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The  simplicity  of  Eq.  (10),  and  a  general  knowl¬ 
edge  of  the  behavior  of  -a//3£  from  Fig.  2,  al¬ 
lows  one  to  Identify  the  underlying  causes  of  ther¬ 
mal  anomalies  in  exotic  systems.  In  VJC  com¬ 
pounds,  for  example,  the  strong  T  dependence  of 
X  correlates  closely10  with  high  superconducting 
T( ,  and  thus  with  large  X,  exactly  as  Eq.  (10), 
with  a  peak  lnN(E)  near  C  and  broadening  propor¬ 
tional  to  x,  suggests.  It  is  also  evident  that  the 
lncluence  of  the  lattice  should  lead  to  an  Isotope 
effect  on  the  critical  temperature  for  itinerant- 
electron  magnetism  distinct  from  that  proposed 
by  Hopfleld.11  Neither  of  these  properties  Is  evi¬ 
dent  In  the  form  given  In  Eqs.  (10')  and  (10"). 

In  terms  of  the  two  questions  posed  at  the  out¬ 
set,  (1)  the  behavior  of  /  is  qualitatively  as  shown 
in  Fig.  2  and  is  of  itself  useful  and  perhaps  neces¬ 
sary  in  Interpreting  thermodynamic  behavior, 
and  (2)  /  has  been  shown  to  arise  simply  and  nat¬ 
urally  In  the  expression  for  x,p.  In  general,  each 
thermodynamic  quantity  must  be  investigated  In¬ 
dividually  for  a  useful  expression  Involving  /  and/ 
or  fv  It  Is  encouraging  that  Lee  and  Yang11  have 
shown  that  thermodynamics  can  be  formulated 
exactly  in  terms  of  /,  although  the  author  is  un¬ 
aware  of  any  application  of  their  very  formal  ap¬ 
proach  to  metals. 

A  contrast  can  be  drawn  between  the  present 
viewpoint  and  that  of  Fermi-liquid  theory.7  The 
latter  approach  describes  the  low -temperature 
thermodynamic  properties  in  terms  of  noninter¬ 
acting  quasiparticles  described  by /0  and  a  re¬ 
normalized  quasiparticle  density  of  states  X. 
Typically  X  Is  a  constant  for  excitations  of  Inter¬ 
est,  and  this  approach  has  been  very  successful 
for  phenomenological  descriptions.  It  is  pro¬ 
posed  here  that  viewing  the  interactions  as  dis¬ 
tributing  the  excitations  over  the  noninteracting 
spectrum  will  prove  a  more  useful  approach  when 
variation  of  N(E)  la  Important.  This  formulation 


provides  a  conceptual  basis  as  well  as  a  compu¬ 
tational  approach  for  the  detailed  understanding 
of  many  classes  of  interesting  compounds. 
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THEORY  OF  THE  NORMAL  STATE  HEAT 
CAPACITY  OF  NbjSn 

Warren  E.  Pickett  and  Barry  it.  Klein 

Condensed  Matter  Physics  Branch 
Naval  Research  Laboratory 
Washington.  D.C.  20375 

The  experimental  a2F  of  Wolf  et  al  is  applied  to  calculate  the  T-dependence  of  the  electronic  heat  capa¬ 
city.  It  is  found  that,  f  the  phonon  spectrum  is  Debye-like  ( which  is  not  the  case  in  'vb^SnJ.  this  T-dependence 
results  in  an  extrapolated  value  of  \(Ef)  O  +  A)  which  is  overestimated  by  an  amount  which  is  proportional 
N(Ef)k.  equal  to  13%  for  \4jSn.  It  is  further  shown  that  the  change  in  slope  of  Cl  T  found  by  Stewart.  Con 
and  Webb  can  be  modeled  by  a  combined  Debye  and  Einstein  spectrum  with  60  «  267  K  and  having  1.5 %  of 
the  acoustic  modes  at  an  Einstein  frequency  Cl  -  40  K. 

I.  INTRODUCTION 

The  specific  heat  of  NbjSn  has  long  been  an  enigma.  Vieland  and  Wicklund  (1968)  fit  the  normal 
state  specific  heat  C„(  D  above  the  superconducting  transition  temperature  Tc  -  18  K  to  the  form 

CJT-y+BT2. 

However,  the  inferred  normal  state  entropy  S,(Tr>  obtained  by  extrapolating  this  form  below  Tt  was  14% 
larger  than  the  measured  entropy  in  the  superconducting  state  Sr(Tr).  strongly  violating  the  thermo¬ 
dynamic  constraint  S,(TC)  -  Sr(  Tc)  for  a  second  order  phase  transition.  Thus,  although  C.  had  assumed 
the  'low  temperature  form*  of  Eq.  (1)  in  the  range  of  their  measurements,  it  was  evident  that  this  form 
must  change  at  still  lower  temperatures. 

Recently  Stewart,  Con  and  Webb  (1981)  (hereafter  referred  to  as  SCW)  have  resolved  this  particu- 
lar  enigma  by  measuring  Cin  fields  up  to  18  T,  which  lowers  Tr  to  around  6.S  K.  The  Vieland- Wicklund 
form  was  found  to  hold  down  to  12-13  K,  whereupon  the  slope  rapidly  changed  to  give  90  -  207  K  and 
an  extrapolated  y  •  35  ±  3  mJ/mole  K2  at  zero  temperature  compared  to  y  -  96  mj/mole  K2  from 
above  13  K.  Although  this  new  extrapolation  satisfies  the  entropy  constraint,  the  interesting  question 
arises  as  to  what  causes  the  rapid  change  in  slope  ('knee'  or  'kink')  between  10  K  and  12  K.. 

In  this  paper  we  initiate  an  investigation  into  this  unusual  and  unexpected  behavior.  We  assume  at 
the  outset,  following  the  argument  of  Grimvall  (1969).  that  even  in  systems  in  which  the  electron  and 
phonon  subsystems  are  strongly  interacting,  the  specific  heat  nevertheless  can  be  separated  into  two  parts 
interpreted  as  due  to  renormalized  electrons  and  renormalized  phonons.  It  is  not  clear  how  this  approach 
is  to  be  reconciled  with  Allen  and  Hui’s  (1980)  observation  that  the  effects  of  electron-phonon  (EP) 
interaction  appear  to  be  double-counted  in  such  an  approach.  Furthermore.  Grimvall’s  analysis  does  not 
take  anharmonicity  of  the  phonon  system  into  account,  while  the  phonons  are  known  to  show  strong 
temperature  variation  in  NbjSn.  For  the  present  study,  however,  we  shall  proceed  to  investigate  the  heat 
capacity  of  NbjSn  in  terms  of  the  usual  electron/phonon  separation. 

II.  ELECTRONIC  HEAT  CAPACITY 

Grimvall  (1969)  has  shown  how  Prange  and  KadanofTs  (1964)  generalized  quasiparticle  expression 
for  the  electronic  heat  capacity 


C„.,(  T)  -  lN(Ef)  f  dE  E  |l  -  T)  illi 


df<E)  3 ,tf(f.  T)  8/(£) 
3  T  3  T  dE 


can  be  derived  from  field  thee-  cal  considerations.  Here  \t  is  the  real  part  of  the  electronic  self-energy, 
/is  the  Fermi  function  and  other  quantities  have  their  usual  meanings.  Using  standard  expressions  for 
M  this  can  be  written  explicitly  as 

c,,!T)  -  (y0  +  yi< nir  s  * ( nr.  (2) 


y0-  (2n2/3)  V<£,  I  kf 


3  r  ...  r  .r  dt'IE)  3 flE)  r  ,  ,  ,  IE-E)2 

y'(T)/yn  T^TFp-  J  d£J  d£  —  —J  dua  F(“'  TTeTZ-  ,J 

where  a2Fis  the  electron-phonon  spectral  function. 

? u ;•? i rcc r. t* cti t  i ; y  In  <J-  f-u.vid  Ratals  1932 
-  1C3"  : 1 ‘ r."  : r :  i rj:..  Xarlcfuh* 
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Equation  (4)  has  been  evaluated  previously  for  an  Einstein  model  (Grimvall  1968).  for  Pb  and  Hg 
(Grimvall  1969)  and  for  amorphous  Pb  <Bergmann  era/  1971)  Here  we  use  the  or/" function  for  NbiSn 
obtained  by  the  inversion  of  tunneling  data  by  Wolf  et  at  <1980)  to  evaluate  Eq.  <4).  In  Fig.  I  we  show 
the  result  for  CnJyoT  -  1  +  y^  Thy 0  in  the  range  0  to  100  K  The  result  is  similar  io  those  of  Pb  and 
Hg  (scaled  to  the  suffer  phonon  spectrum  of  NbiSn)  At  9  K  y,  peaks  at  a  value  10%  higher  than  us 
value  Ay0  (A  -  1.78)  at  T  -  0.  Then  y,/y0  decreases  rapidly,  reaching  us  minimum  value  of  -0  24  at 
approximately  100  K. 

This  drop  in  C„,  below  yoThas  interesting  implications.  The  thermal  mass  enhancement  y,  is  pro¬ 
portional  to  the  magnitude  of  a:F.  i.e..  A,  and  its  shape  depends  only  on  the  shape  of  a:F  As  A  is 
increased.  C„,  will  vanish  at  the  temperature  at  which  y,  has  its  minimum  for  a  critical  value  a,  of  A. 
and  will  become  negative  for  a  range  of  temperatures  when  a  >  A,.  For  a:F  having  the  NbiSn  shape 
Ar  -  7.4.  For  the  Pb  shape  Af  -  8.2  and  the  minimum  in  y.  occurs  near  35  K.  for  an  Einstein  spectrum 
(which  gives  the  sharpest  structure  possible)  Ar  -  3. 7  and  the  minimum  occurs  near  A,r  -  7iflEl„„t,„,2 
This  negative  ’electronic"  specific  heat  causes  no  obvious  violation  of  thermodynamic  principles,  since 
increasing  a  will  lead  to  a  softening  of  the  phonon  spectrum  and  thus  an  increase  in  the  ’lattice’  contri¬ 
bution.  which  should  keep  the  total  heat  rapacity  positive.  (The  condition  C  >  0  can  be  used  to  put  a 
weak  constraint  on  how  hard  the  phonon  spectrum  can  be.)  Nevertheless,  this  behavior  raises  further 
questions  into  the  separation  of  C  into  'electronic*  and  ’lattice*  contributions  in  strongly  interacting  EP 
systems. 

It  is  now  widely  accepted  that  variation  of  .V<£>  on  the  scale  of  phonon  energies  is  a  likely 
occurrence  in  high  T  4  15  compounds.  NbiSn  in  particular  In  such  a  case  the  expression  for  C.r  is 
altered  in  several  wavs.  Following  Grimvall’s  <1969)  derivation  of  Eq.  <1).  variation  of  Vi£)  near  £f 
results  in  an  explicit  dependence  of  C.,  on  the  imaginary  part  T  of  the  electronic  self-energy  (also  see 
Gnmvall  1978);  i.e..  a  breakdown  of  the  generalized  quasiparticle  expression  Eq.  (I)  This  dependence 
will  be  neglected  for  reasons  to  be  discussed  below  Then  Eq.  <1)  is  modified  by  replacing  .V(£*)  by 
V(£)  under  the  integral,  with  a  similar  replacement  in  the  equation  defining  A HE.  T)  (see  Grimvall 
1969).  Finally  one  obtains 

7=6?  *SF  ‘Sf1  /  *  *■"- 

with  the  normalization  yi<0)  -  A.  Equation  (5)  has  been  evaluated  by  Fradin  (1975.  1977)  for  a  model 
AM  £)  function  thought  to  be  appropriate  for  V,Ga.  but  no  comparison  with  experiment  was  attempted 


Nb3Sn 


Nb,Sn  X=1.78 


T  (K) 

ENERGY  (mRy) 

Fig.  1  —  The  variation  o/y/y„  »  1  r  y,/yj  with  tempera¬ 
ture  in  VOj  Sn.  The  difference  between  this  curve  and  unity 
is  proportional  to  A . 

Fix.  2  —  Density  of  states  near 
£,  -  0  0/  VA,Sn  (solid  line)  Irom 
Klein  et  al  (IVtXI.  Also  shown  are 
the  effects  of  delect  scanennx  due  to 
5.  ID  20  and  All  *tfl  cm  of  residual 
resistivity.  Vertical  lines  denote  the 
respective  Fermi  enerytes.  The  III 
1*11  cm  curve  is  used  lor  the  tabula¬ 
tions  described  in  the  text . 
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We  have  evaluated  Eq.  <5)  for  NbiSn  using  the  experimental  data  for  a:F and  the  V(£)  function 
calculated  by  Klein  et  al  ( 1 978).  Since  NbjSn  samples  invariably  have  a  residual  resistivity  of  10  itSl  cm 
(or  morel,  we  have  broadened  .V(£)  accordingly  (Pickett  1982)  These  V(£)  functions  and  the  shift  in 
E,  due  to  smearing  are  pictured  in  Fig.  2.  It  should  be  noted  that  ,V<£)  functions  with  even  sharper 
structure  than  that  calculated  by  Klein  et  al  will  not  lead  to  broadened  V<£)  functions  with  significantly 
more  structure  than  the  '10  cm'  curve  in  Fig.  1. 

In  the  range  of  interest  here  IT  ^  30  K)  the  temperature  variation  of  the  chemical  potential  is 
negligible.  In  this  range  yi  from  Eq.  (5)  differs  from  that  of  Eq.  (4)  by  less  than  0.2%.  This  null  correc¬ 
tion  results  because  ,V<£>  is  essentially  linear  in  the  energy  region  sampled  by  Eq.  (S)  at  these  tempera¬ 
tures.  and  increased  weighting  at  £  <  Ef  is  cancelled  by  decreased  weight  at  £  >  £f.  This  is  also  why 
the  contribution  to  C,,  from  T(£.  T)  (mentioned  above)  is  negligible.  This  contribution  arises  only 
from  changes  due  to  the  energy  variation  of  V(£).  and  for  a  linear  variation  these  changes  are  vanishingly 
small.  Even  if  EF  is  artificallv  placed  at  the  peak  of  the  "10  nfl  cm'  curve  in  Fig.  2,  where  the  effect  is 
largest,  the  correction  is  only  2%  at  30  K. 

It  remains  to  determine  the  effect  of  the  variation  of  y ,  upon  the  interpretation  of  C».  In  Fig.  3 
yi/-)  -  C„JT is  plotted  versus  T1  along  with  the  normal  state  experimental  data  of  SCW.  On  the  scale 
shown  the  peak  in  y(D  is  comparitiveiy  weak  and  has  little  to  do  with  the  knee  in  the  data,  despite  the 
fact  that  both  structures  occur  in  the  same  temperature  range.  However,  above  12  K  the  curve  is 
approximately  linear,  with  slope  $  -  -0.014  mi/mole  K4  Removing  this  linear  contribution  from  the 
linear  fit  to  the  data  above  12  K  results  in  80  -  2b  7  K  in  this  region,  rather  than  the  SCW  value  of  270 
K. 

The  temperature  variation  of  y,  can  lead  to  a  more  significant  correction  to  the  extrapolated  value 
of  y0(l  +  A).  Assuming  the  lattice  contribution  C.,  ~  @P  (i.e.,  extrapolates  to  zero),  the  approxi¬ 
mately  linear  decrease  in  y0  +  y,  vs  P  above  the  peak  in  yt  leads  to  an  overestimate  of 
y (0)  —  ,V(£f)  (I  +  x)  which  is  proportional  to  .V(£f)A.  For  NbjSn  this  overestimate  is  13%.  If  SIEf) 
is  known  independently,  the  overestimate  of  A  is  further  leveraged  by  (I  +  A)/A,  equal  to  20%  for 
NbjSn.  Use  of  incorrect  values  of  V(£f)  can  of  course  lead  to  worse  empirical  estimates  of  A.  This 
overestimate  of  A  (0)  in  itself  can  lead  to  an  apparent  violation  of  the  entropy  constraint  mentioned  in  the 
Introduction.  For  the  actual  data  for  NbjSn  above  12  K,  where  the  phonon  contribution  does  not  extra¬ 
polate  to  zero,  the  overestimate  of  y0(l  +  A)  is  4.7  mj/mole  K2,  a  rather  small  fraction  of  the  difference 
(61  mJ/moie  K2)  between  the  high  and  low  temperature  extrapolations  of  SCW. 
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Fig.  3  —  C„  ,/  T  vs  P  from  0  to  30  K.  The  calculated 
curve  is  normalized  to  35  mJ/mole  K1  at  T  —  0.  Data 
points  (x.  H  —  0;  +  .  7  T:  o,  18  T.  A.  19  T)  are 
normal  state  values  at  various  magnetic  fields  (Stewart 
et  al  1981).  The  lowest  4-7  data  points  reflect  a 
broadened  transition  at  T,  =  6  K  in  18  T.  The  doited 
line  represents  the  extrapolation  of  Stewart  et  al 
(1981). 


Fig.  4  —  (formal  state  data  points  of  Stewart  et  al 
(1981).  with  constant  electronic  contribution  removed, 
compared  with  the  Debye  +  Einstein  model  (plotted 
curves,  see  text). 
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III.  LATTICE  HEAT  CAPACITY 

The  complexity  of  the  lattice  dynamics  of  NbjSn  makes  a  quantitative  evaluation  or  the  lattice  heal 
capacity  C,l  7")  impossible.  "Good"  samples— typically  ihose  with  low  residual  resistivity— undergo  a 
cubic-to-tetragonal  transition  at  T„  —  43  -  5!  K.  as  did  the  SCW  sample.  Definitive  neutron  scattering 
studies  iShirane  and  Axe  1971,  1 978:  Axe  and  Shtrane  1973a.b)  have  uncovered  a  "soft  mode"  type  of 
behavior  of  the  <{.{.0>  7"i  branch  as  T approaches  Tm  from  above,  the  small  J  modes  tend  to  disappear 
under  a  narrow  central  peak  broadened  by  instrumental  resolution,  whereas  for  T  <  Tm  only  broad 
scattering  centered  at  zero  energy  transfer  could  be  observed  for  ;  $  0.2  (ir la).  Thus  the  lattice  dynam¬ 
ics  are  not  only  unknown  but  are  also  expected  to  be  strongly  anharmomc.  at  least  near  Tm. 

We  initially  investigated  two  (spherically  symmetric)  acoustic  phonon  dispersion  models  with 
upward  "kinks*  in  one  TA  branch,  as  found  by  Axe  and  Shirane  11973b)  for  T  >  Tm.  In  one  of  the 
models  the  branches  were  tied  to  the  T-dependent  sound  velocities  at  {  —  0.  Both  models  displayed  a 
smoothly  increasing  Debye  temperature  with  temperature,  qualitatively  similar  to  the  data,  but  neither 
could  be  fit  even  semiquantitatively  to  the  data. 

ft  is  therefore  surprising  that  the  following,  much  simpler,  model  provides  a  quantitative  fit  to  the 
SCW  data  above  8  K.  We  take  a  fraction  6 it  of  the  acoustic  modes  to  be  modeled  by  an  Einstein  fre¬ 
quency  fl  and  assume  the  remaining  1  -  in  to  be  characterized  by  a  Debye  temperature  9.  The  results 
for  0  -  40  K  sm  3.5  meV,  6 n  -  0.015.  9  -  267  K  are  shown  in  Fig.  4  compared  to  the  normal  state  data 
of  SCW  (For  this  comparison  we  have  neglected  the  small  T-dependence  of  y  and  subtracted  out  a  con¬ 
stant  electronic  contribution  of  35  mJ/mole  K2)  The  model  provides  an  excellent  fit  above  T1  -  60  K2. 
with  a  possible  small  discrepancy  in  the  range  100  <  T2  $  150  K2  where  the  data  points  are  sparse. 
SCW  concluded  that  an  'abrupt  change"  occurs  in  this  range,  but  a  few  more  data  points  are  needed  to 
ascertain  the  abruptness.  We  emphasize  that  this  model  of  the  phonon  spectrum  should  not  be  inter¬ 
preted  literally;  it  does  not.  for  example,  reproduce  f)0  -  208  K  in  the  range  8  <  T2  ^  20  as  is  found 
from  data  taken  in  the  superconducting  state  (Stewart  et  at  1981).  Nevertheless  it  removes  much  of  the 
mystery  in  the  NbjSn  data  by  showing  how  an  "excess"  of  soft  modes  can  produce  the  observed  behavior. 

IV.  final  points 

A  question  remains:  why,  SCW  have  asked,  is  the  abrupt  change  in  slope  not  observed  in  the 
superconducting  state  data?  Again,  the  SCW  data  in  the  superconducting  state  are  sparse  in  this  range 
and  it  is  the  authors'  opinion  that  structure  in  this  region,  where  the  electronic  contribution  is  increasing 
rapidly,  should  not  be  ruled  out.  A  definitive  answer  to  this  question  will  require  a  calculation  of  the 
electronic  heat  capacity  in  the  superconducting  state,  as  was  done  by  Daams  and  Carbotte  (1979)  using 
earlier  experimental  data,  in  addition  to  more  experimental  data. 

Finally  we  note  that  abrupt  changes  of  slope  of  the  C/T  curve  at  low  temperature  are  not 
uncommon:  Nb  shows  such  anomalies  at  3  K  and  9.5  K  (see  Leupold  er  al  1977  and  references  therein). 
NbjSn  is  however  unusual  in  the  mattnitude  of  the  slope  change. 

We  are  indebted  to  G  R.  Stewart  for  helpful  conversations  and  for  providing  the  experimental  data 
in  numerical  form. 
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